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FINAL TECHNICAL REPORT 


The NASA Solar- Terrestrial-Theory program, “Research in Solar Plasma Theory” has 
been supported by Grant NAGW-93 at the University of California, Irvine since 1 June 1980. 
During the intervening years, the UCI group has produced 58 refereed publications [listed in 
Appendix I] and three Ph.D. dissertations (#2, #12, #43) acknowledging this grant. The 
results of this research effort have been described in nearly a score of Invited Talks presented 
at scientific meetings and workshops. In addition, we have created five numerical codes for 
the supercomputer simulation of solar activity phenomena, and have significantly adapted 
three others. [Brief descriptions of those codes are given in Appendix II.] 

Since the Abstracts of all of our publications are provided in Appendix III, we will 
concentrate in this report on the main thrusts and significance of our research results. [The 
numbers refer to the publication list provided in Appendix I.] Three recent invited reviews 
of the UCI solar-activity simulation program (#42, #50, #56) are also available in print. 

A. CORONAL STRUCTURE AND DYNAMICS 

We have spent considerable effort in trying to understand the magnetohydrodynamic (MHD) 
structure, evolution and stability of the myriad of coronal forms observed by space-borne 
telescopes. 

A.l. Magnetic Loops 

An early emphasis was to understand the MHD stability of the ubiquitous coronal loops 
so comprehensively observed by the wide spectrum of Skylab instruments. The fact that 
they persisted for days flew in the face of the well-known laboratory instability of localized, 
magnetically confined, dense plasma structures. Earlier work, by the UCI group and others, 
had suggested a qualitative stabilizing influence, magnetic line tying at the loop ends. We 
then set out to make the effects of line tying more realistic physically (#1) and more rigorous 
mathematically (#3). The end result was the ability to compute the instability boundaries, 
in terms of the length-to-radius aspect ratio for realistic coronal loops, and finally to under- 
stand the reason for their MHD stability (#4). A second aspect of loop behavior is their 
magnetic activity, when they make the transition to the resistive- reconnection process of a 
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compact flare. Both the linear and nonlinear (#20) energy- release aspects of such loop flares 
were evaluated. 

We have also recently developed and applied a three-dimensional MHD model to the 
study of the dynamical evolution of coronal loops. This model can self-consistently represent 
the various phases of development of coronal loops, from the quasi-static evolution though 
equilibrium states as a result of twisting photospheric flow profiles, through the linear ideal 
and resistive MHD instability phase, and through the nonlinear behavior after kink insta- 
bilities grow to large amplitudes and cause field-line reconnection. We have used this model 
to study the equilibrium properties and ideal linear MHD kink stability of coronal loops 
(#52). We have found that typical localized-twist profiles produce loops with zero net axial 
current which are unstable to internal kink modes (and not to the faster-growing external 
kink modes). We have also confirmed and extended previous results on the stability of the 
idealized Gold-Hoyle model of a loop (.with uniform twist). A preliminary solution of the 
MHD equations for the nonlinear evolution of the internal kink mode has shown that concen- 
trated current layers develop as a result of the instability (#50). These layers have profound 
consequences in regard to the formation of current sheets in the line-tied solar corona, and 
are intimately related to the coronal heating problem. 

A. 2. Helmet Streamers 

In connection with our mass-ejection studies, to be described in the next section, we 
have investigated helmet streamers, large-scale structures that are seen frequently in the 
solar corona. Observations indicate that coronal mass ejections (CMEs) are often associated 
with helmet streamers and the disruption of overlying streamers by CMEs has been reported. 
These observational studies have convinced us that investigating the stability of streamer-like 
magnetic configurations may yield important insights about CME initiation. 

Previous simulation studies of streamers were two-dimensional (axisymmetry assumed). 
Although a two-dimensional description is useful, observations indicate that coronal stream- 
ers have a finite longitudinal extent, and are thus three-dimensional objects. The 3-D struc- 
ture of streamers could be an important aspect of their stability. Using a transonic Parker 
wind as the initial condition for the density, pressure and velocity (as was done in previ- 
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ous studies) and a non-axisymmetric potential magnetic field (never previously attempted) 
we have performed the first three-dimensional simulations of coronal streamers (#58). A 
streamerlike configuration with a finite longitudinal extent forms. The closed field region 
has a higher plasma density and is bounded by a current sheet. We have also treated open 
and closed fields in the stellar context (#6). 

A. 3. Mass Ejections 

An important manifestation of solar activity is the coronal mass ejection. CMEs add 
mass and magnetic fields to the solar wind and may cause interplanetary shock waves, which 
provide a causal link between solar and geomagnetic activity. We have performed several 
studies of the evolution of CMEs in the solar corona. As most previous simulations used a 
thermal driver as the mechanism for CME initiation, we have examined some of the possible 
effects of a non-thermal driver; specifically, plasmoid ejection. As a simple model of an 
erupting prominence we accelerated a cold, dense, plasma parcel to a given outward velocity, 
then followed the time evolution as the parcel moved through an initially hydrostatic, current- 
free atmosphere. We found that this simple driving mechanism yielded results that are 
consistent with many of the observational characteristics of CMEs, suggesting that the form 
of the driving mechanism in fact plays an influential role in the evolution of CMEs (#48). 
We have also performed studies in which we incorporated a plasmoid-like magnetic topology 
and allowed the force resulting from the current associated with this magnetic structure to 
drive plasma outward (#54). 

B. CORONAL HEATING 

One of the principal unsolved problems of solar physics is the mechanism for heating the 
corona. We have made significant contributions to the elucidation of the Parker “topological 
dissipation” proposal and of the Alfven surface-wave model. 

B. 1. Current Filaments 

We have studied an atmospheric heating model, due to Parker (1972, 1986), in which 
the convective twisting of the coronal magnetic field leads to strong currents. We developed 
detailed solutions to the 3-D ideal-MHD equations describing the coronal response when a 
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sequence of smooth, quasi-random, photospheric flows was applied (#40). These solutions 
show that smooth flows characteristic of photospheric convection can induce a fine-scale 
spatial structure in the coronal magnetic field, leading to the formation of concentrated 
current filaments. The resistive dissipation of these currents can heat the corona (Parker 
1986, van Ballegooijen 1986, 1988). 

The principal results show that smooth flows can be very effective in generating current 
filaments in the corona. In particular, spatial structure is induced at increasingly shorter 
length scales, with an exponential dependence of the shortest length scale on the magnetic 
foot-point displacement. The spatial structure is generated through nonlinearities in the 
MHD equations, and cascades from the driven long wavelengths to shorter wavelengths. 
This mechanism could be an ingredient of a coronal heating model if the cascade of spatial 
structure continues to the short length scales at which the magnetic field can be dissipated 
resistively or by reconnection in the highly conducting corona. In future research it would 
be desirable to determine whether this mechanism satisfies the quantitative requirements of 
a coronal heating model. 

B. 2. Alfven Surface Waves 

We have investigated the coronal-heating potential arising from the propagation and 
dissipation of Alfven waves in a non-uniform plasma embedded in a magnetic field with 
continuous transverse gradients (#27). Such configurations arise in the neighborhood of 
active regions on the solar surface and (with a larger scale and at higher altitude) at the 
edges of coronal holes. The purpose of this study was to understand the energy transport 
to the corona from the photosphere, where the convective motions of the magnetized fluid 
are believed to be the source for the coronal energy input. Although the dissipation of 
conventional, uniform- medium, Alfven waves is insufficient to deposit energy in the coronal 
plasma, we have shown that there exists a normal mode which has the form of a surface 
wave in this continuous medium, that has a damping rate proportional to the magnetic 
field gradient. In the lowest-order approximation, the dissipation is also proportional to the 
square of the wavenumber and independent of resistivity as long as the Lundquist number S 
is sufficiently large. The damping of these waves provides a potentially important mechanism 
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to transport energy from the photosphere to the corona where it can be effectively dissipated. 

To investigate the properties of these waves in the lower atmosphere, where the magnetic 
Reynolds number is not as large as in the corona, we carried out the second-order calculation 
of their damping rates and found that it was proportional to S ~ 0 33 (#29) where S is the 
Lundquist number. Similar results were also obtained in numerical simulations by studying 
the space-time evolution of the wave energy (#19) and the phase-mixing of these propagating 
normal modes (#24), including the effects of viscosity (#30). 

The results of these studies were then applied to evaluate the energy-deposition rate in 
the solar atmosphere (#35). In active regions, with strong magnetic fields and mostly closed 
field lines, these waves suffer fast damping as they propagate upward because of their short 
wavelengths, as limited by the length of the field lines, and because of the small horizontal 
scale width of the active-region fields. Their contribution to coronal heating is furthered by 
the fact that they are trapped in the closed-field regions and easily dissipated within a few 
reflections. In the open-field regions of coronal holes, their dissipation length is much longer. 
Although they deposit less energy in the corona in this case, much of their energy can be 
absorbed into the solar-wind plasma flowing out of these magnetically open regions. 

B. 3. Anomalous Resistivity 

The heating effects of anomalous resistivity and wave dissipation in coronal arches have 
been explored. For steady currents, exciting turbulence is not easy — electron drift velocities 
must exceed 0.1 v te to excite plasma instabilities other than drift waves. Above this threshold, 
the current-driven ion-cyclotron instability seems the best candidate for a resistive heating 
mechanism, since it is easy to excite and difficult to suppress. Detailed numerical calculations 
show that currents can heat arches to ~100 eV. Convective losses then cause the instability to 
“percolate” at marginal stability, maintaining the temperature(#8). Ion-cyclotron waves are 
easier to excite than the ion-acoustic instability and can also transport heat across magnetic 
field lines. A current-carrying core of size a ~ 100 m., set by the scale of self-field gradients 
in a loop, can heat a warm cocoon. Within the cocoon, temperature declines on a scale of 
kilometers, set by energetics. Electron Landau damping of the waves deposits energy in the 
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cocoon(#41). The damping electrons couple through convective cooling to the foot-point 
boundaries. Complications such as magnetic shear can affect these estimates, but cross-field 
wave transport promises a natural way to heat large volumes with small current threads. 

C. FILAMENT FORMATION 

Solar filaments and prominences are believed to form as a result of runaway radiation in a 
magnetically insulated region (Field 1965). Both because of their intrinsic thermodynamic 
interest and because of their function as the most reliable flare precursor (#21, #32), we 
have been engaged in a long-term study of the coupled energetics and dynamics of plasmas 
immersed in sheared magnetic fields. Our work has emphasized parameter regimes applicable 
to the formation of cool and dense structures in the magnetized solar atmosphere such as 
prominences. 

C. 1. Thermal Instability 

In order to clarify the growth and structure of sheared-field condensations, which were 
originally found numerically (Chiuderi and Van Hoven 1979), we first addressed the radiative 
instability in two dimensions (#28). We demonstrated how plasma dynamics can generate 
linear modes localized within a very narrow region of the magnetic shear layer, provided 
merely that the magnitude of the temperature perturbation is constrained to decay in regions 
where the wave vector is parallel to the equilibrium field. In practice this constraint is 
provided by the temperature-leveling effects of parallel thermal conduction, but the mode 
width is not explicitly determined by the conduction coefficient. This paper (#28) also 
defined the accessible regions of growth- rate/ wavenumber space, and reported the results of 
a series of numerical calculations that characterize the mode structure of the long- wavelength, 
sheared-field condensations first described by Chiuderi and Van Hoven (1979). 

Our next effort was devoted to an extension of the completeness of the energy- transport 
model applicable to thermal instability. We first introduced the effects of finite perpendicular 
(to B) thermal conductivity k_l, which bring new unstable modes to the plasma spectrum in 
a way similar to the added dynamic effects of resistivity. In an analytic study (#36) which 
complemented a parallel computational effort, we detailed the dynamics and structure of 


6 


these modes, and obtained their growth rates for arbitrary electron-ion thermal coupling 
(#31). A crucial result is that these modes grow preferentially at an intermediate position 
in the magnetic shear layer where the parallel sound-speed (mass-flow) rate is dominant 
(#53). Thus the radiative cooling occurs under isobaric conditions; that is, with a large 
accompanying density enhancement which drives the cooling at a nonlinearly accelerating 
rate (see, e.g., #47). 

Any dynamic phenomenon which is observable on the sun must represent a large-amplitude 
state, so that one must provide a fully nonlinear treatment. We have attacked the filament- 
formation problem from this perspective, which required the development of a new 2.5- 
dimensional simulation code. Our first results described the accelerating cooling of the long- 
wavelength radiative modes, which eventually drive the local temperature down at 10 3 x the 
initial rate (#34). 

During this period we discovered the crucial dynamic effects of the short- wavelength #Cx 
modes described above. The fact that the sound-speed rate is the fastest one in the problem 
means that a significant parallel (to B) mass flow occurs. To demonstrate this effect, we 
then performed a comprehensive set of supercomputer experiments (#47, #53) with various 
excitations of these modes, including the empirically relevant case of a noise- spectrum input. 

We were able to show, in this series of simulations, that an even more dynamic effect 
arises nonlinearly (#47). After an early phase of isobaric cooling, the local pressure begins to 
drop, leading to a strong, parallel, mass inflow to the condensation. The final state produced 
in these simulations (resulting from a variety of different initial states) provides a very good 
match to the observed properties of a coronal prominence/filament. The density increases 
by 10 2 x from coronal values, while the temperature drops to 10 -2 x the initial value. In 
addition, in a new and unique result, the angle between the filament axis and the local 
magnetic field agrees with observations (#47). These results indicate that we now know 
how filaments form, dynamically and energetically, in the sheared field of the ambient solar 
corona. 
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C. 2. Prominence Condensation 

We have proposed and demonstrated a new mechanism that triggers the formation of a 
coronal prominence in an otherwise thermally stable magnetic arcade. Based on the concep- 
tion that heat in the corona is supplied from the photosphere through various mechanisms, 
such as Alfven waves, and is guided by the magnetic field, we argue that, if the energy supply 
to the apex is curtailed, thermal conduction may not be able to sustain the apex temperature, 
resulting in the initiation of local cooling. By numerical simulation, we have demonstrated 
that a prominence can indeed form naturally in this model (#51) which exhibits strong 
siphon flow from the chromosphere as needed to obtain sufficient mass in the condensation. 
This multi-phase process results from a single causative input, in contrast with other previ- 
ous models which require multi-step events in the corona (Poland and Mariska 1986). We 
further simulated the effect of the bending of the magnetic field due to the heavy condensed 
materials at the apex, where the field lines are essentially horizontal before the condensation, 
by modifying the gravity projected along the field. A magnetic-field “hammock” is formed, 
as a result of the modified geometry, which stably traps the condensed plasma in the form 
of a levitated apex prominence of the Kippenhahn-Schlviter (1957) type. 

D. FLARE ENERGY RELEASE 

The energy which drives a solar flare is believed to be stored in the local, highly stressed, 
active-region, magnetic field. In order to release this energy, one must change the topology 
of (i.e., “reconnect”) this field, using a resistive process. The only known truly dynamic 
mechanism to accomplish this is by an instability called the tearing mode which has been a 
long-time focus of this program. 

D. 1. Magnetic Reconnection 

Among our early efforts in evaluating the application of resistive tearing to the flare- 
release mechanism was a quantitative study of the high-conductivity (high Lundquist number 
S) growth and scaling of this instability (#5). 

We followed this linear study with a series of nonlinear simulations of the energy output of 
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the resistive tearing instability of a plane field-shear layer. These computations, which held 
the high-5 record (at 10 6 ) for a number of years, showed that more than 30% of the stored 
magnetic energy could be liberated (#15, #18, #26). Comparable results were obtained in 
a different set of simulations performed with a cylindrical-pinch (coronal loop) field (#20). 

A second nonlinear focus of our flare/reconnection studies involves the possible speed- 
up effects of MHD turbulence , whether through modifications of the resistivity or of the 
dynamics. Our first effort diagnosed the temporal evolution of a (time-averaged) anomalous 
resistivity proposed by Biskamp and Welter (1983). We demonstrated (#38) that the sign 
and magnitude of the “resistivity”, claimed to depend on the turbulent average values of 
(kinetic energy)— (magnetic energy), quickly oscilllated and decayed without a significant 
effect on the tearing mode. 

An alternative proposed mechanism for the modification of magnetic reconnection in- 
volves the effect of co-existing MHD fluctuations. Matthaeus and Lamkin (1986) have 
claimed that turbulence can drive reconnection, and we have completed an investigation 
of this question. Our method was to apply varying levels of wide-spectrum noise input to a 
tearing-mode-unstable magnetic field, similar to that used by Matthaeus and Lamkin, but to 
follow the temporal evolution of the system for much longer times. Using either our periodic 
or our ADI 2.5-D codes (see Appendix II), which have both been previously used to study 
nonlinear tearing, we can compute for periods commensurate with the saturation time. We 
find that, when starting with low levels of noise, the expected fastest- growing tearing-mode 
excitation appears, and grows to its usual nonlinear saturation or stored-energy level. As 
the turbulence amplitude is progressively increased to the natural stored magnetic-energy 
level, near that chosen by Mattheus and Lamkin as an initial excitation, the tearing mode 
is obscured and modified (buffeted about, actually), but the degree of energy release is un- 
changed. The results of these simulations are reported in Deeds and Van Hoven (#43). 
These two simulations formed the core of the Ph. D. dissertation of D. Deeds (#45). 

D. 2. Radiative Effects 

In parallel with these nonlinear studies, we began an investigation of the effects on resis- 
tive tearing of the addition of those energy-transport contributions which are important in 
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the solar application; namely, ohmic heating, thermal conduction and (predominantly) radi- 
ation loss (#16). In particular, this effort produced a very complete (and clear) treatment 
of this complex of effects in an analytic boundary-layer calculation (#10). 

We then expanded this study, motivated by the close spatial and temporal relation be- 
tween filaments and flares, to a series of studies on the interaction between radiative and 
reconnection instabilities (#7). Computational investigations of this coupling, which arises 
from the temperature dependence of collisional resistivity, produced a number of publications 
detailing the growth rates and mode structures (#7, #11, #14, #25), and the additional 
effects of compressibility and viscosity (#22). These studies also resulted in the completion 
of a Ph.D. dissertation for T. Tachi (#12). 

The crowning result of this series of energy-transport modelling efforts was the discovery 
and elucidation of the radiative-tearing instability , a very fast magnetic-reconnection mode 
(#13, #17). This instability’s high speed arises from the effect of radiative cooling on the 
Coulomb resistivity (77 oc T~?) which produces a very effective resistive-diffusion reconnec- 
tion layer. 

D. 3. Observable Diagnostics 

As a result of the physical complexity of the solar flare event, we have early-on sought 
to characterize and learn from the relative simplicity of the preflare state. To this end the 
PI served as a group leader in the COSPAR sponsored Flare Build- up Study, devoted to 
precursors and onset (#21, #32), with particular attention to the connection (both spa- 
tial/magnetic and temporal) with pre-existing filaments. A second diagnostic effort sought 
to evaluate the prospects for flare x-ray polarimetry to distinguish between thermal and 
non-thermal models of the impulsive phase (#33, #39). 

We have previously investigated a loop-flare model based on a resistive-pinch configu- 
ration with a force-free magnetic field (Mok and Van Hoven 1982). Under the line-tied 
boundary conditions presented to the loop field by the dense photosphere, the fundamen- 
tal unstable mode is localized at the apex of the loop where the local magnetic field is 
nearly azimuthal. The induced electric (DC) field, as a result of magnetic reconnection, 
can accelerate charged particles in a very efficient manner along the local, closed, field lines 
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like a synchrotron (#9). This picture is consistent with observations indicating a compact, 
partially polarized, microwave source at the loop apex during a flare. 

All reconnection flare models naturally lead to thick-target x radiation as the energetic 
electrons drift across the field lines and are directed down toward the photosphere. Since 
electron beams going through a plasma are often unstable and thus subject to scattering 
by a variety of plasma oscillations, their transmission through the solar atmosphere is an 
important question. Based on their energy spectrum, deduced from hard x rays, we have 
studied a potentially unstable plasma wave (#23) driven by these fast electrons through the 
anomalous Doppler resonance (u> — k • v — fi e = 0). The instability condition was shown to 
be u> pe /O e > 1.9 for a power-law distribution with a wide range of hardnesses, where u pe is 
the local plasma frequency and fl e is the electron gyrofrequency. This instability threshold is 
often crossed when the fast electrons pass the transition region where the local density and 
plasma frequency rise. The slowing down of the beam, as a result of pitch- angle scattering, 
can cause intense heating at a location different from the one predicted by conventional 
models (#23). There are also a few unexplained aspects of the UV observations, which may 
be attributed to this instability. 

Another flare-diagnostic study (#55) explored a model for impulsive energization of the 
arch foot points by bombarding them with slugs of magnetically-driven plasma. These 
plasmoids are released when magnetic field lines reconnect at the arch peak. The plasma in 
the resulting magnetic pockets then accelerate as the field lines retract, slamming plasma 
into the dense footpoint layers and heating them rapidly. This can account for impulsive 
x-ray emission. 

E. SPECIAL TOPICS 

An early effort in this solar plasma theory program was devoted to the evaluation of the influ- 
ence on radiative opacity of plasma shielding in the solar interior. Some effects which could 
reduce the central temperature of the sun, and thereby lessen the neutrino-flux discrepancy, 
were described in the Ph.D. dissertation of R. 0. Hunter (#2). 

At the other extreme of the heliosphere, we have performed two studies of plasma flow 
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past Jupiter’s satellite Io. In the first, we verified that 3-D MHD simulations could exhibit 
the formation of “ Alfven wings” , a region of standing Alfven-wave perturbations generated 
by the interaction of the Jovian magnetic field with Io’s ionosphere. We also discovered that 
standing slow-mode perturbations were present, a feature that had not been quantitatively 
modeled previously(#37). 

In a second study, we examined the effects that mass pickup can have on the plasma 
temperature in an MHD approximation. Although previous studies have claimed that pickup 
can heat the torus plasma, we showed that such thermal pickup effects are controlled by 
the local sonic Mach number. The process is nonlinear because pickup can slow the flow 
appreciably. Thus, less heating than anticipated, or perhaps even cooling, may occur locally 
near Io as a result of ionization (#46). 
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APPENDIX II: UCI STT Numerical/Simulation Codes 
Developed Under NASA Solar Terrestrial Theory Grant NAGW-93 

A.II.A. Inverse-Power MHD Linear-Instability Code 

Computational studies of the behavior of ideal and resistive plasma instabilities in a 
sheared magnetic field often encounter a problem of conflicting length scales; in particular, 
the characteristic width of an unstable perturbation is often much smaller than the 
characteristic length (e.g., the shear scale a), of the equilibrium. One strategy is to retain 
second-order accuracy in the finite-difference expressions for first- and second-order 
derivatives by employing an appropriate coordinate transformation (e.g., Van Hoven and 
Mok 1984). We make the change of variable r = sinh -1 [6(j/ — y c )/a\, where y c marks the 
region of high resolution and b is an adjustable parameter which controls the magnitude of 
the resolution enhancement. Away from the region y « y c , r effectively represents a 
logarithmic scale. This transformation is more flexible than the one used by Van Hoven 
and Mok (1984) since it remains well behaved for y < y c . 

The full linearized system of MHD equations may be expressed as 

dX/dt + C-X = 0 (Al) 

where X is a vector whose elements are the perturbation variables and C is the 
corresponding matrix differential operator (in the transformed variable r). Working on an 
N-point spatial grid and using three-point spatial difference formulae, we write equation 
(Al) in finite-difference form at each grid point (i = 1,2 , . . . , A, ), implementing boundary 
conditions at the first and last grid points. To obtain exponentially growing solutions, we 
use an inverse-power method (Conte and de Boor 1980, p. 193); the growth rate of a 
desired solution is treated as an eigenvalue and estimated to be vq. Using a block 
tridiagonal algorithm (e.g., Roache 1976, pp. 345-349) to solve the set N x 7 simultaneous 
MHD equations, iteration is then performed according to the relation 

X N+1 = {£ + v 0 l)~'X N 
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where I is the identity matrix. The vector is found in general to converge rapidly to 
the eigenvector X, with corresponding eigenvalue u, if u Q is initially chosen to be 
sufficiently close to v. 

A.II.B. Nonlinear 1-D Thermal and Gravitational Hydrodynamics Code 

It is well known that the solar coronal plasma is locally unstable to rapid cooling and 
condensation due to energy losses from- optically thin radiation. Throughout most of the 
corona this effect is mitigated by thermal conduction; thus the corona is mostly hot and 
tenuous. However, cool, dense plasma is observed to form in highly localized regions 
associated with sheared magnetic field configurations. An understanding of the thermal 
stability of plasmas under such conditions is complicated by the complex, 
three-dimensional, magnetic geometry, and by the thermal connection to cooler, denser 
regions by means of magnetic field lines that thread both the corona and the 
chromosphere. Theoretical progress can be made by considering several related model 
problems with simpler geometry and boundary conditions. For several years we have 
studied radiative instabilities in two-dimensional sheared magnetic geometry. These studies 
have shed light on the plasma dynamics transverse to the magnetic field, and on the role 
played by the sheared field in thermally insulating the cool plasma from the surrounding 
hot corona. Such models say nothing about the thermal connection to the cool and dense 
chromosphere provided by dynamics parallel to the magnetic field. 

Recently we have developed a fully nonlinear, one-dimensional, gravitational 
hydrodynamics code to study these Thermal-connection effects. As this model considers 
only parallel flows, magnetic field evolution is ignored, as are curvature effects. The ends of 
the field line are considered anchored in the chromosphere, where a variety of boundary 
conditions (e.g., thermal insulation, inflow/outflow, constant temperature, exponential 
density, etc.) can be applied. In addition to the usual adiabatic terms, the energy equation 
includes thermal conduction (« ~ T 1 *), optically thin radiation, and variable coronal 
heating. The flows respond to pressure and (parallel) gravity forces. A simple artificial 
viscosity (proportional to the minimum mesh spacing) is used to prevent nonlinear 
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numerical instability at small wavelengths. All physical terms are advanced by using a 
leap-frog algorithm on a non-uniform Eulerian grid. The equation of motion is modified 
with semi-implicit terms (proportional to the timestep, Harned and Schnack 1986) which 
allow for numerically stable computations to be performed on time scales dictated by 
interesting physical processes, rather than by the usual CFL condition. 

A.II.C. Nonlinear 2.5-D ADI Incompressible MHD Magnetic- Reconnection 
Code 

A pervasive phenomenon in solar and space physics is magnetic reconnection: the 
spontaneous annihilation of magnetic field and its corresponding conversion into plasma 


kinetic and internal energy. These processes are of interest because they occur on 

timescales much shorter than that indicated by simple resistive dissipation, yet much i 

longer than the timescale dictated by Alfven or sound waves. Numerical studies of the -fl 

nonlinear evolution of reconnection events thus require that special computational methods 

be used. Over the past few years we have developed a 2.5-dimensional resistive MHD code 

suitable for these studies. This code employs a coordinate transformation (App. II. A.) in a 

the direction transverse to the primary magnetic field to allow for the enhanced spatial 5 

resolution required by astrophysical values of the Lundquist number. The assumption of 3 


incompressibility eliminates sonic and magnetosonic waves, while retaining the essential 
reconnection physics (Furth et al. 1963, Van Hoven and Cross 1973). Joule heating, 
thermal conduction, and optically thin radiation are also included in the model. Potentially 
unstable spatial oscillations at the mesh size, arising from quadratic nonlinearities in the 
MHD equations, are controlled by a simple spatial filter (Schnack, Baxter, and Caramana 
1984). The time advance uses the Alternating-Direction-Implicit (ADI) method (Schnack 
and Killeen 1980) to compute on the long time scales defined by the reconnection process. 
Time centering is accomplished by Newton iteration of the nonlinear terms (Finan and 
Killeen 1981). In this way the time step is determined by the requirements of accuracy and 
efficiency rather than by that of numerical stability of the fastest waves. 
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A.II.D. Nonlinear 2.5-D Compressible and Resistive MHD Thermal-Instability 
Code 

The structure and operation of this code are described at some length in the Appendix of 
Sparks, Van Hoven and Schnack (#47). 

A.II.E. Nonlinear 2.5-D Resistive MHD Spectral Code 

This code was used for the first nonlinear simulation of the resistive tearing mode (Van 
Hoven and Cross 1973) and is described in that publication. It has since been partially 
vectorized and a predictor-corrector time advance installed (#38). 

A.II.F. Modified 3-D MHD Code in Spherical Geometry 

Our time-dependent, three-dimensional MHD code solves the nonlinear, compressible MHD 
equations in spherical coordinates using a two-step Lax-Wendroff method. This algorithm 
is useful for problems where advection of flow is important, including the capturing of 
shocks and discontinuities (Richtmyer and Morton 1967). The code was originally 
developed to study three-dimensional plasma flow past Io (#37 gives a brief description), 
but has been adapted for solar problems. We have considered 2-D problems with azimuthal 
symmetry (#48) and full 3-D applications (#58). Our code can be used either for 
problems where boundaries are imposed at arbitrary values of r, 6, and <f>, or for global 
problems where the only boundaries occur at the solar surface and an outer radial location 
(solution at the poles of the spherical coordinate system is then accomplished by 
differencing a limit form of the MHD equations, similar to a form of the gas-dynamic 
equations derived by Bohachevsky and Mates 1966). For future problems, we plan to add 
semi-implicit terms to the code (Harned and Schnack 1986) to allow us to also follow the 
slow build-up of energy in the magnetic field while using large time steps. These terms are 
automatically turned off during the dynamic portions of the computation. 

A.II.G. Modified 3-D MHD Code in Cartesian Geometry 

The structure and operation of this code, which was used in #40, are described in the 
Appendix of Mikic et al. (1988). 
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A.II.H. 3-D MHD Code in Cylindrical Geometry 

We have recently developed a 3-D cylindrical MHD code which is similar to the Cartesian 
code described by Mikic, Barnes, and Schnack (1988). It has a Fourier representation in 
the azimuthal coordinate, with finite differences in the r and z coordinates. It has been 
applied to the dynamics and stability of coronal loops (#52). 


APPENDIX III: Abstracts of UCI STT Publications 
Supported by NAGW-93 

(see following pages) 
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Stability of a diffuse linear pinch with axial boundaries 

Giorgio Enaudi 

latmto di Flnca. Unimrtud di Pisa. Pita. Italy 
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A formulation of the stability behavior of a finite- length pinch is prnented. A general initial perturbation is 
es pre s sed as a uniformly convergen t sum over a complete discrete A set. A variational calculation is then 
p erfo r m ed, based on the ener gy principle, in which the end-boundary conditions appear as constraints. The 
requisite Lagrange multipliers mutually couple the elemental periodic excitations. The resulting extended 
form of 6W still admits a pr op er seoood-variadon treatment so that the mnunuxarioe and stability 
consi d erati ons of Newcomb remain applicable Comparison theorems are discussed as is the relevance of this 
end-effect model to the stability of solar coronal loops. 


I. INTRODUCTION 

Q is well known 1 that the study of the hydromagnetie 
stability of an equilibrium configuration can be reduced 
to the solution of the linearized equation of motion 



where i is the Lag rang tan displacement of a fluid ele- 
ment from its equilibrium position. The force operator 
F can be written in the form 

FU)»VF/> l (V*t)+T(fT)F, 

t<FKQ)xB, + (FJCB,)icQl, 

where p«, P%, and B,, respectively, are the equilibrium 
density, pressure, and magnetic field, 7 is the ratio 
of specific heats and Q= V* (in B») is the first-order 
perturbation in the magnetic field. By solving (1) with 
the proper boundary and initial conditions, it is possi- 
ble, in principle, to follow in time any small motion 
about an equilibrium state. However, one can deter- 
mine if a system is stable or not in an easier way by 
using an energy principle, 1 * 1 if F is a self-adjoint op- 
erator. To do so, one uses the variation in kinetic 
energy UK (a volume integral of j^Pe I € I a ) the 
change in the potential energy, which is a quadratic 
form 

aJF=-|/ d T «.F ( t) (2) 

not involving f because of the form of (1). From (1) and 
(2) it follows that 

0 K= l - f dr[|- F({)- {• F({)] 

=«»;({,«)- 6 tF.U.i), 

where 

f do A* {*• [(«' •?)/>„ + >#**{'-B- O'] 

+ B{{* O')} 

is the surface contribution to 6 W when { = {', S is the 


boundary of the plasma, and ii is the outward unit 
normal to the surface. It is clear that the requirement 
for the operator F to be self-adjoint depends upon the 
boundary values t has to satisfy. The perturbations we 
have chosen provide this condition in the form 0 w, = o. 

Newcomb 1 investigated the stability of a cylindrical^ 
symmetric and infinitely long configuration of plasma 
and magnetic field (a diffuse pinch) bounded on the out- 
side by a perfectly conducting wall. In this case 0W t is 
obviously zero. We went to investigate the stability of 
a similar structure, except that we introduce an in- 
homogeneity along the axis; that is, we consider the 
equilibrium density to have the form 

Jp t , |*| *L (I), 

P#( * )ss U»Pt, M> L (ID, (3) 

where 2L is the length of the region (I) of interest. 

Since the Alfvfa speed v 4 is inversely proportional to 
pj^ 1 , it is much larger in region I than in region II; it 
follows that, on a short time scale, we can consider the 
perturbation { to be zero in regions II, and thus require 
it to satisfy specific boundary conditions at * =* L. In 
order to solve the problem using an energy principle 
these conditions have to be such that 6W, is zero at 
the "ends”. A possible choice is 


{ i (-L)={ i (L)=0, 

(4a) 


(4b) 


(4c) 


ii and it being, respectively, the components of { per- 
pendicular and parallel to B». Notice that, since condi- 
tions (4) have to hold at every time, the boundary condi- 
tions on ( are the same. It ia easy to see the physical 
meaning of conditions (4) by considering the perturbation 
in the magnetic field Q. Because of (4a), Q t vanishes at' 
* * ±L and, therefore, the normal component ofthe mag- 
netic field across the boundary between fluids I and □ is 
continuous even if Q vanishes in fluid IL It follows that 
condition (4a) is the necessary and sufficient condition 
to model the inertial anchoring of the lines of force of 
the magnetic field at the surfaces z = ± L. This an- 
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A boundary-coupled generalization of the Newcomb stability criterion 

Alak Ray* 1 and Gerard Van Hoven 
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In magnetohydrodynamic problems in which there are boundaries in the axial or azimuthal 
symmetry directions, the energy principle involves multiple nonorthogonal trial perturbations. 

The subsequent Euler-Lagrange coupling of these allowed excitations necessitates a 
generalization of the Newcomb necessary-and-sufficient stability criterion. This extension, in 
which the absence of a conjugate point (instead of a simple zero) provides stability, is described in 
this paper. An important aspect of the proof involves the treatment of certain asymmetries which 
arise in the presence of the anisotropy caused by the magnetic field. The general method described 
here has applications to finite-length axial laboratory pinches and to astrophysical plasmas with 
rooted magnetic fields. 


I. INTRODUCTION 

The Newcomb criterion 1 provides a straightforward 
test, based on the Energy Principle, 2 for determining the sta- 
bility of a periodic, or infinitely long, cylindrical pinch. Sig- 
nificant simplification is obtained, in Newcomb’s theory, by 
expressing the variational trial functions in terms of the orth- 
ogonal complete set 

|*(r) = %\r\m,k )exp[i(mfl + kz)), 111 

which fits the pinch symmetry. One then determines ideal 
magnetohydrodynamic stability by the absence of zeros of 
the properly initialized Euler-Lagrange solution £ r (r,m,k ) 
for all m and k, a test which can be performed by a straight- 
forward numerical procedure. 3 

If, however, such a linear pinch has a limited extent in 
either symmetry direction, with specified boundary condi- 
tions, then at least two test functions of the form of 1 1 ) are 
required, and they become coupled in the minimization pro- 
cedure. 4 Situations of this kind, which provide improved sta- 
bility 5 by restricting the allowed perturbations, arise for fin- 
ite-length line-tied pinches 6 and for arcade-like half pinches 7 
each of which models structures in the solar corona. 

In this paper, we generalize the primary Newcomb cri- 
terion 1 to the case of a bounded, magnetized plasma. The 
boundary conditions both quantize the wavenumber k and 
destroy the energy-integral orthogonality of the members of 
the allowed trial-perturbation set, including that given in 1 1 ) 
and other alternative complete sets. 6 This situation then 
leads to coupled Euler-Lagrange equations, whose solutions 
provide a test for stability. Such problems are known in the 
Calculus of Variations le.g., Gelfand and Fomin 8 ), but the 
available criteria do not apply to the situation in which heli- 
cal magnetic fields thread the plasma. In this case, the resul- 
tant anisotropy destroys the ik.k 'I— i — k, — k') symmetry 
of certain coupling coefficients, which is taken as an essential 
assumption of previous proofs. 

Section II briefly describes the energy-principle formu- 
lation required for bounded plasmas. Section 111 provides 
the statement and a proof of the necessary generalization of 
the Newcomb criterion Finally, Sec. IV discusses applica- 
tions and conclusions. 

‘ Present address Tata Institute of Fundamental Research. Homi Bhabha 
Road. Bombay 400005. India. 


II. ENERGY PRINCIPLE FORMULATION 

In this section we will follow the treatments of Bern- 
stein et al. 2 and Newcomb, 1 substantially using the notation 
of the latter. The energy principle states that a plasma is 
stable if the (potential) energy integral 

<5 2 ^= -yJ<f 3 r|- F( |) (2) 

is positive, where F is a self-adjoint operator representing the 
ideal magnetohydrodynamic force 2 , and is a function of 
equilibrium quantities. Newcomb applied this method to the 
linear pinch and, after expressing the displacement £(r) in an 
orthogonal set as in ( 1 ), he minimized 5 2 If' with respect to 
and £, . The resulting form 

aiK'-y [ d r[f(i' k ) 2 +gii], (3) 

given in terms of \r) = $ r \r,m,k ) [where £, la) = 0 = f,lb)] 
and $ l = dg r /dr, is minimized by the solutions of the 
Euler-Lagrange equation 

(A;)'-#* =o (4) 

for each value of m (suppressed) and k. Newcomb then 
showed, 1 in a very careful way, that 

6\W>0 cb £ t lr)^0 (5) 

for 0 < a < r < b (or in certain independent subinterval.; ther- 
eof), with $ k (r) having proper initial values at r — a. This 
allows the determination of stability to be made from a sim- 
ple numerical integration 3 of (4). 

When a finite plasma is considered, however, with 
specified axial or azimuthal boundary conditions, the New- 
comb formulation can no longer be used. One way of seeing 
this is that the functions / and g are not even in k and m. 
Thus, a resonant “standing wave” solution [£ 4 = | - *. in the 
notation of ( 1 )], which will satisfy end conditions on a finite- 
length pinch, for example, cannot be constructed from (4). 
There are also phase relations' among the components of 
which prevent their vanishing simultaneously on specified 
surfaces. 

To circumvent this problem, one must either constrain 
the variation of (2) as described by Einaudi and Van Hoven 4 
or devise special forms of the minimizing test functions 
which have the boundary values built into them. 6 
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THE STRUCTURE, STABILITY AND FLARING OF SOLAR 
CORONAL LOOPS 


Gerard VAN HOVEN 
Department of Physics 
University of California 
Irvine, Calif 92717. USA 


Abstract 

A review is given of recent progress in the theory of the 
magnetohydrodynamic behavior of coronal loops, beginning with 
a brief characterization of the observations. The equilibrium 
magnetic field is described, along with the consequences of 
the empirical requirement for short-term, or infinite-conduc- 
tivity, stability which is shown to be dominated by the 
end-effect influence of the quasi-rigid photosphere. A new 
loop-flare model is then developed, which takes account of the 
finite loop length. The primary resistive-sausage-mode insta- 
bility exhibits the necessary threshold behavior, and produces 
a number of spatially and energetically distinct flare-release 
manifestations . 
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The growth of the tearing mode: Boundary and seating effects 
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The linear development of the resistive tearing instability in a sheet pinch is investigated 
numerically. Particular emphasis is placed on effects which differentiate magnetic tearing in 
astrophysical situations from that in laboratory devices. These include extreme values of the 
parameters determining the mode growth and a variety of boundary conditions. Eigenfunction 
profiles for long and short wavelengths are computed and the applicability of the “constant V" 
approximation is investigated. Nearby conducting walls tend to validate this condition and 
reduce the growth rate, especially for the long wavelength modes which, otherwise, disturb a 
larger region of the plasma than do short wavelength modes. Finally, the growth rate p is 
computed for values of the magnetic Reynolds number 5 up to 10 12 and of the dimensionless 
wavelength parameter a down to 10~ 3 . The results demonstrate, without approximation, theS 2/5 
scaling ofp at large a (constant IP ) and the S 2/3 scaling at small a (nonconstant* IP). The a and 5 
variation of the growth maximum, which would provide the dominant excitation in the absence of 
nearby boundaries, is specified for both single- and multiple-tearing layers. The growth maximum 
is shown to occur in a parametric regime where the constant-lP approximation is not valid. 


I. INTRODUCTION 

The study of magnetic reconnection has a long history, 
dating from Dungey’s proposal 1 of such a process to explain 
the energy release of solar flares and other astrophysical ex- 
plosions. Two primary reconnection mechanisms have been 
developed; one a true instability 2 and the second a static, 
externally driven, resistive diffusion. 3 

The former mechanism, resistive magnetic tearing, has 
been treated most often with the original analytic boundary- 
layer theory of Furth, Killeen, and Rosenbluth 2 (cited here- 
after as FKR) and has been used most extensively for labora- 
tory-fusion applications, often in cylindrical geometry. 4 
Many of these calculations depend upon certain ordering 
assumptions 3 or on a long-wavelength approximation. 4,6 

Previous investigators using numerical computa- 
tions 7-10 have not particularly investigated the adequacy of 
these analytic simplifications, nor have they covered a wide 
range of parameter space. Many recent workers 11,12 have 
concentrated on the study of nonlinear effects 13,14 over a 
practically limited span of physical conditions. 

We have several near- and far-term motivations in the 
present work. We wish to assemble a complete and self-con- 
sistent numerical model of resistive magnetic tearing in or- 
der to (a) verify and relate the results of the principal ap- 
proximations used in analytic analyses, and (b) to investigate 
the solutions and their growth-rate scalings 2,4,6 over a large 
range of the primary parameters which include parametric 
values applicable to the solar atmosphere. 13 The computa- 
tions described here cover the linear behavior for a variety of 
boundary conditions, as a first step to a series of numerical 
calculations probing the nonlinear evolution and eventual 


“Permanent address: Department of Physics. University of California, Ir- 
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saturation behavior of linearly unstable modes. 13,14 We are 
also, ultimately, interested in astrophysical and energy-re- 
lease applications of the tearing mode, 13 for which the mag- 
netic Reynolds numbers and wavelengths are larger than, 
and the boundary conditions and relevant diagnostics differ- 
ent from, those used for laboratory purposes. 

In the present paper we employ a model which uses 
several of the same assumptions as FKR. 2 Our approach to 
obtaining a solution to the dynamic equations differs from 
that used by FKR, however, in that we solve the time-depen- 
dent equations numerically throughout the entire region and 
do not follow a boundary-layer approach in which solutions 
applicable to two different regions must be matched at an 
appropriate interface. We, therefore, do not require either a 
"constant- IP” 2 or a long-wavelength approximation, 4 al- 
though we examine solutions covering the range from where 
constant V would be expected to be valid (relatively short 
wavelengths) to where it should not be valid (long wave- 
lengths). By computing simulations over a wide span in S 
(10 2 <5< 10 12 ) and a (10 -3 <ar < 1, where a is inversely pro- 
portional to the wavelength) we are able to verify various 
scaling laws for the growth rate which have appeared in the 
literature, 2-4,6 without benefit of a boundary-layer approach 
and any related subsequent approximations. 

II. THE RESISTIVE EQUATIONS AND THEIR 
NUMERICAL SOLUTION 

The basic equations 2 and the numerical procedure 16 
used to solve them are well-established, so they will just be 
briefly reviewed here. We assume that collisional magneto- 
hydrodynamic (MHD) theory 17 is applicable, that the plas- 
ma is incompressible and inviscid with constant isotropic 
resistivity t), and that gravitational effects are negligible. The 
necessary equations can then be written in cgs units as fol- 
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ABSTRACT 

We report on a numerical study of the interaction between a thermal wind and a global dipole 
field in the Sun and in a giant star. In order for closed field lines to persist near the equator (where a 
helmet-streamer-like configuration appears), the coronal temperature must be less than a critical 
value T e . which scales as M/R. This condition is found to be equivalent to the following: for a static 
helmet streamer to persist, the sonic point above the helmet must not approach closer to the star than 
2. 2-2. 6 stellar radii. Implications for rapid mass loss and X-ray emission from cool giants are pointed 
out. Our results strengthen the case for identifying empirical dividing lines in the H-R diagram with a 
magnetic topology transition locus (MTTL). Support for the MTTL concept is also provided by 
considerations of the breakdown of magnetostatic equilibrium. 

Subject headings: hydromagnetics — starsxoronae — starsrwinds 


I. INTRODUCTION 

The main question which we wish to examine here is 
he following: What condition must be satisfied in order 
hat dosed magnetic field lines can persist in steady 
’■ate in a stellar atmosphere? 

Recent empirical data provide motivation for this 
study. Thus, across a certain dividing line among giant 
stars in the H-R diagram, mass loss seems to become 
tforc rapid (Stencel and Mullan 1980), while X-ray 
'mission essentially disappears (Ayres etal. 1981). Since 
s°lar data indicate that strong X-ray emission is associ- 
J tcd with closed magnetic loops, while mass loss is 
Permitted only along open field lines, it is tempting to 
scribe the empirical dividing line in the H-R diagram 
10 a change from closed to open magnetic field topology 
'Mullan 1981. 1982). 

The fundamental process in the present discussion is 
’he mechanism bv means of which fields evolve from 
d°sed to open. Fields begin their evolution in a stellar 
atmosphere as closed lines, with footpoints rooted in the 
'tar. but if mass loss is to occur, the link with the star 
must eventually be severed. The severing process proba- 
bly involves more than a single stage. Thus, a loop may 
undergo MHD evolution at first, while finite resistivity 
J nd plasma effects may not play important roles until 
'atcr in the evolution. In the absence of detailed knowl- 
c d?e of the latter effects in stellar atmospheres, we here 
Cn nfine attention to the first stage of evolution, in which 
'vleu! MHD modeling turns out to be adequate to allow 


us to draw an important conclusion concerning the 
ability of closed field lines to persist. 

A preliminary answer to the main question posed 
above can be deduced on the basis of Pneuman's ( 1 968) 
analytic work on static helmet streamers. This work is 
summarized in § II. A more genera] MHD treatment of 
the problem is undertaken here. The method is de- 
scribed in § III. Results for stars of two different 
gravities are described in § IV. We find some numerical 
differences from the preliminary conclusions, but the 
preliminary conclusions are not radically altered. Impli- 
cations for stellar mass loss are briefly enumerated in 
§ V where we confirm the earlier suggestion that the 
empirical dividing line involves a transition in magnetic 
topology. In § VI we examine the dividing line from a 
different standpoint (namely, breakdown in magneto- 
static equilibrium), and we find corroboration for the 
magnetic topology transition interpretation. Conclusions 
of the work are presented in § VII. 

ii. pneuman's work on static helmet 

STREAMERS IN THE SUN 

The magnctohvdrostatic equilibrium of a helmet 
streamer (consisting of closed field lines below, sur- 
rounded by open field above) was studied analytically 
by Pneuman (1968). Various simplifying assumptions 
were made, including isoihermalitv and a total area 
function (equal to the sum of closed and open field 
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ENERGY DYNAMICS IN STRESSED MAGNETIC FIELDS. 

THE FILAMENTATION AND FLARE INSTABILITIES 
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ABSTRACT 

The thermal and tearing instab ilities are believed to be the two primary temperature modification 
mwhaniatm in sheared astrophysical magnetic fields. The former gives rise to the formation of cool 
filam ent* and the latter to the release of magnetic energy. It has long been known that these 
processes are interrelated, most conspicuously in the case of the solar corona where prominences 
often precede flares wi thin the same magnetic structure. It is also clear, from first principles, that the 
energy transport undertying the thermal instability should have a strong effect on the resistivity 
which facilitates magnetic taring and that the energy release of the latter should affect the 
temperature drop of the former. This paper describes some results of the first calculations which 
attempt to unify the dynamic treatment of these two coexisting instabilities. Growth rates as a 
function of resistivity, and examples of the primary mode structures are provided, along with a 
discussion of some critical aspects of the interaction of these two astrophysical energy flux 
mechanisms. 

Subject heading j: hydromagnetics — Sun: flares 

I. INTRODUCTION 

It is a familiar fact, as observed in the solar example, that increasing magnetic Geld shear or stress gives rise to the 
formation of filaments (Chiuderi and Van Hoven 1979) and to flares (Van Hoven 1979). Theories of the former 
mechanism (Field 1963), which is driven by a radiation output that decreases with temperature, have usually ignored 
the resistive magnetohydrodynamic effects of the resulting, very collisional. relatively low-temperature plasma. 
Treatments of the latter instability (Furth, Killeen, and Rosenblutb 1963), which is catalyzed by finite conductivity, 
have usually ignored the energy flux consequences of the resulting magnetic energy release. 

Since, in astrophysical situations, radiation is a strong effect and the relevant resistivity is the temperature-depen- 
dent Coulomb value, these two dynamic processes should be treated in a unified way. This paper represents the initial 
step in such an investigation, which will eventually describe the nonlinear development of the dynamics of a resistive, 
energy- transforming, plasma region subjected to magnetic stress. The present work concentrates on the numerical 
solution of the linear instability problem, which has also been treated by an analytical method (Steinolfson 1983) that 
is able to predict approximate growth rates. 


It. FORMULATION 

We consider a current-carrying, stationary volume of plasma whose magnetic Geld is described by the model form 

flo[e T tanh(y/a)+e.sech(>-/fl)] - So**. (I) 

which is force-free and thus consistent with an initially uniform temperature and density. We describe the development 
of this system by using the resistive MHD equations, relating flow velocity to magnetic induction, along with the 
incompressible energy transport equation. 

Kp%-(y-')[H(p)+l(T)J 2 -f> 2 *(T) + VK n (T)e h erVT}, ( 2 ) 

specifying the variation of temperature. Here the notation is the same as that used in Chiuderi and Van Hoven (1979). 
except that K is Boltzmann's constant per unit (average) mass, an Ohmic heating term is added, and the temperature- 
dependent Coulomb values (Spitzer 1962) of resistivity tj and (dominant) parallel thermal conductivity k„ are used. 
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TURBULENT RESISTIVE HEATING OF SOLAR CORONAL ARCHES 
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ABSTRACT 

Field-aligned currents in coronal arches can heat the plasma through instabilities if the electron 
drift velocity is large enough. We consider the heating as a time-dependent problem in a plasma 
cylinder, including heating rates for ions and electrons, longitudinal conduction and convection, 
radiation, and cross-field transport, all using both Coulomb and turbulent effects. Electron drift 
velocity is held constant, resembling an arch with high inductance. The electrostatic ion cyclotron 
instability proves the dominant effect at all times, starting from an initially isothermal plasma at 
T - 10 eV. Ion acoustic modes never become unstable, since T e /T j never exceeds 3. Rapid heating to 
temperatures of 100-1000 eV occurs, after which a slow “percolation” sets in, with turbulence 
switching on and off rapidly, keeping the plasma temperature in the - 100 eV range. This steady 
tradeoff between weak microturbulence and macroscopic loss processes determines the long time 
heating level. Turbulence is described with a resonantly broadened saturation model which enjoys 
some experimental support. Radiation losses are never important. High density arches heat more 
slowly (-10 s). Finite plasma pressure effects and turbulent cross-field transport have little influence 
on the long term temperature. Cases with a turbulent Ohm’s law (low inductance case) also yield 
strong heating rates. Despite its intrinsic unsteady behavior, resistive heating seems capable of 
providing efficient arch heating. 

Subject headings: hydromagnetics — plasmas — Sun: corona — turbulence 


I. INTRODUCTION 

Coronal heating models depend either upon energy 
deposition by waves traveling up from the base of a 
coronal loop, or upon currents running from one foot of 
a loop to another. This paper studies the latter possibil- 
ity, using anomalous Joule heating by electrostatic ion 
cyclotron waves. 

In the wave model it is difficult to allow for boundary 
losses and convective effects. In the current-generated 
turbulence picture, however, careful balancing of losses 
with turbulent heating rates can be done. (Tucker 1973; 
Rosner. Tucker, and Vaiana 1978; Nolte, Petresso. and 
Gerassimenko 1979; Vlahos 1979). Earlier work often 
assumed that Buneman or ion-acoustic modes were ex- 
cited in the loop (Hinata 1979). Typically these waves 
heat the coronal plasma very quickly, produce high 
resistivity, and lead to large changes in transport phe- 
nomena in a very short time ( < sec). This may be 
relevant to strongly driven currents and flares, but it 
seems doubtful that quiescent heating follows this pic- 
ture. 

Here we study current-driven heating without making 
any assumptions about which mode is dominant. In- 
stead, we begin with a plasma in a cylinder with T e = Tj, 
and turn on an electron drift at time t = 0. The ends and 
sides of the cylinder are held at a constant temperature 


(usually 10 5 K). We allow currents along z. the direction 
of the magnetic field, and further allow the plasma 
cylinder to expand radially (due to finite plasma pres- 
sure). We assume as a basic condition that the self-fields 
from these heating currents are negligible. A specific 
picture of how this could be so, Figure 1, portrays 
currents which reverse direction on scales of about 100 
m. This minimizes the energy lodged in self-fields. Our 
calculation includes turbulent transverse transport, 
longitudinal conduction, convection, and radiative losses, 
based on existing laboratory and theoretical work (see 
§ III). Hinata (1980) wrote down similar equations and 
speculated on possible solutions. 

In the absence of inhomogeneities in the equilibrium, 
there are no unstable drift modes. Because the plasma in 
the loop starts out isothermal, the first instability excited 
by relative electron-ion drift is the electrostatic ion 
cyclotron mode, denoted IC. The critical electron drift 
speed for exciting this instability is u c = 0.3 (T t /T e ) 3/2 v e , 
which is always less than the electron thermal velocity. 
In all our solutions, the equivalent electric field driving 
the currents never exceeds the runaway Dreicer field, so 
we can use drifting Maxwellian distributions throughout 
the temporal evolution. The ion cyclotron instability 
differs from the more virulent ion acoustic and Bune- 
man modes in that it is in principle self-regulating, since 
u c is proportional to T} n /T t . Typically, in times less 
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MICROWAVE SIGNATURES FROM A RECONNECTING PLASMA 
PINCH. WITH APPLICATION TO LOOP FLARES 1 
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ABSTRACT 

We calculate the microwave signature of a cylindrical plasma pinch undergoing magnetic recon- 
nection. a process which occurs in many astrophysical situations, such as solar flares. Depending on 
the viewing angle and the average energy of the accelerated electrons, the microwaves from this 
betatron-like source show various amounts of circular polarization. The degree of polarization is 
shown to be frequency dependent, and the sense of polarization is sometimes reversed. The power 
spectrum is predicted to have several interesting properties, which can be compared with high-resolu- 
tion measurements. 

Subject headings: polarization — Sun: flares — Sun: radio radiation 

I. INTRODUCTION 

Microwave radiation is observed in a variety of astrophysical situations. Among many possibilities, such emission 
can come from energetic electrons traveling in a magnetic field, producing gyrosynchrotron radiation in the microwave 
range. The origin of these electrons, however, is less well understood. Depending on the physical conditions, there can 
be a number of acceleration mechanisms, due to turbulence. DC electric fields, etc. Various types of these energization 
processes have been used to provide microwave emission models. 

Generally speaking, for an observed gyrosynchrotron source, a particular energetic-electron model can explain 
certain characteristics of the microwave spectrum. A more complete theory, however, must also consider the origin of 
these electrons and the mechanism that accelerates them in the plasma environment and magnetic structure derived 
from measurements. Since field reconnection instabilities occur in numerous astrophysical conditions, it seems natural 
to investigate how they give rise to energetic electrons and to calculate their microwave signatures, so that they may be 
compared with observations. 

In the present work, we examine a specific example based on magnetic field reconnection in a magnetized plasma 
pinch (Mok and Van Hoven 1982). In this configuration, the instability produces a quasi-DC electric field (Van Hoven 
1979) in the azimuthal direction, which happens to be parallel to the (azimuthallv) closed magnetic field lines. The 
acceleration is, therefore, extremely efficient so that the effective radial size of the plasma column does not have to be 
large. In fact, this can be one of the most compact microwave sources in the domain of astrophysics. 

To cite a specific application, we choose a simple bipolar magnetic loop, with physical parameters corresponding to 
those observed in solar active regions. The triggering of a solar flare is considered as the onset of the magnetic 
reconnection instability. The instability criterion and a detailed analysis of such a physical process have been published 
(Mok and Van Hoven 1982). In this paper, we will confine ourselves to the emission of microwaves as a consequence of 
the reconnection process. Although we use a solar flare as an example in the present computation, the principles and 
methods can be applied to various astrophysical objects under appropriate conditions. 

Microwave emission in solar flares has been widely observed and theoretically interpreted (Marsh. Zirin. and 
Hurford 1979: Marsh etal. 1980, 1981; Marsh and Hurford 1980: Alissandrakis and Kundu 1978: Svestka 1976; 
Scalise, Basu. and Marques Dos Santos 1971; Tanaka and Enome 1970; Sturrock 1980; Emslie and Vlahos 1980: Dulk. 
Melrose, and White 1979; Matzler 1978; Bohme etal. 1977; Kundu, Velasamy. and Becker 1974; Takakura and Scalise 
1970: Holt and Ramaty 1969: Ramaty 1969). In this work, we calculate the microwave signatures of a reconnecting 
magnetized plasma pinch and suggest a plausible alternative to the existing loop flare models based on physical 
conditions in the solar corona. When compared with the location of the microwave source, the polarizations, and the 
radiated power, our results seem to be in agreement with observations. Since high (frequency) resolution measurements 
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Energetics and the resistive tearing mode: Effects of Jouie heating and 
radiation 
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Hydromagnetic theory is used to investigate the influence of Joule heating and temperature- 
dependent radiation on the linear development of the incompressible, resistive tearing instability 
in both force-free and neutral-sheet magnetic fields in a slab geometry. A temperature-dependent 
resistivity couples the energetics and the dynamics of the instability. Analytical expressions are 
derived for the growth rates utilizing, in separate analyses, the constant - V and the long- 
wavelength (nonconstant- V) approximations. The solutions indicate the occurrence of several 
modes in addition to the usual tearing mode. A number of the modes have relatively slow, 
complex growth rates. However, at large values of the magnetic Reynolds number 5 (S £ 10 7 , the 
actual value depending on the relevant parametric values), there are at least two modes with 
purely exponential growth (real, positive growth rates) when the radiative loss decreases with 
increasing temperature. One is the resistive tearing mode, and the other is a more rapidly growing 
mode with a growth rate resembling that of a thermal instability. If the radiation is neglected, the 
Joule heating alone also results in two modes with real, positive growth at large 5. One mode is 
again the tearing mode, but, in this case, the second mode is a more slowly growing Joule-heating 
mode. Another point of departure from previous theoretical results is that, below a particular 
value of S’ (SS 10 6 , the value once again being determined by the parameters), all the modes are 
generally stabilized. 


I. INTRODUCTION 

Resistive tearing instabilities have been invoked as a key 
constituent in several models of various naturally occurring 
phenomena in solar-flare physics' -3 and magnetospheric 
physics 4 (e.g., substorms, geomagnetic tail) and have been 
observed in numerous laboratory devices, particularly in to- 
kamak discharges. 3 The common occurrence of this instabil- 
ity and its potential importance in understanding phenom- 
ena in which it may occur have prompted a large theoretical 
effort, initiated by the work of Furth, Killeen, and Rosen- 
bluth 6 (referred to hereafter as FKR). The tearing instability 
has the potential of occurring at locations where a compo- 
nent of the magnetic field reverses direction. It results in a 
reconnection of the field lines about the initial location of 
field reversal with a subsequent conversion of magnetic ener- 
gy in the initial field to heating and particle acceleration. 

Despite the aforementioned and well-recognized 2 ’ 3 
close connection between the dynamics and energetics of re- 
sistive tearing, the linear dynamic evolution of the mode, 
which is of interest here, has previously been uncoupled 
from the thermodynamics. This is accomplished by assum- 
ing either that the resistivity is constant'- 7-10 or that only 
convective effects in the energy equation are important. 6 Di- 
biase and Killeen" do include a variable resistivity and a 
complete energy equation in their numerical study, but they 
do not include radiation and do not attempt to isolate the 
effect of the other energy term in which we are interested: 
Joule heating. Studies of the thermal instability 12 ' 3 do in- 
clude the energy equation with the heating and radiation, 
but they assume zero resistivity, and hence, do not allow for 
the competing development of the tearing instability. 


We have undertaken a study to investigate the contribu- 
tion of energy flux to the dynamics of magnetic field recon- 
nection, with the aim of unifying the resistive tearing' and 
thermal instability 13 mechanisms in a sheared magnetic 
field. We are pursuing this study from two separate, but 
complementary, approaches — the primary one being nu- 
merical' 4- ' 3 and the other, as in the present paper, analytical. 
Our model includes the effects of resistivity in Ohm’s law, 
and Joule heating, radiation, and thermal conduction in the 
energy equation. The numerical approach solves the model 
equations without any further approximations beyond lin- 
earization. The analytical approach, on the other hand, re- 
quires additional approximations and, as presented herein, is 
only valid for negligible thermal conductance. 

A temperature-dependent, Coulomb-like resistivity is 
considered in the present paper, thereby providing the cou- 
pling between the dynamics and the energy equation. In 
thermal instability studies the coupling is accomplished by 
compressibility. However, as is commonly done in tearing- 
mode analyses, the plasma is here assumed to be incompres- 
sible. We use magnetohydrodynamic (MHD) theory with 
Joule heating and temperature-dependent radiative loss in- 
cluded in the energy equation. 

Our object is to determine the influence of Joule heating 
and radiation on the linear development of the tearing insta- 
bility in slab geometry. The analysis closely parallels that of 
FKR for modes which satisfy the constant- approxima- 
tion. For long-wavelength modes which do not satisfy con- 
stant V, the approach of Coppi et al. 9 (hereafter cited as 
CGPRR) and Pritchett et al . 10 is followed. The differenti- 
ation between constant- and nonconstant - 'P solutions 
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A study is made of the effect of a temperature-dependent Coulomb-like resistivity on the planar 
tearing mode. The local evolution of the temperature is described by an energy equation which 
includes Joule heating and optically thin radiation. The resulting system of coupled linear 
magnetohydrodynamic equations is solved numerically, and eigenfunctions and growth rates are 
obtained. In the absence of radiation, there are two distinct solutions above a critical value of the 
magnetic Reynolds number S, a tearing-like mode and a Joule-heating mode. Below this point, the 
growth rates coalesce into a conjugate-complex pair. When stable radiation {dR /dT > 0) is added, 
the heating mode disappears and a modified tearing excitation exists to much lower values of S 
before its growth is cut off by Ohmic heating. Examples are given for solar coronal parameters, 
and for those characteristic of fusion-research devices. The introduction of an effective value for 
the resistivity, in the presence of energy transport, allows a simple qualitative discussion of the 
different modes. 


I. INTRODUCTION 

The resistive tearing mode has been a topic of consider- 
able interest since it was first proposed by Furth, Killeen, 
and Rosenbluth 1 (FKR). It has been the subject of linear 
analysis, showing two principal regimes of interest, the 
“constant- P ” case 1 and the small-wavenumber limit , 1,2 and 
an intermediate range where the growth rate is maximum . 3 
Connections between these regimes have been made 4-5 which 
show the relationship between the planar and cylindrical 
forms, and their nonlinear behavior has been examined . 6-9 

The essence of the physical cause underlying magnetic 
tearing is the spontaneous generation, in a sheared field, of a 
very narrow layer where accelerated resistive diffusion, or 
field reconnection, can occur. The rate of this process de- 
pends directly on the value of the resistivity rj in the tearing 
layer, scaling as 17 v where 0.4 < v < 0.67. It is a fact that the 
relevant resistivity, in most astrophysical and laboratory si- 
tuations, is the Coulomb value 10 i j = t/ c T ~ 312 , and that Oh- 
mic diffusion heats the ambient plasma. Thus, the question 
arises whether the localized Joule heating, caused by the re- 
connection, can quench the instability. When one examines 
previous attacks on this problem, one discovers instances 
where a variable resistivity has been included in the formula- 
tion , 18 - 9 but no cases where its dynamic effects have been 
carefully studied. Almost all previous calculations assume 
constant resistivity. 

In the present paper we examine the consequences of 
the use of a space-time-variable Coulomb resistivity to cata- 
lyze magnetic tearing. In order to describe the necessary 
temperature evolution, we use an incompressible (slow) form 
of the energy-transport equation, which includes the effects 
of Ohmic heating and its principal energy-loss competitor, 
optically thin radiation. The emission laws used have posi- 
tive temperature slopes (e.g., bremsstrahlung), so that we 
need not consider radiatively unstable modes . 1112 Finally, 
we do not include thermal-conduction losses , 1 1 since they 
are negligible at the relatively small wavenumbers which are 
relevant for magnetic tearing. A recent publication on the 


normal tearing mode , 5 a brief report on the full energy-trans- 
port formulation including unstable radiative modes , 12 and 
an analytical treatment , 13 which predicts the growth rates 
found here, serve as an immediate background to the present 
work. The latter paper 13 uses the constant- IP and small- 
wavenumber approximations to show the existence of the 
tearing, radiative, and Joule-heating modes. 

In the following sections we describe the formulation 
and the numerical methods (II), and the results of the com- 
putations (III), including both solar-coronal and fusion-re- 
search examples. Section IV contains a discussion of the en- 
ergy dynamics in terms of an “effective” value for ij and the 
conclusions derived from this study. 

II. EQUATIONS 

We consider a stationary distribution of plasma with 
uniform density and temperature throughout, the latter 
maintained against radiation by a constant heat input. The 
equilibrium plasma is thus force-free (has no pressure gradi- 
ents) and its magnetic field is described by 

B 0 (.v) = 5 0 [e, tanh {y/a) + e, sech(y/a)], ( 1 ) 

where e x (e, ) is a unit vector in the x(z) direction and a speci- 
fies the magnetic shear scale. 

It is then assumed that the magnetohydrodynamic 
equations properly describe the evolution of the plasma in 
terms of its magnetic field B, its flow velocity v, and its tem- 
perature T, and that the plasma is incompressible, inviscid, 
and free from gravitational effects. 

The equations in cgs units are given by 

\. y = V-B = 0, (2) 

iM. = VX(VXB) - 4 - VX MnVxB)]. (3) 

dt 4 tt 

v(VXv)= -Vx(v-V)v+ -L-VX[(VXB)XB], (4) 

at 4irp 0 

^f = ^[nJ 2 -nl0[T) + H{n o )], (5) 

dt 2k„n 0 
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ABSTRACT 

Two energy modification mechanisms which are known to occur in sheared magnetic fields are the tearing 
and thermal instabilities. These processes can be studied separately with formalisms incorporating just the 
effective driving mechanism of interest (finite resistivtty for the tearing mode and unstable radiation for the 
thermal mode). A model which includes both effects, and a temperature-dependent resistivity, indicates that 
modified forms of these two instabilities may coexist for identical physical conditions. When they are isolated 
computationally, one can show that their limiting growth rates are approximately those of the uncoupled 
instabilities. The spatial structure and energy content of these two new hybrid processes are then individually 
examined and are found to differ considerably from those obtained from separate treatments of the driving 
mechanisms. The faster radiative instability, which has a hydromagnetically scaled growth rate like the 
condensation mode of the thermal instability, is shown to involve a substantial amount of magnetic field 
reconnection. This can be partially explained by a large temperature drop (or resistivity rise) at the X-point. 
The island width of the Coulomb-coupled radiative mode is 30 u „ of that produced by a comparable level of 
the slower tearing instability. In additioa the perturbed magnetic energy in the radiative instability is S times 
that of the perturbed thermal energy, indicating an appreciable modification of the initial magnetic structure. 
Subject headings: hydromagnetics — instabilities — Sun: flares 


I. INTRODUCTION 

The resistive tearing instability (Furth. Killeen, and Rosen- 
bluth 1963) has long been recognized as a potential 
candidate for initiating the release and conversion of magnetic 
energy in sheared astrophysical fields. (See. for instance. 
Sturrock 1968; Coppi and Friedland 1971; Van Hoven and 
Cross 1973). For the solar flare application of interest here, 
one shortcoming of the original tearing instability is that the 
rate of energy conversion, based on observed large-scale 
quantities in the solar atmosphere, is somewhat slow. The 
predicted time scale can be shortened by using small shear 
scales (below the observable spatial resolution) or unexpected 
combinations and values of the thermodynamic quantities 
and magnetic field strength. Additionally, shorter time scales 
can be obtained from long-wavelength disturbances (Van 
Hoven 1976; Steinolfson and Van Hoven 1983), or by 
hypothesizing the existence of either multiple, interacting, 
tearing layers (Pritchett. Lee. and Drake 1980). or the onset 
of micro-instabilities near the tearing layer (Coppi and 
Friedland 1971). 

We suggest that the time scale of the instability and. 
therefore, of the energy release may be reduced by approxi- 
mately two orders of magnitude in solar atmospheric con- 
ditions by recognizing that the conversion may occur in a 
radiatively modified variation of the (dynamic) tearing 
instability introduced by Furth, Killeen, and Rosenbluth 
(1963). We propose that the energy conversion takes place in 
an instability which involves tearing but which has been 
substantially altered by energy transport, in particular, by 
Joule heating and unstable (emission varies inversely with 
temperature) optically thin radiation. This instability, called 
the radiative instability, is one of two purely growing modes 


that coexist when the energetics are coupled to the dynamics 
of the tearing mode. The existence and structure of the two 
modes have been described by Van Hoven. Steinolfson. and 
Tachi (1983) and their growth rates given in analytic form by 
Steinolfson (1983). As demonstrated in these studies, the 
growth of the faster radiative mode is similar to that of the 
condensation mode of the thermal instability (Field 1965) in 
the ideal magnetohydrodynamic limit, in a sheared field 
(Chiuderi and Van Hoven 1979). while the growth rate of the 
slower mode is similar to that of the usual tearing mode 
(Furth. Killeen, and Rosenbluth 1963: Coppi et al. 1976; 
Steinolfson and Van Hoven 1983) as obtained without 
inclusion of the energetics. Although the growth rate of the 
radiative mode is determined by the unstable radiation, it has 
not previously been established that it also involves appre- 
ciable reconnection at the same rate. It could have been that 
the radiative instability just produced alternating hot and cold 
spots (relative to the background), without any substantial 
reduction in the magnetic energy, as in the conventional 
thermal instability. 

Previous studies of the radiative instability have concen- 
trated on determining when purely growing modes exist, 
their growth rates, and their eigenfunctions. We now extend 
this study and make a detailed comparison of the spatial 
structure and the energy content of the two instabilities 
(tearing and radiative). In order to accomplish this in a 
meaningful manner, the disturbance amplitude of each mode 
must be selected such that the two modes are effectively equal 
in some sense. The relative amount of tearing and recon- 
nection in the two modes depends, of course, on the method 
used to select the amplitudes of the disturbances. We compare 
the modes when the total contribution from the nonlinear 
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RADIATIVE AND RECONNECTION INSTABILITIES: FILAMENTS AND FLARES 
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ABSTRACT 

An investigation of the energetically coupled tearing and radiative instabilities in sheared magnetic fields has 
been made. It has been reported elsewhere that, if a temperature-dependent Coulomb-like resistivity is used in the 
normal tearing-mode analysis, instead of a constant resistivity, it is possible to study the local evolution of the 
temperature in terms of an energy equation. When the energetics further include Joule heating and optically thin 
unstable radiation, their effect on the dynamics of the instability produces a pair of purely growing unstable modes. 
One is the usual tearing mode, and the other is the radiative mode, which exhibits growth similar to the well- 
known condensation mode of the thermal instability in ideal compressible MHD. The modes are distinctly 
separated at high values of the magnetic Reynolds number S. but, as this parameter is lowered, a mutual 
coupling between the modes appears, and the two growth rates finally merge at a critical value of S. below which 
no purely growing mode exists. This quenching of the unstable modes is shown to be strongly wavenumber 
dependent. A study of the eigenfunctions and their derived characteristics shows that the radiative mode exhibits 
a large temperature drop in the reconnection region, due to the overwhelming dominance of radiative cooling 
over Joule heating, whereas the tearing mode exhibits a much smaller drop in the central layer. The strong 
X-point resistivity increase, connected with the radiative-mode temperature perturbation, helps to explain the 
significant magnetic reconnection and field perturbation provided by this mode. For typical solar coronal 
conditions, the growth rate of the radiative mode is approximately two orders of magnitude larger than the 
growth rate of the tearing mode. Thus, the time scale of the energy conversion due to magnetic reconnection can 
be greatly reduced by the radiative mode. 

Subject headings: hydromagnetics — instabilities — Sun: corona — Sun: flares 


i. introduction 

The energy source of a solar flare is generally believed to 
reside in the strong, stressed, coronal magnetic fields existing 
above photospheric active regions. The magnetic tearing 
instability (Furth. Killeen, and Rosenbluth 1963) has been 
proposed as the mechanism for the release of this stored energy 
(Sturrock 1968; Coppi and Friedland 1971: Van Hoven 1976, 
1979). This resistive process allows the magnetic energy to be 
dissipated directly into particle acceleration and heating by 
collision-dominated magnetic field reconnection. Despite the 
fact that the dynamics and energetics of the instability are thus 
coupled together, the assumption of a constant resistivity in an 
incompressible plasma, which has been used in most tearing- 
mode analyses (Furth. Killeen, and Rosenbluth 1963; Cross 
and Van Hoven 1971; Coppi et al. 1976; Pritchett. Lee, and 
Drake 1980). does not allow a simultaneous study of the two 
important aspects of the instability, namely, the dynamics of 

reconnection and the energetics of Joule heating and radiation. 
Whereas the use of a Couiomb-like collisional resistivity 
(Spitzer 1 962) allows us to study the dynamics and energetics of 
the tearing mode, the particular temperature dependence of 
the resistivity suggests the possibility that local heating may 
introduce a self-quenching effect on the instability, thus raising 
a question as to its adequacy in modeling a solar flare. 

The first consistent study of the effect of the collisional 
temperature-dependent resistivity has been approached 
numerically by Van Hoven. Steinolfson. and Tachi (1983) and 

analytically by Steinolfson (1983). These reports briefly 

described the existence of a number of unstable modes, along 


with a comparison of their growth rates and eigenfunction 
profiles, which appear in the different cases where the energy 
equation contains (1) Joule heating only. (2) Joule heating and 
stable (d/dT > 0) optically thin radiation. (3 ) Joule heating and 
unstable ( d/dT < 0) radiation, and (4) the additional influence 
of the parallel (to B ) component of thermal conduction. Cases 
1-3 were treated with computational and analytical methods, 
while case 4 was studied numerically. The analytic approach 
considered both force-free and neutral-sheet magnetic con- 
figurations, while the computational approach concentrated 
on force-free equilibria. 

An extensive numerical study of cases 1 and 2 has been 
carried out by Tachi. Steinolfson. and Van Hoven (1983). who 
showed that, in the absence of temperature-dependent radiative 
power loss, there exists a pair of unstable modes (tearing-like 
and Joule-heating) which, below a certain low value ofS ^ 10 i. * * * * 6 . 
coalesce and become (slowly) overstable as a result of a mode 
coupling, so that no purely growing unstable mode remains. 
When stable (d/dT > 0) radiation is added, for the particular 
set of initial parameters chosen for the study, the effect is to 
remove the Joule-heating mode and to allow the tearing-like 
mode to persist to a much lower value of S % 10 : . This change 
results from the suppression, by the very effective radiative 
energy loss, of the temperature rise due to Joule heating. The 
mutual coupling and quenching of the mode, at low 5. appears 
when Joule heating overwhelms radiation. 

The results from the study of cases 3-4 with unstable 
(d/dT < 0) radiation were nearly identical since, in the latter 
case, only a long-wavelength limit was considered, and the 


Nonlinear evolution of the resistive tearing mode 
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The nonlinear behavior of the tearing instability is investigated with numerical solutions of the 
resistive, incompressible, magnetohydrodynamic equations. Simulations have been completed for 
values of the Lundquist number S from 10 2 to 10 6 and wavelength parameter a = 2-rra/k from 
0.042 to 0.5. The initial state for the nonlinear computations is provided by the linear instability. 
In all cases, the nonlinear mode initially evolves at the linear rate, followed by a period of 
considerably reduced growth. At high Sand low a, secondary-flow vortices, opposite in direction 
to the linear vortices, generate a new magnetic island centered at the initial X point. The nonlinear 
growth rate of the one constant-0 (in the linear regime} solution considered is approximately an 
order of magnitude less than that of a comparable nonconstant-0 solution over the same time 
period. The nonconstant-0 computations indicate a reduction of from 8% to 27% of the initial 
shear-layer magnetic energy, with the larger reductions occurring for the longer wavelength 
disturbances. The constant-0 simulation shows a reduction two orders of magnitude smaller. The 
island width of the nonconstant-0 solution becomes larger than twice the scale of the initial shear 
layer. For all cases, the electric field parallel to the magnetic field grows, at the end of the run, to 
about 10“ 3 times the Driecer field. 


I. INTRODUCTION 

Ambient plasma and magnetic field configurations, un- 
dergoing relatively slow temporal changes, are often sus- 
ceptible to purely growing instabilities, on a much shorter 
time scale, in response to small-amplitude disturbances. Lin- 
ear theory can frequently be used to determine the parameter 
space in which a particular instability exists as well as the 
initial growth rate. As useful as the linear results are, they 
are seldom capable of establishing whether the instability 
may ultimately be responsible for the large-amplitude, ob- 
servable consequences. Without relying on questionable ex- 
trapolations of the linear computations, the only way to ob- 
tain this information is through an examination of the 
subsequent nonlinear evolution. 

We consider the nonlinear phase of the resistive tearing 
mode. 1-3 This instability may be important in astrophysical 4 
and magnetospheric 3 phenomena, as well as in laboratory 
plasmas. 6 Our interest is primarily in solar applications, 
such as the flare. Previous nonlinear treatments have con- 
centrated on the energy output, 7 special analytic limits, 8 or 
have been limited either to restricted ranges of parameter 
space (which are not relevant to the flare) or to multiple lay- 
ers. 7-9 

Possible changes in the nonlinear structure and evolu- 
tion of the tearing mode over a range of parametric values 
are of particular interest. The influence of the linear condi- 
tions 3 (e.g., the constancy of 0) on the nonlinear growth char- 
acteristics is investigated. We also compute the astrophysi- 
cally relevant quantities of global energy release and local 
electric fields, which can provide observable diagnoses of 
such magnetic reconnection processes. 

We study the nonlinear evolution in slab geometry us- 
ing resistive magnetohydrodynamic (MHD) theory, with the 
main assumptions being those of incompressibility and neg- 
ligible viscous dissipation. The initially stationary plasma, 
with constant thermodynamic properties, is embedded in a 


force-free magnetic field. Our approach is to use a linearly 
evolving mode that has attained its maximum linear growth 
as the initial state for the nonlinear computation. This mini- 
mizes the nonlinear calculations, ensures that the nonlinear 
phase evolves from a truly linear disturbance, and sets the 
spatial scale for the simulation. Computations have been 
carried out for values of the Lundquist number 5 from 10 2 to 
10 6 and of the wavelength parameter a from 0.042 to 0.5. 

A principal result of our computations is that the non- 
linear evolution generally differs from one region of param- 
eter space to another, and, hence, a prototypal characteriza- 
tion of the evolution in the nonlinear regime is not possible. 
As an illustration of this, constant-0 1-3 linear modes evolve 
much more slowly and convert considerably less of the 
stored magnetic energy in their nonlinear phase than do 
those with nonconstant 0 2-3 (long wavelengths). Similarly, 
high-S’ ( 2 10 3 ), low-a ( S 0. 1 ) modes undergo a change in con- 
figuration space with the formation of a new magnetic island 
centered at the original X point, while the remaining solu- 
tions do not. 

Some general characteristics that apply to all solutions 
are: (1) The nonlinear spatial distributions of the physical 
variables differ substantially from the linear behavior; (2) 
once nonlinear effects become important, the growth slows 
considerably from the linear rate; (3) the growth rate is a 
poor predictor of the nonlinear performance of a particular 
mode in terms of magnetic energy conversion; (4) the major- 
ity of the stored magnetic energy is resistively converted to 
thermal energy ( > 80%) for the long wavelength modes with 
a < 0.5; and (5) maximum electric fields parallel to the mag- 
netic field are about three orders of magnitude less than the 
Driecer field. 10 

The computational procedure is outlined in the follow- 
ing section, and the numerical results are presented in Secs. 
Ill and IV. The final section includes some elaboration on 
our results, particularly with respect to solar flares, as well as 


1207 


Phys. Fluids 27 (5), May 1984 


0031 -91 71 /84/051 207-08S01 .90 


© 1 984 American Institute of Physics 


1207 


The Astrophysical Journal. 281:854-861. 1984 June 15 

* 1984. The American Astronomical Society. All rights reserved. Primed in U.S.A. 
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ABSTRACT 

The inclusion of parallel and perpendicular thermal conduction in the analysis of the tearing and radiative 
instabilities produces modified versions of these two modes. This improvement of the energy-transport descrip- 
tion also introduces a large number of new instabilities whose growth rates lie between those of the thermally 
conductive tearing and radiative modes. With the exception of the modified radiative mode, all of the modes with 
conductive heat transport contain spatial temperature oscillations, perpendicular to the tearing surface, on a scale 
comparable to the width of the resistive tearing layer. These thermal ripples arise in a region where perpendicular 
thermal conduction and unstable, optically thin radiation are the dominant energy-transport mechanisms. This 
region extends from the edge of the inner resistive layer out to a distance sufficiently removed from the tearing 
surface that parallel thermal conduction once again dictates the energetics in a force-free field configuration. 
The oscillating temperature in the new, thermally dominated modes produces small-amplitude (relative to the 
temperature) velocity oscillations. The magnetic field perturbation, however, remains unaltered from its behavior 
without conductive heat transport and is essentially identical to that obtained in the usual tearing-mode analysis. 
Subject headings: hydromagnetics — instabilities — Sun: corona 


1. INTRODUCTION 

Sheared magnetic fields, in which one component of the field 
reverses direction, are susceptible to development of the 
resistive tearing instability. This has been established with 
incompressible, resistive magnetohydrodynamic (MHD) 
theory, which considers the dynamic magnetic field-velocity 
interaction (Furth, Killeen, and Rosenbluth 1963). In applying 
this theory to solar coronal plasmas, optically thin radiative 
energy loss may strongly influence the resistively driven 
reconnection. For instance, it has been shown that the same 
sheared, force-free, magnetic field configurations which are 
amenable to resistive tearing also create a unique spatial 
location for the occurrence of a radiatively driven thermal 
instability (Field 1965; Chiuderi and Van Hoven 1979). 
Resistivity (and, therefore, reconnection) was omitted from 
these compressible MHD analyses, the essential physics being 
provided by the velocity-energy (temperature) interaction. 

More recently, the competing effects of field reconnection 
and radiative loss were combined into a unified theory by Van 
Hoven. Steinolfson. and Tachi (1983) and Steinolfson (1983). 
Two unstable, purely growing modes resulted from this study. 
The faster radiative mode is predominantly radiation driven 
and has a growth rate typical of the uncoupled condensation 
mode of the thermal instability. The mainly reconnection- 
driven tearing mode grows like the usual resistive tearing 
instability. Both, however, differ from their respective un- 
coupled modes, and each has characteristics of the other. This 
has been demonstrated by Steinolfson and Van Hoven (1984) 
who show that the radiative instability involves a level of 
reconnecting fields comparable to that in the tearing branch. 

Thermal conduction has generally been neglected in theories 
similar to those discussed above — especially when the field is 
force-free. Parallel conduction was omitted on the assumption 
that the force-free field would thermally isolate x-points and 
magnetic islands (also assuming insignificant energy transport 
to external boundaries parallel to the plane in which reconnec- 
tion occurs ). Parallel thermal conductivity exceeds perpendicu- 


lar conductivity by approximately 12 orders of magnitude in 
the solar atmosphere; hence, in regions where they compete, 
one would expect parallel transport to dominate (Chiuderi and 
Van Hoven 1979). These contributions to the thermal conduc- 
tion become comparable when the cross-field temperature 
gradient has a scale of 10 cm (for the physical conditions 
used later in this paper). This is much smaller than the resistive 
tearing layer width (Furth, Killeen, and Rosenbluth 1963) for 
the same conditions, which is 10 3 cm. Based on this argument, 
the perpendicular component should be inconsequential in 
instabilities involving reconnection since the tearing-layer 
width usually dictates the smallest spatial scale. However, in a 
force-free magnetic field, parallel conduction near the field- 
reversal location (y = 0 in the slab geometry discussed later) 
does not directly contend with the perpendicular component 
since it is effectively nullified by the large (infinite) distance to 
the boundaries. As a result, perpendicular conduction then 
competes with the other prominent energy transport 
mechanism— radiative loss. These two become commensurate 
for temperature gradients over a distance of 1.8 x 10 4 cm. 
which is larger than the tearing-layer width discussed earlier, 
thereby implying that perpendicular heat flux cannot be 
neglected. At some larger distance from the tearing layer, the 
usually dominant parallel conduction again governs the evolu- 
tion. Thus, contrary to usual expectations, the force-free field 
configuration creates an environment, near the field reversal, 
where perpendicular thermal conduction may influence energy 
transport. 

Both parallel and perpendicular thermal conduction are 
considered in this paper, in addition to finite resistivity and 
radiative loss. After presenting the equations and the model in 
the following section, the spatial singularity that arises with 
consideration of only the parallel heat conduction (no 
perpendicular component) is discussed in § III. The removal 
of this singularity and the formation of temperature oscillations 
(thermal ripples) by inclusion of the perpendicular heat-flux 
component are demonstrated in § IV. Some elaboration on the 




Resonant absorption of phase-mixed AlfvSn surface waves in ideal 
and resistive magnetohydrodynamics: initiai-vaiue problem 
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Numerical solutions of the linearized magnetohydrodynamic equations are used to investigate 
resonant energy absorption in the continuous frequency spectrum. Energy in the phase-mixed 
surface waves relatively dissipates with the absorption time and width scaling as resistivity to the 
— 1/3 and 1/6 powers, respectively. 


The ph enome n on of Alfven surface waves in a m e d i u m 
with noounifbrm Alfven velocity is well known and has been 
studied extensively for some time (see, e.g., Hasegawa and 
Uberoi 1 and references therein). Early interest in the surface 
wave centered on low-frequency heating of laboratory plas- 
mas. 2 The highly structured nature of the solar atmosphere 
subsequently generated interest in these waves as a candidate 
for the elusive coronal-heating mechanism. 2 

Energy transfer from the surface waves to the plasma is 
generally considered to involve phase-mixing. As back- 
ground. it is worth noting that surface waves propagating on 
a discontinuity are not damped and that the medhtm. in this 
case, sup por ts a discrete sp e ctrum of modes. In add iti o n. 
distinct namni modes on a continuous magnetic field or 
density profile must oscillate indefinitely (imaginary fre- 
quency) in the ideal m ag ncs o hydr ody nam ic (MHD) limit. 4 
Surfaco-wave damping (resonant ab sorption) occurs as a 
consequence of the generation of a continuous frequency 
spectrum by disturbances in a non uniform, but continuous, 
plasms. The normal modes of the continuum phase mix and 
eventually result in s reduction of wave amplitude. The 
damping is centered at the spatial location (resonance) where 
the continuum frequency matches the local Alfven velocity. 

R e s on ant absorption in the continuous spe c trum has 
been treated analytically with the linearized equations of 
ideal MHD. 1 The solution "" "«■»» of noncollective oscilla- 
tory m odes of the continuous sp ec trum and collective oscil- 
latory modes that are exponentially damped (the resonantly 
absorbed modes). The further addition of resistivity to this 
model does not alter the afaeorption rate of ideal MHD in the 
limit of small resistivity.* Although the wave decay can be 
computed with ideal MHD, the actual energy absorption 
must occur by nonideal pro cesses such as the con v ersi on of 
surface waves to kinetic Alfven waves 7 or by fluid dissipative 
processes (e.g., resistivity). 

In the present letter we consider an initial-value ap- 
proach to the study of Alfven surface waves, in which the 
linearized MHD equations are solved numerically in time 
and sp ac e. A disturbance of a particular wavelength is excit- 
ed in a non uniform plasma, and its temporal evolution is 
s i mu la ted with numerical solutions of the ideal and, sepa- 
rately, resistive MHD equations in geome tr y. This pro- 
cedure generates a continuous spec trum therefore, sim- 

ulates phase-mixing and subsequent resonant absorption. 
Resistive d is sipa ti on provides the dominant energy absorp- 


tion mechanism in this study, rather than mode conversion 
to kinetic Alfven waves. 

A slab geometry is used with an ambient magnetic field 
directed along the x axis and non uniform in the x direction. 
All quantities are assumed independent of the third dimen- 
sion, and the initial density is taken to be uniform. The lin- 
earized, incompressible momentum and induction equations 
can then be written as [in dimensionless form and for a dis- 
crete disturbance wavenumber in the x direction, exp(ikr)] 

+ aFV'B. (1) 

“ .« ofb.+i-V 2 ^. (2) 

where V 2 —^Vcbt 2 — a 2 . The x components, u, and B, , are 
given by the solmoidal equations. The quantity F represents 
the initial magnrtic field, which is assumed to be uniform 
outside a region of width 2a centered at x « 0. All dattnees 
an referenced to this distance a. The initial field dec reases 
monoton ically from the magnetic field reference value B 0 
(F«* l)forx< — I to a smaller value B^ \F = B^ /B 0 ) for 
x> 1. In the nonumform region F vanes as a fifth-order po- 
lynomial of odd powers with zero derivatives at x = ±1. 
Time is referenced to the Alfven time r» and velocities to a/ 
r*. 

The tiro parameters in Eqs. (1) and (2) are the wave- 
length parameter a «■ Iva/X (waven um ber k = 2v/X ) and 
the Lundquat number S ™ r,/r», when r, is the resistive 
time. The two times an given by r* ™ «(4«p) in /2? 0 . 
r, — 49 bVci 7, when p, e, and if an the uniform density, 
light speed, and uniform resistivity, respectively. 

Equations (1) and (2) an solved numerically using a fin- 
lte-difference method that has been successfully applied in 
the study of resistive tearing m odes.* The time step must be 
small enough to follow the temporal oscillations. A non uni- 
form grid is necessary to resolve the spatial oscillations in the 
phase-mixed waves and to compute their asymptotic (large 
|x|) behavior with a reasonable total number of grid points. 

A velocity disturbance of the form v, (x) s exp( — a|x|) 
initiates the surface waves. This particular disturbance was 
selected since it excites just the surface wave, while numer- 
ous other attempted forms (e.g., a broader or s uniform dis- 
turbance) also trigger transverse Alfven wsves in the uni- 
form medium. The boundary conditions require that the 
perturbed variables in the time-dependent computation have 
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A cylindrical axisymmetric tearing mode model for solar flares is investigated numerically. Large 
magnetic energy release only occurs when there are at least two mode rational surfaces in the 
current-carrying plasma. 


I. INTRODUCTION 

It has been proposed that solar flares may be produced 
by reconnection of magnetic flux loops. 1 However, it has 
been recently pointed out that the finite length of magnetic 
arch is strongly stabilizing, 2 precluding all but sausage-type 
(m = 0) tearing modes. In this paper we investigate the non- 
linear development of m ** 0 tearing modes in the cylindrical 
solar flare model of Molt and Van Hoven. 1 We find that for 
substantial release of magnetic energy it is necessary to have 
more than one mode rational surface within the current- 
carrying plasma. 

In a study of linear m =» 0 tearing stability. Molt and 
Van Hoven found that fast reconnection with the appropri- 
ate timescale for solar flares occurred when a mode rational 
surface was located at the edge of the current channel. The 
nonlinear development of this fast r ec o nn ec ti o n is quite dif- 
ferent, depending on the number of mode rational surfaces 
present If there is only one rational surface, the released 
magnetic energy is less than 1 % of the magnetic energy of 
the initial equilibrium. When there are two mode rational 
surfaces, the equilibrium current is larger and much more 
free energy is available during reconnection. The magnetic 
energy released is a large fraction, about 23%, of the energy 
stored in the initial equilibrium, as in the planar force-free 
case. 4 

As pan of the nonlinear development shon -wave- 
length magnetic islands coalesce to form islands with the 
longest available wavelength. This suggests that the final re- 
connected state depends primarily on the magnetic forces 
and is not sensitive to the form and strength of the resistivity. 

II. EQUATIONS OF MOTION AND UNEAR STABILITY 

The solar flare equilibrium model of Mok and Van Ho- 
ven is a cylindrical plasma of length L, with a force-free 
magnetic field. Within the current channel of radius r t , the 
field is given by the well-known Bessel function model of 
Lundquist. 1 For radii greater than r , , the plasma current is 
assumed to vanish. In our computations, we have caused the 
current to fall off smoothly to zero in a layer of width 5. We 
have also assumed the presence of a fixed conducting shell of 
radius r„ >r t , for numerical reasons, rather than allowing 
the fields to extend to infinity. 

We have studied the linear and nonlinear, axially sym- 
metric, resistive modes associated with this model using the 
incompressible equations of motion** 7 in cylindrical coordi- 
nates r, 6, z, with resistivity ij and viscosity p. 


BmVV'xVS +IV8, 

(1) 

f.VI/XW + KV0, 

(2) 

— (Vv~ 2 pVU) = V(r- 2 d •<f')X'V'P-V9 
dt 


+r- 4 /-|/-pr- 4 firi : Ar 
vZ dz 

- 

+ $ VpxViP-V6 + pr~*A *(ri 

•CM. (3) 

dt 

(4) 

— (/>- J ) - V(/>“ 1 )xVlP-Vff 4- rjr-'A •/. 
dt 

(5) 

pl-K-VI XVIP-V0 + pA *K. 
dt 

(6) 

where 


dt dt 

(M 

d \ d P 
A • *r — — — + — . 
dr r dr d* 

(8) 


The equations may be marie dimensionless by normalizing 
lengths to r„ V to Br*, and time t o p xn r~ ' B ~ '. A key 
parameter for tearing modes is the ratio of resistive decay 
time to the shear Alfven period, which we will denote as 

J-J.tOM-'V. (9) 

The energy is the sum of kinetic energy KE and poten- 
tial energy PE, with 

KE - jp JJ dz dr r~ *[(Vf/ f + K *], 

?£ = i- jjdidrr-'W'Tf + I*), 

j-E- -^JJdrdrr-'Uri-f'P-HV/) 1 ] 

-pjjdzdrr-'WU) 1 + {?/:)*). ( 10 ) 

The energy is conserved, except for dissipation due to 
resistivity and viscosity. The integration extends over the 
plasma volume rOo, 0<zC L. 
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ABSTRACT 

He report oo the progreae of e search for precursors chat have direct physical connections 
to the start of subsequent solar flares. The discussion includes recant results at radio, 
visible, ultraviolet, and x-ray wavelengths, which are relevant to the pro-impulsive 
(onset) phase. He also relate the aspects of a theoretical scenario, based on magnetic 
recoeaeccioo with transport-coefficient phase changes, for explaining flare onset. The 
pertinent tine scales for pre-iapulaiva teapots! developments are discussed. 

ihtkodu crion 

Solar activity encompasses a raoga of tiaaaeaias. Overall levels follow the veil-known 
eleven-year activity cycle, and large-scale aagaetic features evolve on a tiaascale of 
aootba. The life cycle of an individual active region ranges aver eeny days to a few 
weeks. Particularly in its growth and nature phase, eagnatic flux is often sene to change 
oo a tiaascale of bours. This flux evolution plays an important role in the 
destabilisation of the quasi-static aagaetic fields in the corona, and the subsequent 
release of energy on a tiaascale of seconds to tens of seconds in the fora of solar flares. 

Since the 1970's chare has been a growing appreciation of an into res diets tiaascale of 
active-region changes, those changes occurring on a scale of aiautes to perhaps tana of 
aiaucea iaawdiataly preceding the start of the impulsive flare energy release. This aspect 
of solar activity, sometimes denoted as the “onset" phase of a solar flare, is worthy of 
study for a nuabar of reasons, aaong which is the insight it can give into the plasaa 
instabilities which result in the impulsive energy release. 

The pre-SMI status of our knowledge of solar flares end active regions has been reviewed in 
the eooographa which suaaarised the results of the Sky lab Solar Horkahops / 1,2/. However, 
one unavoidable weakness of Che f ila-liaicad Sky lab observations was the paucity of data oo 
preflare cooditioas. In part to remedy this, the flare Buildup Study was organised as part 
of the Solar Maximum Year to help coordinate the observational coverage of preflare 
conditions. One aspect of the flare Buildup Study/Solar Maximum Analysis program was 
concerned with flare Precursors and Onset (fPO). This paper, which suamMrises the results 
of the fPO aspect of the program, is baaed on the deliberations of a study-group (0) 
meeting at Big Bear Solar Observatory in 1983, and the subsequent efforts and reports by 
its aenbera. 

This paper is not intended to be a comprehensive review of what is known about flare 
precursors and onset. Rather it is intended to highlight some of tbs recent 
developments, for a review of the observational status prior to the present solar auxinum 
the reader in referred to reviews by Van Hoven at al. / 3/ and Martin Ikl. Another curreot 
discussion will appear in the proceedings of the 1983-84 Solar Maximum Mission workshops. 

Osiog the terminology suggested by Martin /4/, this group restricted its deliberations to 
the "distinct" category of preflare events, namely those for which direct physical 
associatioo with the subsequent flare is implied. Studies such as this are, of course, 
bedeviled by the issue of what is preflare and what is flare. Our choice has been to 
consider chose active-region phenomena which occur oo a timescale significantly longer than 
the impulsive-phase energy release, significantly shorter than the overall flux-evolutioo 
process (which is dealt with elsewhere in this issue), and which have a plausible, direct, 
physical connection with the subsequent energy release. Under these rules the "flare" 
begins with the start of the rapid, impulsive-phase, energy release. Put another way, we 
are concerned with the preflash phase of the flare (or equivalently the onset phase in 
Sturrock'a /!/ terminology). 
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RADIATIVE AND RECONNECTION INSTABILITIES: 
COMPRESSIBLE AND VISCOUS EFFECTS 


T. TaCHI. R. S. STEINOLFSON. and G Van HOVEN 
Department of Physics. Unntrstty of California. Irrine. C.1 92717. V. S i 

(Received 4 July; in revised form 16 September. 1984) 

Abstract. Filaments and Hares are prominent indicators of the magnetic fields of solar activity. These 
instability phenomena arise from the influence of weak transport effects (radiation and resistivity, 
respectively) on coronal magnetodynamics and energy flow. We have previously shown that the filament 
and flare i tearing or reconnection ) mechanisms are resistively coupled in sheared magnetic fields of the kind 
existing in active regions. The present paper expands this treatment to indude the effects of compressibility 
and viscosity, which are most prominent at short wavelengths. The results show that compressibility affects 
the radiative mode, including a modest increase of its growth rate, and that viscosity modifies the tearing 
mode, partially through a decrease of its growth rate. A comprehensive discussion of the mode structures 
and flows is presented. The strongest effect found is a reversal, at very long wavelengths, of the radiative 
cooling of the resistive interior layer of the tearing mode, caused by compressions! heating 


1. Introduction 

Filaments and flares are observed to occur in sheared magnetic field regions in the solar 
atmosphere (Martin, 1973 ; Van Hovea et al., 1980). The former are primarily driven by 
unstable radiation (Field, 1963) and the latter by finite resistivity (Furth et al.. 1963). 
Considerable progress has been made on the quantitative description of these mecha- 
nisms : the filamentation (thermal) instability has been treated in a nonuniform magnetic 
field (Chiuderi and Van Hovea, 1979), and the tearing (reconnection) instability has 
been extended to extreme astrophysical conditions (Steinolfson and Van Hoven, 1983). 

Recently, a unified formulation of the interacting radiative and reconnection processes 
has been achieved (Van Hoven et al ., 1983 ; Steinolfson, 1983 ; Van Hoven et al.. 1984), 
which appbes in the incompressible and inviscid limit. These works established the 
existence of two unstable modes which coexist in a force-free solar coronal environment. 
They also examined the structure of the individual modes in terms of the perturbed 
physical variables and of the derived ratios between their amplitudes and widths. The 
latter quantities led to a d efinit ion of the physical characteristics which clearly distinguish 
the two modes. 

The radiative mode resembles the thermal instability of ideal MHD. with a growth 
rate similar to that previously found by Chiuderi and Van Hoven (1979) in the 
compressible limit. This results from the fact that the resistive effect on the energetics 
of the plasma is small and, as a result, the inertia term is balanced by the unstable 
radiation, so that the growth rate is found to be directly proportional to the amount of 
radiative cooling available. The only significant difference between the ideal and resistive 
modes is the varying degree of magnetic reconnection in the radiative mode ( Steinolfson 
and Van Hoven, 1984), which is a strong function of the wavenumber. In short, the 
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RESISTIVE WAVE DISSIPATION ON MAGNETIC INHOMOGENEITIES: NORMAL MODES 

AND PHASE MIXING 

R. S. Steinolfson 

Department of Physics. University of California. Irvine 
Received 1984 July 20: accepted 1985 February 13 

ABSTRACT 

Numerical solutions of the linearized, resistive MHD equations indicate the presence of two characteristic 
forms of wave decay on magnetic inhomogeneities. Solutions for relatively long-wavelength disturbances and 
large values of the resistivity have the behavior of decaying normal modes. In the opposite limits, phase 
mixing becomes important, with the eventual build-up of large spatial gradients in which the resistive terms 
dominate. This composite behavior is shown to be conceptually consistent with analytic results, which predict 
that the complete solution consists of the sum of collective and noncollective contributions. The numerical 
simulations go beyond analytic theory by defining where each contribution prevails and computing wave 
decay when phase mixing is effective. The dispersion relation in the normal-mode regime is also determined 
analytically, with some approximations, and is in good agreement with the numerical solution. The decay time 
of the normal modes varies with Lundquist number as S 1/6 , while the phase-mixed decay time scales as S I/3 . 
Subject headings: hydromagnetics — Sun: corona 


I. INTRODUCTION 

Of the numerous methods that have been proposed for 
hejitng the solar corona (for reviews see Kuperus, lonson, and 
Speer 1981: Wentzel 1981), dissipation of magnetohydro- 
rirumte tMHDl waves propagating on magnetic inhomoge- 
*tties remains one of the more promising candidates. 
Application of the theory to coronal loops suggests that the 
*i'es may supply enough energy to attain the required heat 
“rut t lonson 1978: Heyvaerts and Priest 1983; Sakurai and 
Ltjmk 1984: Hollweg and Sterling 1984). 

!n general, wave-dissipation studies divide into two separate 
theories depending on the wave polarization. Most attention 
An locused on disturbances lying in (and propagating in) the 
of the initial field and inhomogeneity. Heyvaerts and 
^**1 (1983). however, considered waves with a displacement 
P er I*ndicular to both the initial field and inhomogeneity 
k^ar Alfven waves). The former (referred to here as plane- 
prized waves) may decay by virtue of an imaginary com- 
in the frequency spectrum, while the latter lose energy 
®*°ugh a phase-mixing process. 

riie actual energy-transfer mechanism is poorly understood 
•present, but it has been suggested that it occurs by collision- 
Processes arising through the excitation of kinetic Alfven 
*A'es (Hasegawa and Chen 1976) or a Kolmogoroff turbulent 
®rodc to high wavenumbers (Heyvaerts and Priest 1983). The 
j** a| on taken in the present paper is that the dissipation 
u ‘ ts from collisional MHD dissipative processes (e.g.. 


''usiiv 


'<>')• Shear Alfven waves are not considered, although 
mixing will be shown to be important in planar polariz- 
*' , °n as well. 

, v 11 C;ln he shown, using the energy or variational principle. 
,.' JI 'hstinct normal modes on a continuous magnetic field or 
Profile must oscillate indefinitely in the ideal MHD 
^ " 'Haieinan 197-Si. Despite this well-known result, an appar- 
Vom Ue ^ ecav ' n ‘he above case can be obtained analytically 
^,'heal MHD by assuming that the frequency (for a real 
'•auniheri contains both real (oscillating) and imaginary 


(decaying) components (Tataronis and Grossman 1973). We 
attempt to shed some light on these seemingly contradictory 
results by showing, first of all. that the imaginary frequency 
from ideal theory agrees well with that computed numerically 
using resistive MHD. To illustrate that this agreement may not 
be fortuitous, the addition of resistivity to the analytic, ideal 
theory is shown to produce a small correction (depending lin- 
early on resistivity) to the nonresistive solution. The picture 
thereby produced for this normal-mode behavior (and sup- 
ported by numerical computations) is that the dissipative terms 
are indeed small at each instant throughout the disturbance, 
relative to the other terms, yet they produce a cumulative elfect 
resulting in wave decay. 

In addition to decaying normal modes, the resistive solution 
also contains a regime where a phase-mixing process domin- 
ates. This disparity of responses is interpreted to be analogous 
to the behavior predicted analytically from ideal MHD. The 
ideal solution consists of noncollective oscillatory modes (at 
the local Alfven frequency) coupled with collective oscillatory 
modes (at an average Alfven frequency), which are exponen- 
tially damped at a uniform rate (Sedlacek 1971: Tataronis 
1975). The amplitude of the noncollective oscillations decreases 
inversely with time, which will be shown later to be significant 
in comparing these results with the numerical computations, 
including resistivity, presented herein. 

One could also envision that the phase-mixing process in 
ideal MHD occurs when normal modes of the continuum 
phase mix (Grad 1969). Similarly, the imaginary frequency or 
“damping rates" computed without dissipation could be inter- 
preted as coupling rates of these normal-mode motions to the 
continua. 

Two separate approaches to the study of wave decay on 
inhomogeneities are considered in this paper. First of all. ideal 
and resistive MHD are used to analytically derive the collec- 
tive (normal-mode) behavior, referred to above, for long- 
wavelength disturbances. The resistive MHD equations are 
then solved numerically in time and space usine an initial- 
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PROPAGATION OF ENERGETIC ELECTRON STREAMS 
IN SOLAR FLARES 


YUNG MOK 

Physics Department, Vmwersuy of California. In in*. CA 9271 7. U.S.A. 

(Received 3 February; in Anal form 4 Sept em be r . 1984) 

Abstract. The microscopic stability of an electron stream flowing down to the photosphere from the corona 
is examined. It is found that, while a power-law distribution is stable in the low-density corona, it is unstable 
against the generation of magnetized electron plasma waves in the high-density photosphere. The scattering 
of these energetic electrons may alter their radiation signatures. 

1. Introduction 

The generation of energetic electrons, streaming down to the photosphere along a 
magnetic loop, has been a major feature of many flare models. Although the acceleration 
mechanism is an open question, observational evidence from X-ray measurements 
seems to indicate the existence of an electron population with a hard spectrum, leading 
to thick-target radiation (Jordan, 1981 ; and references therein). In addition, the energy 
deposition by these electron streams (beams) has been studied by various authors using 
different energy spectra, mostly power-law distributions corresponding to observational 
data (Ricchiaza and Canfield, 1983; Brown etaL, 1978; Lin and Hudson, 1976; 
Shmeleva and Syrovatskh, 1973; Brown, 1973). However, the question of the propa- 
gation of these energetic electrons also deserves careful attention. As a group of energetic 
electrons streams through a plasma, various interesting phenomena may occur, and 
eventually affect their propagation characteristics, and their ability to penetrate deep into 
the photosphere. At this time, it is unclear how such effects will change the electron’s 
ultimate radiation signature. 

In the present work, we will study how well such an electron stream can propagate 
through the solar atmosphere along a magnetic loop. In particular, we will examine their 
capability for exciting collective instabilities which can significantly change their velocity 
distribution, giving rise to possibly different radiation signatures. In Section 2, we will 
study a particular instability which, we believe, is most likely to be excited by an electron 
stream in this environment. In Section 3, we use a special distribution function model 
to illustrate this possible scenario. In Section 4, some consequences of this instability 
are discussed. 

2. Microscopic Instability 

A common feature among all electron-stream (beam) flare models is an anisotropic 
velocity distribution in the low atmosphere, since an isotropic spectrum cannot cany 
a net flow. Generally speaking, anisotropic distributions resulting from injection of 
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A UNIFIED TREATMENT OF THE FILAMENT AND FLARE INSTABILITIES* 


Gerard Van Hoven 
Department of Physics 
University of California 
Irvine, California 92717 (U.S.A.) 


ABSTRACT 

Filaments and flares occur in sheared magnetic 
structures as a result of radiative cooling and resistive 
reconnection, respectively. A new integrated theory of 
these two unstable processes is described, which includes 
the relevant effects of magnetohydrodynamics and energy 
transport. The normally dissociated thermal and tearing 
phenomena are coupled together by a temperature-dependent 
Coulomb resistivity. As a result, the fi lamentation and 
flaring instabilities of a sheared field may coexist, as 
is familiar from the solar example. 

The growth rates and spatial structures of these two 
modes are detailed here. The much faster radiative 
instability is shown to provide significant magnetic 
reconnection, particularly at shorter wavelengths. The 
long-wavelength reconnection mode is found to be abetted by 
the resistivity increase caused by the dominance of cooling 
at the X point, in contrast to its nonradiative behavior. 
Implications of these results for the development of coronal 
activity are described. 


INTRODUCTION 

It is well known, as exemplified by the development of 
solar activity, that increasing mag ne t i c- f i e Id shear or 
stress gives rise to the formation of filaments (Chiuderi 
and Van Hoven, 1979) and to flares (Van Hoven, 1979). 
Theories of the former mechanism (Field, 1965), which is 
driven by a radiation ouput that decreases with temperature, 
have usually ignored the resistive magneto- hydrodynamic 
effects of the resulting, very collisional, relatively low- 
temperature plasma. Treatments of the latter 

Invited paper, to appear in. the Proceedings of International 
Astronomical Union Symposium 107, Unstable Current Systems 
and Plasma Instabilities in Astrophysics (Reldel) . 


NONLINEAR EVOLUTION OF THE RESISTIVE TEARING MODE 


R. S. Steinolfson and G. Van Hoven 
Department of Physics 
University of California 
Irvine, California 92717 


ABSTRACT 


Numerical solutions of the MHD equations are used to investigate 
the nonlinear behavior of the tearing instability. The mode evolves 
from a linearly growing excitation, followed by a period of greatly re- 
duced nonlinear growth. Constant- Y solutions evolve much more slowly 
than comparable nonconstant- ? modes with orders of magnitude less con- 
version of the stored magnetic energy. The nonconstant- Y computations 
indicate a reduction by approximately 20.X of the energy in the initial 
shear layer. For long-wavelength solutions, secondary-flow vortices, 
opposite in direction to the linear vortices, generate a new magnetic 
island centered at the initial x-point. 


COMPUTATIONAL PROCEDURE 

The nonlinear phase of the tearing mode (Furth et al., 1963) is 
studied in slab geometry using incompressible, constant-resistivity, MHD 
theory. The initially stationary plasma, with uniform thermodynamic 
properties, is embedded in ■» force-free, nondissipating magnetic field. 

A linear mode, at its maximum linear growth, provides the initial state 
for the nonlinear computation. We present results for a magnetic 
Reynolds number S (ratio of the resistive time to the hydromagnetic time) 
of 10 4 and values of the wavelength parameter a(2rra/X , where a is the 
shear scale and X the disturbance wavelength) of 0.05, 0.13, and 0.50. 

[A larger parameter range and additional computational results are con- 
sidered by Steinolfson and Van Hoven (1983b).] The mode with a = 0.5 
is a constant- 1 }' solution in the linear regime, while the other two are 
nonconstant- Y , and the or = 0.13 mode corresponds to maximum linear 
growth (Steinolfson and Van Hoven, 1983a). 

The linear theory predicts a chain of x-points and islands in the 
magnetic field lying along the tearing surface (x-axis in our geometry) 
at y = 0. We isolate one wavelength of this initial disturbance and do 
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Resistive decay of Alfven waves in 
a non-uniform plasma 

By Y. MOK and G. EINAUDIt 

Department of Physics, University of California, Irvine CA. 9271? 

(Received 16 October 1984) 

The effect of resistive dissipation on the propagation of an MHD disturbance in a 
non-uniform plasma is examined. The present analysis, baaed on a boundary- 
layer technique, shows the existence of resistive normal modes with complex 
eigenfrequencies. The real part of the eigenfrequency is associated with an 
oscillatory behaviour and defines the location in space of the layer where 
resistivity is important. The dissipation mechanism is responsible for the 
damping of the wave, in contrast with previous works in which the ideal MHD 
theory was used. 


1. Introduction 

The properties of magnetohydrodynamic (MHD) waves propagating in a 
non-uniform plsama have received much attention in the past several years (see 
Haaegawa ft Uberoi 1982, and references therein). This interest arises from the 
fact that the Alfvdn wave is a potential candidate for radio-frequency heating of 
a fusion plasma and that it can also be responsible, in the highly structured solar 
atmosphere, for the heating of the magnetised corona. 

In the presence of a time-constant, spatially non-uniform, background 
magnetic field, the infinite-conductivity hydromagnetic wave equation admits 
solutions which do not have a well-defined spectrum of discrete modes. Thus, no 
dispersion relation analogous to the unif orm case can be found. The reason for 
this behaviour is that the local variation of the Alfv6n-wave velocity gives rise 
to a singularity in the equations at the point where the wave-phase velocity equals 
the local Alfvdn velocity. The presenoe of such singularities leads to a continuous 
spectrum of Alfv6n-wave frequencies, with associated eigenfunctions exhibiting 
logarithmic behaviour at the singular point. Well-behaved solutions can be 
obtained only by integrating the singular modes over the entire spectrum. In this 
way, it is possible to obtain the asymptotic behaviour of non-collective oscillatory 
modes of the continuous spectrum, which present a (dissipationless) damping 
proportional to the inverse power of time (Banton 1964). Using an initial-value 
approach (Laplace transform), rather than a normal mode approach, it is 
possible to show that the solution also exhibits collective modes of oscillation 
(the surface eigenmodes) with position-independent frequency and exponential 
damping (resonant absorption) (Sedla£ek 1971). The damping rate found in this 
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THE PHYSICS OF THERMAL INSTABILITY IN TWO 
DIMENSIONS 

L. SPARKS andG. VAN HOVEN 
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(Received 10 January, 1985) 

Abstract. Previous studies of a thermal (radiative) instability in a sheared magnetic field have shown that, 
under solar coronal conditions, cool condensations can form in a small neighborhood about the shear layer. 
Such results have served to model the formation of solar filaments (or prominences (observed to occur above 
photosphenc magnetic polarity-inversion lines. A surprising conclusion of these studies is that the width 
of the condensation does not depend on the thermal conductivity ( *, ). By examining the mass-flow patterns 
of two-dimensional condensations in the absence of thermal conduction, we demonstrate that local plasma 
dynamics and the constraints imposed by boundary conditions are together sufficient to explain the size 
of the condensation width. In addition we present the results of a senes of numerical calculations which 
illustrate the characteristic mode structure of shcared-field condensations. 

1. Introduction 

Solar prominences, interstellar clouds, and condensations in planerary nebulae are 
examples of astronomical objects which owe their existence to a condensation process 
of nongravitational origin. It has been proposed (Parker, 1953; Field, 1965) that such 
phenomena occur as a consequence of an instability in the thermal equilibrium of a 
diffuse medium. The condensation mechanism relies on optically thin radiation whose 
dependence on thermodynamic variables (e.g., density and temperature) is such that a 
cool, dense perturbation loses more energy through radiation than it gains through 
adiabatic and non-adiabatic heating processes and thermal conduction. As it cools, the 
pressure decreases, and the resulting mass inflow increases the density further. 

Cool condensations in the solar corona (called filaments when viewed against the disk 
and prominences when viewed above the limb) are often observed to form above a 
magnetic neutral (polarity-inversion) line in regions of increasing magnetic field shear 
as indicated by photosphenc magnetograms (Martin, 1973; Leroy. 1978). Local heat 
conduction, which is sharply collimated in a direction parallel to the magnetic field, 
would dominate radiation and other forms of energy transport if no field were present, 
thereby suppressing the thermal instability. Thus, one can expect the equilibrium field 
structure to exert a strong influence over the formation of prominences. 

These empirical considerations have motivated previous computational studies of the 
dynamics of the thermal instability in a sheared magnetic field (Chiuderi and Van 
Hoven, 1979; Van Hoven and Mok, 1984). A surprising result of these studies is that 
the width of the condensation in a direction perpendicular to the shear layer does not 
depend on the thermal conductivity ( k„ ) but only on the wavenumber and on the growth 
and radiation rates. Contrary to what might be expected, the condensation width does 
not correspond to those points at which the radiation loss is roughly balanced by parallel 
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Resistive Alfven normal modes in a non-uniform plasma 

By G. EINAUDIf and Y.MOK 

Department of Physics, University of California, Irvine, California 92717 
(Received 8 May 1985) 

Resistive normal -mode solutions of the MHD equations are found numerically 
in a smooth, non-uniform, magnetic field. The (a, 3) boundary within which 
normal-mode solutions exist is explicitly computed, where a is the normalized 
wavenumber and 3 the Lundquist number. As an extension of our previous 
analytic results (Mok & Einaudi 1985), the damping rate of these modes is 
computed to a higher accuracy, and is found to have an a + bS~l dependence, 
where a and 6 are independent of S. 


1. Introduction 

The dissipation of Alfvdn waves in a uniform plasma is generally believed to be 
weak. Recently, more attention has been turned to the damping of these waves 
in non-uniform configurations related to laboratory experiments and astro- 
physical conditions. The study of these waves has a wide range of applications. 
In particular, they may be responsible for the heating of solar coronal plasmas 
near active regions where the magnetic field is strongly inhomogeneous. Because 
they are easily generated in the low atmosphere by photospheric fluid motions, 
they are among the more plausible candidates for transporting energy to the 
corona. They may also be partly responsible for the heating of solar wind plasmas 
emerging from open magnetic regions where the field is far from uniform. In 
laboratory experiments, they are considered to be one of the potential candidates 
for heating fusion plasmas. 

Studies of Alfvdn waves in inhomogeneous media have been well documented 
(Haaegawa & Uberoi 1982, and references theirein). In the ideal magneto- 
hydrodynamic (MHD) framework, these waves form a continuous spectrum, in 
contrast to the discrete normal modes which appear in the case of homogeneous 
media (Baraton 1964; Sedlacek 1971). For the latter, small perturbations can 
oscillate at their local Alfvfin frequencies and propagate at their local Alfven 
velocities. Unfortunately, these normal-mode solutions usually possess logarith- 
mic singularities in non-uniform media, which appear at the positions where 
their eigenfrequencies match the local Alfvln frequencies. In order to obtain a 
mope nearly complete physical picture, it is necessary to take into account 
other phenomena, such as resistivity, viscosity or kinetic effects, which become 
important near the singularity of the solution. 

t Permanent address: Scuola Normals Superiors, 56100 Pisa, Italy. UCI Technical 
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VISCOUS NORMAL MODES ON CORONAL INHOMOGENEITIES AND THEIR ROLE AS A 

HEATING MECHANISM 
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ABSTRACT 

Viscous damping of Alfven surface waves is examined both analytically and numerically using incompress- 
ible MHD. Normal modes are shown to exist on discontinuous as well as continuously varying interfaces in 
Alfven speed. The waves experience negligible decay below the transition zone. High-frequency waves damp 
just above the transition region, while those of lower frequency lose energy further out. A comparison of dissi- 
pative decay rates shows that wave damping by viscosity proceeds approximately two orders of magnitude 
faster than by resistivity. 

Subject headings: hvdromagnetics — Sun: corona 


1. INTRODUCTION 

The generation of magnetohydrodvnamic (MHD) waves by 
photospheric and convection zone motions and the wave 
transfer of this mechanical energy into the corona continue to 
be viable mechanisms in some coronal heating theories (see, 
e.g.. reviews by Holiweg 1983: Priest 1982: Kuperus. lonson, 
and Spicer 1981; Wentze! 1981). An important aspect of the 
heating problem, which is receiving some attention currently, 
concerns the process by which the waves dissipate. 

Waves may. of course, propagate throughout the corona, 
but those whose existence is intimately associated with inho- 
mogeneities in the local Alfven velocity appear to offer some 
heating advantages. Suggested processes by which the corona 
may absorb energy from these waves include the excitation of 
kinetic Alfven waves (Hasegawa and Chen 1976). phase mixing 
(Pritchett and Dawson 1978: Heyvaerts and Priest 1983; 
Browning and Priest 1984). a Kolmogoroff turbulent cascade 
(Holiweg and Sterling 1984: Nocera and Priest 1984), and 
resistive normal-mode decay (Steinolfson 1985; Mok and 
Einaudi 1985). Resonant absorption of such waves has been 
considered by. for example, lonson (1978). Pritchett and 
Canobbio (1981). Rae and Roberts (1981). and Sakurai and 
Granik (1984). 

The propagation of surface waves as normal modes on mag- 
netic or density discontinuities has been established using ideal 
MHD (e.g.. Rae and Roberts 1981: Hasegawa and Uberoi 
1982). When the discontinuity is replaced with a smooth 
profile, the ideal equations contain a singularity leading to the 
possibility of resonant absorption. A particularly lucid descrip- 
tion of the singular problem is given by Lee and Roberts 
(1986). The presence of dissipation, however, removes the sin- 
gularity. For instance, it has been shown that decaying normal- 
mode solutions occur on a smooth mhomogeneity when a 
finite resistivity is introduced (Steinolfson 1985). 

The present paper considers the effect of viscosity on surface 
waves for both discontinuous and smooth variations in the 
Alfven speed. The viscous tensor given by Braginskii (1965) is 
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used. A related problem was studied by Gordon and Holiweg 
( 1983) for the discontinuous density case. They used the surface 
wave solution from ideal MHD and computed the potential 
heating rate due to wave fluctuations, while we self- 
consistentiy determine the viscous normal-mode behavior. In 
addition, their analysis only allowed velocity and magnetic 
field fluctuations perpendicular to the initial magnetic field. We 
include parallel fluctuations since it has been shown previously 
that they may be important in the dissipative process 
(Steinolfson 1984. 1985). Despite these differences in approach, 
the two studies are complementary with generally similar 
results and conclusions. 

After establishing the existence of decaying viscous normal 
modes, we compare the viscous dissipation rate with that due 
to resistivity. The viscous decay rate is always about two 
orders of magnitude faster than the comparable resistive rate. 

It. LINEAR EQUATIONS 

Single-fluid theory is used, and the corona is assumed to be 
incompressible and initially stationary. For a simplifying slab 
geometry, the ambient magnetic field is taken to coincide with 
the r-axis (perpendicular to the solar surface). Inhomogeneities, 
which in some cases we treat as a discontinuity, in either the 
magnetic field or density (i.e.. the Alfven speed) occur along the 
x -direction. Despite these variations in the initial quantities, we 
assume that the dissipation coefficients remain constant. 

Since we are mainly interested in outward-propagating dis- 
turbances. the wave-propagation vector is required to be paral- 
lel to the magnetic field. A perpendicular component of wave 
propagation would not be expected to produce any new physi- 
cal effects. For this assumed parallel orientation, there are no 
changes in the dependent variables along the y-axis. 

Of the five ion viscosity coefficients given by Braginskii 
(1965). one (v 0 ) is several orders of magnitude larger than the 
others for typical solar coronal conditions (v 0 = 0.61 g cm' 1 
s' 1 , v , = 2.8 x 10' 11 gem' 1 s' 1 , v, = 4v,.v 3 = 3.9 x 10' 6 g 
cm' 1 s' 1 , v* = 2v 3 for number density n = 10 9 cm' 3 . tem- 
perature T = 2 x 10* K. magnetic field B = 3.725 G). When 
the four smaller coefficients are neglected, a considerable sim- 
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THE GROWTH OF RADIATIVE FILAMENTATION 
MODES IN SHEARED MAGNETIC FIELDS 


Gerard Van Hoven 
Department of Physics 
University of California 
Irvine, CA 92717 


ABSTRACT 

Observations of prominences show them to require well-developed magnetic shear 
and to have complex small-scale structure. We show here that these features are 
reflected in the results of the theory of radiative condensation. We have studied, 
in particular, the Influence of the nominally negligible contributions of 
perpendicular (to B) thermal conduction. We find a large number of unstable modes, 
with closely spaced growth rates. Their scale widths across B show a wide range of 
longitudinal and transverse sizes, ranging from ouch larger than to ouch smaller 
than the magnetic shear scale, the latter characterization applying particularly in 
the direction of shear variation. 


INTRODUCTION 

Coronal prominences owe their existence to a condensation process which occurs 
as a consequence of an instability in the thermal equilibrium of a diffuse medium 
(Parker, 1953; Field, 1965; Hildner, 1974). The condensation mechanism relies on 
optically thin radiation whose dependence on thermodynamic variables (e.g., density 
and temperature) is such that a cool, dense perturbation loses more energy through 
radiation than it gains through adiabatic and non-adiabatic heating processes and 
thermal conduction. 

Prominences and filaments (as seen on the disk) in the solar atmosphere are 
often observed to form above a magnetic neutral (polarity-inversion) line in regions 
of Increasing magnetic field shear as indicated by photospheric magnetograms 
(Martin, 1973; Leroy, 1978). Local heat conduction, which is strongly attenuated in 
directions perpendicular to the magnetic field, would dominate radiation and other 
forms of energy transport if no field were present, thereby suppressing the thermal 
instability. Thus, one can expect the equilibrium field structure to exert a strong 
Influence over the formation of prominences. 

These empirical and physical considerations have motivated previous 
computational studies of the dynamics of the thermal instability in a sheared 
magnetic field (Chiuderl and Van Hoven, 1979; Van Hoven and Mok, 1984; Van Hoven et 
al., 1984; Sparks and Van Hoven, 1985). One result of these theoretical studies is 
that a sheared background field is necessary for the existence of a true localized 
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abstxact 

Wo describe recent progress in the study of fiara praeursors tad onset. Ha provide some 
nav thaoratical raaulta cm (ilaaaat format ion aod eruption, including especially tba 
dynamical af facta of magnetic fields. Us find aple evidence that energetic procaaaas ara 
alraady at work ia tba oaaat phase, a f aw miautaa bafora tba rapid riaa of tha hard- 
radiation iapulaa. la particular, tha prevalence of soft aad bard X-ray prahaatiag 
(to < 10 s X) ia datailad, along with a coaaactioa to tba launch af coronal aaaa 
a j action*. A possible iatarpratat ion of tba empirical time prof ilea osar onset is that 
"prahaatiag" signifies that tba flare has slowly started, aad tba transition to tba 
impulsive stage then represents a change of phase ia tba flare-iastability process. 
Finally, tha problai of tba iatarpratat ion of microwave preflare data ia delineated. This 
represents tba final report of tbia part of tha Flare Build-up Study. 


HfTBODOCTIOH 

By tba onset phase of a solar flare, we mean those phe nom en a that occur in tha faw minutes 
prior to tba impulsive phase aad have a direct, physical and perhaps irreversible 
connection to tha subsequent rapid conversion of magnetic energy. Tba existence of 
phenomena on such a timescale, iatarmadiata between tba bourn to days of active-region 
evolution aad tba aacoada to tana of seconds of the impulsive phase, has been established 
beyond question over tha last solar maximum. A wide variety of such manifestations has 
bean observed , ranging from X-ray aad WV emission, through filament activation and Ha 
brightening!, to microwave polarisation changes. 

It might be argued that studying tha onset of a flare is a waste of time, because, given 
the mreratable magnetic field, something is bound to trigger it. However, tha constrained 
rate of tha initial energy growth, as contrasted with the rapid subsequent release, 
indicates that a different regime of plasma physics is relevant, which can be observed 
without distractions from the overwhelming manifestations of the moin flare. Furthermore, 
a physical understanding of the onset phase is vital, if short-term flare prediction III is 
to become a reality. 

Tha extensive observations of onset-phase p he n o m en a during the past solar maximtai are a 
result of the application of tha observing philosophy that lay at tho heart of the Flare 
Buildup Study III', namely, one should observe active regions as continuously as possible 
until (not after) they flare. Tha preflare phase has therefore been tba subject of study 
of the Flare Precursor aad Onset (FPO) group of the Flare Buildup Study/ Solar Hasimun 
Analysis progran. Preliminary results were reported previously HI, This paper simmarinss 
the deliberations of tho FPO group at a meeting at Sac ranento Peak Observatory ia August 
1983, aod the observational and theoretical progress nade ia tha last two years. For more 
complete revieve of earlier work, tho reader is referred to tho papers by Vso Hoven /A/, 
Martin /5 / and Priest at al. /6/. 

We will organise our discussion of flare indicators on a temporal basis by beginning with 
early filament evolution aod ending with the late onset omission riso. 

F1LAMZHT P0IHATI0H 

One of the first manifestations of active-region development leading to s flare is an 
increaae of nagastic shear aod tho subsequent condensation of a filMsnt / 7/ . The cool 
filament appears above the polarity-inversion line, where the flare will later occur, and 
arises for a similar magnet ic-topology reason / 8, 31 . 
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IMPULSIVE PHASE SOLAR FLARE X-RAY POLARIMETRY 
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ABSTRACT 


The pioneering observational work in solar flare X-ray polarimetry was 
done in a series of satellite experiments by Tlndo and his collaborators in 
the Soviet Union; Initial results showed high levels of polarization in X- 
ray flares (up to 401), although of rather low statistioal significance, 
and these were generally Interpreted as evidence for strong beaming of 
suprathermal electrons in the flare energy release process. However, the 
results of the polarimeter flown by the Columbia Astrophysics Laboratory as 
part of the STS-3 payload on the Space Shuttle by contrast showed very low 
levels of polarization. The largest value - observed during the impulsive 
phase of a single event - was 3.4% + 2.2%. At the same time but 
independent of the observational work, Leach and Petrosian (1983) showed 
that the high levels of polarization in the Tindo results were difficult to 
understand theoretically, since the electron beam is isotropized on an 
energy loss timescale - an effect which substantially reduces the expected 
levels of polarization, although not to zero. A subsequent comparison by 
Leach, Emslie, and Petrosian (1985) of the impulsive phase STS-3 result and 
the above theoretical treatment shows that the former is consistent with 
several current models and that a factor of ~3 improvement in sensitivity 
is needed to distinguish properly among the possibilities. In addition, 
there is reason to expect stronger polarization effects at higher 
energies: There may have been a strong thermal component to the flare at 
the energies seen by the STS-3 instrument (predominantly below 10 keV), and 
in adddition the preponderance of Y-ray ( > 300 keV) events on the solar 
limb (Rieger et al. 1983) suggests that beaming must be important at 
sufficiently high energies. 
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NONLINEAR RADIATIVE CONDENSATION IN A SHEARED MAGNETIC FIELD 
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. ABSTRACT 

The thermal instability, in which plasma coaling and condensation axe driven by increasing radiation losses, is 
believed to underty a number of ftlamcntanon processes in astrophysics. Most often, however, the local region 
must be threaded by a sheared magnetic field, so that the strong temperature-equalizing influence of electron 
thermal conduction can be suppressed. This fact, which is well known from the solar example, leads to 
difficulties in treating the dynamics that accompany the filamentation energ etics. Previous linear calculations 
have this radiative instability and have shown the ««*«»«**» of magnetically channeled heat and mass 

flows. This Letter describes a unique. well-resolved, two-dimensional, nonlinear, numerical simulation of the 
formation of a filament in a force-free field. The condensation is initiated by a linearly unstable mode and 
widens until it is slowed by parallel (to B) thermal conduction. During the nonlinear evolution, the minimum 
temperature falls from 10* K to 10* K. and eventually r ea c h e s a state of local thermal equilibrium in about five 
e-foldiag times. 

Subject headings: hydromagnetics — magnetic fields — plssmas — radiation mechanisms — Sun: prominences 


1. INTRODUCTION 

Nongravitatiooal condensations, such as solar p rominences , 
owe their exi s t ence to the fact that thermal equilibrium in a 
plasma can be unstable in the presenoe of a radiation function 
that decreases strongly with temperature (Parker 1953). The 
fust complete dynamic tr e atm ent of this thermal instability in 
a uniform magnetized "w<i"m was given by Held (1963) who 
detailed the linear growth under various conditions. Nonlin- 
ear dynamic computations of radiative condensation in a 
uniform m ag ne tic field have been attempted by Raju (1968), 
using a one-dimensional model, and by Hildner (1974) in two 
dimen s i ons, although he was not able to carry bis computa- 
tions to local thermal M piiKhwim of b oundar y inter- 

actions. 

It is s crucial empirical fact, however, that filaments form 
in ntmuniform regions of magnetic polarity inverses (Martin 
1973). Such configurations serve to suppress the thermal 
conduction, which normally (in a high-temperature field-free 
region) provides s larger energy-flux contribution than radia- 
tion, by increasing the parallel (to B) to the 

boundaries. Chiuderi and Van Hoven (1979) analyzed fila- 
ment condensation in a sheared and produced the first 
localized two-dimensional unstable mode. This effort was 
extended by Van Hoven and Mok (1984) and Van Hoven, 
Sparks, sued Tachi (1986) to include the essential effects of 
perpendicular conduction. 

In this Letter, we report the first results of a two-and-one- 
half— dimensional (vector fields have com ponents in the 
ignorable direction) nonlinear simulatio n of unstable radiative 
condensation in a sheared, force-free, magnetic field, for 


conditions applicable to the solar c oro na. The treatment 
inchlrirs a reasonably complete energy-transport model and is 
therefore quite general and relevant to a number of as- 
tropbyacal situations, as was emphasized by Field (1965). 

D. PHYSICAL MODEL 

The situation of interest in the thermal instability can be 
described by the equations of ideal (infinite electrical conduc- 
tivity) magixetofaydrodynamics (Field 1965; Hildner 1974), 
relating the plana pressure p, density p, velocity e, and 
magnetic field B. To these are added the ideal gas law and a 
suitable tom tor the energy equation 

- -yfiV *c + (y - l)[ff 0 - P 2 *(T) + V • « • V7 - ], 

(1) 

where H 0 is an equilibrium heating function (which, in the 
absence of an accepted parametric model, is chosen to bal- 
ance the radiative losses at t “ 0), ♦(T) specifies the temper- 
ature variation of the radiation loss, and a is the anisotropic 
thermal conductivity which satisfies « A « s„ (Spitzer 1962). 
The T-dependence of the radiation law, due to optically thin 
nnitsioo from variably ionirari trace constituents, can be 
approximated by expressions of the form $(T) - RT r , where 
R and r are piecewise constant (r - 0.5, T > 10 7 K; r - 
- 1.0, 10* * < T < 10 7 ; r - - 2.5, 10 5 5 <T< 10 s *; r - 0. 
lO** < T < 10 5J ; r - 1.8, 10* 2 < T < 10* *; and r - 7.4, 
T < 10 42 [HQdner 1974]). The global equilibrium magnetic 
L91 
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ABSTRACT 

Dissipative Alfven waves in a nonuniform plasma are studied by using a normal mode analysis. Specific 
magnetic geometries are used to simulate various magnetic structures in the solar atmosphere. We have com- 
puted the real and imaginary parts of the eigenfrequency for each of the configurations and found that these 
eigenmodes can deposit a substantial amount of energy to the corona because of their short damping dis- 
tances. 

Subject headings: plasmas — Sun: corona 


I. INTRODUCTION 

The possiblity of heating the solar corona by the dissipation of Alfven waves has been explored by many authors. A number of 
aspects of such theories have been discussed in various review articles (e g., Kuperus, Ionson, and Spicer 1981). A wave theory for 
coronal heating must answer two fundamental questions. The first one is whether the energy requirement of the corona can be met 
by the observed wave amplitudes. The second one concerns the capability of the waves to dissipate a considerable part of their 
energy at the coronal level. 

As far as the first point is concerned, there is ample observational evidence indicating the existence of nonthermal motions in the 
chromosphere and coronal base. These motions, which are unresolved in space and time, are likely to be fluid motions or 
propagating hydromagnetic waves, or both. Although there are constraints on the amount of energy being carried by these waves, 
nonthermal motions of the order of 10-30 km s -1 have been observed (e.g., Bonnet 1978; Cheng, Doschek. and Feldman 1979; 
Doschek and Feldman 1977; Feldman, Doschek, and Tousey 1975). If these observed motions are present in a large volume of 
plasma at the coronal base, it seems that they can provide sufficient energy to heat the corona (Hollweg 1983). 

The great difficulty has been considered to be the low efficiency of viscous or ohmic dissipation, which is very weak in a uniform 
plasma and therefore unable to deposit a significant amount of energy within the distance scale of the corona. However, since the 
corona is far from being a uniform medium, it is important to know the properties of propagation and dissipation of Alfven waves in 
nonhomogeneous plasmas in order to properly understand the role of the waves in coronal heating. 

It is well known that no dispersion relation analogous to the uniform case can be found for a nonuniform plasma (e.g., Barston 
1964; Sedlacek 1971) if the effects of a dissipation mechanism are neglected. This is due to the singular behavior of the solutions of 
the ideal hydromagnetic equations in the neighborhood of resonant surfaces, where the phase velocity of the mode is equal to the 
local Alfven velocity. Modes with real frequencies, as in the framework of the ideal MHD equations (Lee 1980), exhibit a logarithmic 
singularity at this surface (Goedbloed 1983), making the number of boundary conditions to be imposed insufficient to determine the 
eigenvalue, thus indicating a continuous spectrum. The difficulty arising from this ideal MHD singularity can be circumvented by 
simply introducing a complex eigenfrequency, leading to a discrete spectrum of damped regular modes (Tataronis and Grossman 
1973). As pointed out by Lee and Roberts (1986), the decay rate of these modes must not be interpreted as a dissipation rate, because 
of the absence of true dissipation in the ideal MHD framework, but rather as a mode conversion of the collective disturbances into 
local oscillations within a layer around the resonant surface. One of the effects of the local oscillations is to increase the spatial 
gradients of the wave excitations, making it necessary to consider dissipative processes within the resonant layer. 

Some of the mathematical properties of the equations with resistivity added as a dissipation mechanism have been studied by 
several authors (Riedel 1986; Sakurai and Granik 1985; Pao and Kerner 1985; Dewar and Davies 1984; Heyvaerts and Priest 1983; 
Kappraff and Tataronis 1977). An explixit solution has been obtained by Mok and Einaudi (1985, hereafter Paper I), by introducing 
the electron resistivity or ion viscosity, or both, as true dissipation mechanisms. This solution has a mode-damping rate which is 
independent of the magnitude of the dissipation coefficients in first approximation. Interestingly, the damping rate is the same as the 
mode-conversion rate found in the ideal MHD framework, meaning that all the energy of the collective disturbances which is 
converted into local oscillations is dissipated. The property of the dissipative layer to be a perfect absorber of Alfven waves without 
reflection occurring has been shown by Bertin, Einaudi, and Pegoraro (1986). As a consequence of this local behavior of the solution 
around the resonant surface, the eigenfrequency of the mode depends only on the equilibrium parameters and on its wavelength, but 
not the detail of the dissipation mechanism. 

In this paper we consider several magnetic profiles that model different coronal regions and compute the real frequency as well as 
the damping rate of the Alfven waves traveling along the magnetic field. The possible rate of energy deposition by the waves into the 
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Radiative instabilities in a sheared magnetic field 
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The structure and growth rate of the radiative instability in a sheared magnetic field B have 
been calculated analytically using the Braginskii fluid equations. In a shear layer, temperature 
and density perturbations are linked by the propagation of sound waves parallel to the local 
magnetic field. As a consequence, density dumping or condensation plays an important role in 
driving the instability. Parallel thermal conduction localizes the mode to a narrow layer where 
k t = k'B/|B| is small and stabilizes short wavelengths k > k„ where k e depends on the local 
radiation and conduction rates. Thermal coupling to ions also limits the width of the unstable 
spectrum. It is shown that a broad spectrum of modes is typically unstable in tokamak edge 
plasmas and it is argued that this instability is sufficiently robust to drive the large-amplitude 
density fluctuations often measured there. 


I. INTRODUCTION 

Radiative cooling plays a fundamental role in the forma- 
tion of a number of astronomical objects, including solar 
prominences, interstellar clouds, and planetary nebulas. 1 - 1 
The literature on solar prominences is particularly exten- 
sive. 1 Prominences are regions of high density and low tem- 
perature that form in the solar corona. They can be up to a 
hundred times denser and a hundred times cooler than the - 
ambient background with scale sizes of lOMO* km. 

More recently, the importance of radiative cooling in 
the edge region of laboratory tokamak discharges has been 
recognized. The onset of high density disruptions, which 
limit the maximum density that can be confined in a dis- 
charge, is triggered by impurity radiation in the edge plas- 
ma. 4-4 A new phenomenon, the marfe, has also been docu- 
mented in a variety of tokamaks. 10-14 Somewhat below the 
density limit, a torot dally symmetric region of cold, high 
density plasma forms in the edge region on the small major 
radius side of the plasma. The peak density in the marfe can 
become comparable to that in the center of the main body of 
the plasma. 14 It has been proposed that the marfe, like the 
solar prominence, arises from a radiative condensation pro- 
em s l ’- 17 

Short scale length fluctuations in the edge density of 
tokamak plasmas have been measured for more than a dec- 
ade. 1 *" 11 The spectrum of these fluctuations is largely two- 
dimensional, perpendicular to the local magnetic field. Un- 
like the marfe, they therefore track the local magnetic field. 
The level of fluctuations can be substantial, of the order of 
the background density ( h/n — I ). 

It has been proposed that rippling modes, driven by the 
local gradient in parallel resistivity, are the source of these 
fluctuations. 11 However, low current in the plasma edge 
combined with diamagnetic propagation and parallel ther- 
mal conduction tends to weaken the rippling mode. 14 There- 
fore it does not appear that this instability is sufficiently ro- 
bust to drive fluctuations to the high levels often seen 
experimentally. 


Finally, there is recent experimental evidence that the 
amplitude of the fluctuations increases rapidly as the marfe 
density threshold is approached. We propose that the short- 
wavelength edge-density turbulence is also driven by the 
thermal condensation instability. 

The most primitive radiative thermal instability arises 
from the amplification of asimple flute (B-V = 0) perturba- 
tion of the local temperature: The growth rate is given by' 


r- 


(i) 


J_dL 
3 n dT 

where L(n,D — n*g( T) is the radiation rate. Ifa drop in the 
local temperature T increases the radiation rate (dL / 
dT <0), the temperature will fall even farther so that the 
system is unstable When the parallel wave vector 1, is non- 
zero the sound wave can flatten the pressure along the mag- 
netic field line and perturbations of the temperature and den- 
sity are coupled. The local growth rate of the thermal 
condensation instability is then given by 


r — L( 

T in \T dT 1 1 / 


( 2 ) 


where K t , is the co efficient of electros thermal conduction. 
The first term in Eq. (2) is always destabilizing and it arises 
from the density perturbation. Because of parallel pressure 
balance, a local decrease in the temperature causes an in- 
crease in the local density and thus an increase in the radi- 
ation rate. The growth of perturbations arising from the in- 
crease in the density (in contrast to the decrease in the 
temperature alone) is referred to as the thermal condensa- 
tion instability. Parallel thermal conduction is stabilizing. In 
an initially homogeneous plasma in a uniform magnetic 
field, the parallel wave vector can be made arbitrarily small 
and instability is always possible as long as dL /dT <1L/T. 

In configurations of interest in both laboratory and as- 
trophysical applications, however, the geometry of the sys- 
tem is often not so simple. There is observational evidence 
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AN MHD SIMULATION OF PLASMA FLOW PAST 10: ALFVEN AND SLOW MODE PERTURBATIONS 
Jon A. Linker 1,2 , Margaret G. Kivelson 1 2 and Raymond J. Walker 
Institute of Geophysics and Planetary Physics, University of California, Los Angeles 


Abstract. We have studied the flow of plasma past Io 
using a time dependent, three-dimensional magnetohy- 
drodynamic (MHD) simulation. In addition to observing 
the “Alfven wing”, a standing Alfven wave perturbation 
expected from analytic theory, we have found that the 
other MHD modes contribute important perturbations. 
In particular, standing slow mode perturbations also are 
present in the flow. 

Introduction 

The interaction of Io, Jupiter’s innermost Galilean sat- 
ellite, with the Io torus plasma, has long been a subject 
of interest. Much of the discussion has been formulated 
in terms of perturbations carried by Alfven waves. It was 
Drell et al. [1965], in studying drag effects on satellites in 
the earth’s ionosphere, who first demonstrated that the 
presence of a conductor in a flowing plasma results in a 
standing Alfven perturbation in the rest frame of the con- 
ductor. An Alfven wave is launched when a conducting 
obstacle is placed in motion relative to a magnetised fluid. 
For the Alfven wave in a uniform plasma, the group veloc- 
ity is along or opposite to the magnetic field in the plasma 
rest frame. In the conductor’s rest frame, the Alfven per- 
turbation is also carried downstream by the background 
flow. The leading edge of the perturbation region occurs 
along the Alfven characteristic, at an angle tan~ l M A to 
the background magnetic field, where M A is the Alfven 
Mach number. This region, where the perturbation cur- 
rents flow, is known as the “Alfven wing”. 

Subsequent work has applied the concept of an Alfven 
wave interaction to the flow of the torus plasma past Io 
[e.g., Gocrtz, 1980; Nevbaucr, 1980; Southwood et al., 
1980; Wright and Southwood, 1987]. Although these pa- 
pers give a good description of the Alfven perturbations 
far from Io, they consider only the Alfven mode and do 
not model perturbations in the flow close to Io. Wolf- 
Gladrow et al. [1987] modeled numerically the magnetic 
field and current perturbations close to Io. However, they 
neglected the slow mode and other finite plasma beta ef- 
fects as well. 

In this paper we present results from a three-dimension- 
al MHD simulation of the flow of the torus plasma past 
Io. MHD simulations have been used previously to model 
the interaction of the solar wind with the earth’s magneto- 
sphere [e.g., Lebocuf et al., 1981; Wu et al., 1981; Brecht et 
al., 1982; Ogino, 1986; Fedder and Lyon, 1987] and with 
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comets [e.g. Schmidt and Wegmann, 1982; Ogino et al., 
1988]. In our simulations we solve the time-dependent, 
compressible MHD equations, so that all of the MHD 
wave modes are included. We have found that in addition 
to the Alfven mode, other MHD modes can contribute im- 
portant perturbations. In this paper we describe the slow 
mode perturbations caused by the interaction of Io with 
the torus plasma. 

Description of the Simulation 

The details of the numerical techniques used in our sim- 
ulation are discussed by Linker [1987] and will be pre- 
sented in a forthcoming paper; here we briefly outline the 
simulation model. We solve the following form of the nor- 
malized MHD equations: 

^ + V.(pv) = 0 (1) 

= (*) 

! 0*0 + V. (pirn +(/> + ^) 

= -i-V-W (3) 

^-Vx (vxb) = -Vx (Rj) (4) 

In the above equations, p is the density, v is the veloc- 
ity, P the plasma pressure, B the magnetic field, and 
J = V x B is the current density. W is a viscous stress 
tensor used for numerical stability purposes. R, is the 
fluid Reynolds number, R is the resistivity, and 7 = 5/3. 
For the results shown in this paper, R, = 50 and Rm, 
the magnetic Reynolds number (^j^, where V A is the 
Alfven speed and the length scale L is an Io diameter) 
is 100. The dissipative terms are included to help damp 
out short wavelength ripples generated by numerical dis- 
persion, while leaving the longer wavelength phenomena 
mini mall y affected. We have performed simulations with 
magnetic and fluid Reynolds numbers in the range of 10- 
200. The results remain qualitatively the same in this 
range, but higher levels of numerical fluctuations occurred 
in the higher Reynolds number cases. 

The equations are solved as an initial value problem in 
spherical coordinates by using a two-step Lax-Wendroff 
finite difference scheme. There are two boundaries in the 
simulation: at Io’s surface and at the outer boundary of 
the simulation region (r = 10 Io radii from the center of 
Io). Figure 1 shows the coordinate system for the simula- 
tion. Io is centered at r = 0, the center of the coordinate 
system. In our simulation, we assume that Io has no in- 
trinsic magnetic field. If this assumption is correct, in 
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resistivity and the tearing instability 
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Biskamp and Welter (1983) have defined an anomalous resistivity due to short- 
wavelength turbulence. They reported that this resistivity can be of either sign, 
and that negative anomalous resistivity in particular can affect the growth of 
the tearing instability. We use a spectral numerical-simulation code and 
ancillary diagnostics to analyse the behaviour of resistive magnetic tearing in 
the presence of turbulence of the sort postulated by Biskamp and Welter. We 
find that, in general, the ‘anomalous resistivity’ tends to return quickly 
towards zero even when artificially supported away from zero, and that its 
effect on tearing-mode behaviour is not consistent with its interpretation as a 
resistivity. We investigate analytically the behaviour reported by Biskamp and 
Welter, and the behaviour we observe. We also argue that, while not 
meaningful as a true resistivity, the ‘anomalous-resistivity’ parameter is a 
useful diagnostic showing the energy balance of the system - a property we 
refer to as Alfvenicity - illustrating, for example, the onset of nonlinearity in 
the tearing process. 


1. Introduction 

Magnetic reconnection, or tearing, is believed to be an important - perhaps 
the fundamental - mechanism for the release of magnetic energy to other forms 
of energy in both astronomical and laboratory plasmas (Priest 1985). A number 
of hydromagnetic simulations have been performed to investigate the linear 
and nonlinear dissipative evolution of reconnection under various initial 
conditions and other assumptions (Furth, Killeen & Rosenbluth 1963; Van 
Hoven & Cross 1973 ; Matthaeus & Montgomery 1981 ). The essence of magnetic 
tearing is the resistive diffusion of a sheared field toward a more energetically 
favourable configuration, allowing the conversion of magnetic energy to bulk 
motion and heat. 

Of significant interest is the efficiency of tearing : how fast will it proceed, and 
how completely ? This will certainly depend on the various parameters and 
conditions used ; a major objective of research into tearing has been to achieve 
theoretical rates and yields consistent with (primarily) solar-activity ob- 
servations (Van Hoven & Cross 1973). That there is no simple connection 
between the energetic requirement of magnetic release and the empirical output 
levels suggests the complexity of the processes involved. 

One concept that has been explored in the search for faster tearing is 
turbulence (Biskamp & Welter 1983 ; Matthaeus & Lamkin 1986). The premise 
here is that short-wavelength hydromagnetic excitations couple noniinearly 
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Abstract. The history of solar flare X-ray poiarimetry is reviewed and it is shown that as yet, there is no 
experimental evidence for such polarization. The present experimental limits are at the level of a few percent 
but these results may be biased by a large thermal component at low energies which may decrease the 
apparent polarization. To avoid this difficulty it will be necessary to make observations at higher energies 
where thermal emission is less important. 

The theoretical estimates of the polarization expected in the solar flare are also reviewed. The best present 
theoretical estimates are in the range of a few percent and are consistent with the present experimental limits. 

In this paper we discuss a new satellite instrument that has sufficient sensitivity at high energies to detect 
the polarization that is predicted by the present theories. The instrument sensitivity for a moderate (M class) 
event approaches polarization levels of 1 % in each of 7 energy bins spanning the 10 to 100 keV range for 
integration limes as short as 10 s. Comparable results can be obtained for an X class flare in 1 s. 


1. Solar Flare X-Ray Poiarimetry 

The idea that X-ray emission from solar flares might be linearly polarized and that 
polarization measurements could, therefore, provide a strong flare diagnostic was first 
discussed by Korchak (1967) and Elwert (1968). Subsequent theoretical investigations 
(Elwert and Haug, 1970, 1971; Haug, 1972; Brown, 1972; Henoux, 1975; Langer and 
Petrosian, 1977; Bai and Ramaty, 1978; Emslie and Brown, 1980) have resulted in 
polarization predictions for a variety of models. There are two extreme classes of models 
under investigation, termed ‘thermal’ and ‘non-thermal’, whose physical difference lies 
principally in whether the electrons responsible for the bremsstrahlung are part of a 
relaxed distribution or of a suprathermal tail. Although some form of hybrid model (e.g., 
Emslie and Vlahos, 1980) is probably appropriate for actual events, the basic compo- 
nents differ significantly in their polarization predictions: the thermal models predict 
polarizations of at most a few percent, due to either photospheric back-scatter of 
primary photons (Henoux, 1975), or an anisotropy in the source electron velocity 
distribution, caused by the presence of a field-aligned thermal conductive flux (Emslie 
and Brown, 1980). The beamed or linear bremsstrahlung models, on the other hand. 
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ABSTRACT 

It has been suggested that the solar corona is heated by the dissipation of electric currents. The low value of 
the resistivity requires the magnetic held to have structure at very small length scales if this mechanism is to 
work. In this paper we demonstrate that the coronal magnetic held acquires small-scale structure through the 
braiding produced by smooth, randomly phased, photospheric flows. The current density develops a filamen- 
tary structure and grows exponentially in time. Nonlinear processes in tbe ideal magnetohydrodynamic 
(MHD) equations produce a cascade effect, in which the structure introduced by the flow at large length scales 
is transferred to smaller scales. If this process continues down to the resistive dissipation length scale, it would 
provide an effective mechanism for coronal heating. 

Subject headings: hydromagnetics — Sun: corona 

I. INTRODUCTION 

The physical processes which heat the solar and other stellar 
coronae are not well understood (Priest 1982, pp. 206-234). 

However, it is widely accepted that the magnetic field in the 
solar atmosphere plays a prominent role in the energy- 
deposition process. Based on this concept, it has been sug- 
gested that the corona is heated through the dissipation of 
electric currents. In one version of this picture, the shearing 
and compression of a large-scale field whose footpoints are 
anchored below the surface can occur as a result of convective 
motions in the photosphere. This deformation induces electric 
currents to flow in the corona. The length scale on which the 
footpoints are shuffled by the convective motions corresponds 
to the granulation scale length (/„ ~ 10 6 m). If the coronal 
heating rate is to be accounted for by Joule dissipation of 
electric currents, the magnetic field is required to have struc- 
ture at scale lengths several orders of magnitude smaller than 
due to the large value of the conductivity in the high- 
temperature corona. Thus, in order to explain coronal heating 
by the above mechanism, it is required to demonstrate how the 
magnetic field can acquire structure at length scales l <8 
The proposed explanation of the mechanism by which the 
magnetic field acquires small-scale structure is somewhat con- 
troversial. Parker (1972, 1983a, 1986) has noted from a pertur- 
bation analysis that, in general, footpoint displacements lead 
to a loss of equilibrium, inducing current sheets to form in the 
corona. He has argued that the equilibrium MHD equations 
require the small-scale field to be invariant along the direction 
of the large-scale field. In situations in which the footpoint 
motions produce braiding of the magnetic field (which are 
expected to typify quasi-random photospheric flows), it is not 
possible to meet this requirement. According to Parker, the 
field would then necessarily require the presence of current 
sheets. The rapid resistive reconnection of magnetic field at 
these sheets would provide enhanced dissipation over that pro- 
duced by normal Jouie dissipation of the large-scale field, and 
would allow the field to reach equilibrium with a much simpler 
topology. This process has been called “topological 
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dissipation” (Parker 1972). In Parker’s hypothesis, then, the 
structure at l 4 l, in the corona is formed trivially through the 
creation of current sheets. 

In contrast, Rosner and Knobloch (1982) have noted that 
“the finite-amplitude behavior of a nonlinear system can be 
quite different from that predicted on the basis of small- 
parameter expansions.” Van Ballegooijen (1985, 1988a), Anti- 
ochos (1987), and Zweibel and Li (1987) have noted that 
well-behaved, continuous photospheric flows produce coronal 
fields which are free of discontinuities in the absence of initial 
neutral points in the field. However, van Ballegooijen (1985) 
hypothesized that random smooth flows (which are character- 
istic of the flows produced by subphotopheric convection) 
produce structure in the magnetic field at arbitrarily small 
length scales in a “cascade” process. In a statistical analysis 
van Ballegooijen (1986) determined the rate at which the struc- 
ture at /„ is expected to cascade to short length scales. For an 
idealized random flow profile van Ballegooijen (1988a) conjec- 
tured that the electric current density in the corona ought to 
build up exponentially in time. Indeed, a realization of the first 
few steps of this process (van Ballegooijen 19886), involving the 
numerical solution of the equilibrium MHD equations, indi- 
cates that such flows rapidly introduce fine structure in the 
coronal field. 

In this paper we describe dynamical calculations which 
address the issue of the creation of small-scale structure in 
magnetic fields. We solve the time-dependent three- 
dimensional ideal MHD equations in a geometry which is an 
idealization (Parker 1972) of the solar corona. The MHD 
equations are solved in the limit of vanishing gas pressure. We 
find that a sequence of smooth, randomly phased flows (like 
those used by van Ballegooijen 1988a. b) generates structure in 
the field at increasingly shorter lengths. This structure is evi- 
denced by an exponentially growing current density (with step 
number in the flow sequence) and a cascade of spatial structure 
to smaller scales. The rate of generation of the electric current 
density is consistent with van Ballegooijen's (1988a) predic- 
tions. These simulations are dynamical generalizations of van 
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ABSTRACT 

Solar coronal arches heated by turbulent ion-cyclotron waves may suffer significant cross-field transport by 
these waves. Nonlinear processes fix the wave-propagation speed at about a tenth of the ion thermal velocity, 
which seems sufficient to spread heat from a central core into a large cool surrounding cocoon. Waves heat 
cocoon ions both through classical ion-electron collisions and by turbulent stochastic ion motions. Plausible 
cocoon sizes set by wave damping are in roughly kilometers, although the wave-emitting core may be only 
100 m wide. Detailed study of nonlinear stabilization and energy-deposition rates predicts that nearby regions 
can heat to values intermediate between the roughly electron volt foot-point temperatures and the ~100 eV 
core, which is heated by anomalous Ohmic losses. A volume of 100 times the core volume may be affected. 
This qualitative result may solve a persistent problem with current-driven coronal heating; that it affects only 


small volumes and provides no way to produce the 
ments. 

Subject headings: hydromagnetics — Sun: corona — 
i. introduction 

Coronal heating models often invoke anomalous Joule 
heating by current-driven instabilities (Tucker 1973; Heyvaerts 
1974, Vlahos 1979; Hinata 1979). Of late, the electrostatic ion- 
cyclotron mode has emerged as a prime candidate for such a 
role, since it can be excited by lower electron-drift velocities 
and heats ions preferentially in the plane transverse to the 
magnetic field, lessening convective losses (Papadopalous 
1977). 

The electrostatic ion-cyclotron mode is particularly suited 
for cross-field transport because it results from a trade-off 
between harnessing the electron drift energy along the mag- 
netic field, versus the ion damping at the cyclotron frequency, 
Q,. Thus the waves store most of their energy transverse to 
B(k ± P k t ) and have frequency to & 1.1D,, just far enough from 
Q, to avoid ion damping. The resonance broadening model 
seems to describe at least qualitatively the nonlinear dynamics 
(Benford 1976) in which the linear dispersion relation still 
describes the observed wave properties. 

Benford (1983, hereafter Paper I) presented a detailed calcu- 
lation for a coronal loop heated by small-scale current- 
carrying filaments (~ 100 m in radius). He found that only the 
ion cyclotron mode appeared, and the system quickly settled 
down to a quiescent regime in which convective losses caused 
the instability to flicker on and off, “ percolating” the system to 
a steady temperature of ~ 100 eV. 

Benford’s calculation neglected transverse (to B) wave trans- 
port from the ion cyclotron modes, however. Here we make 
some estimates of the importance of this effect and argue that 
large annular regions transverse to a current-carrying zone 
. may be heated by the instability. We describe in detail the 
process of nonlinear wave transport and saturation in the sur- 
: rounding cocoon of initially cooler, roughly electron volt 
i plasma, which is then warmed to intermediate temperatures. 
This effect may be important in describing how coronal energy 
moves from the hypothesized small current-carrying zones. 

Although theory has found plausible ways to heat small 
regions to high temperatures, conventional thermal- 


extended warm structures perceptible to existing mstru- 
wave motions 

conduction mechanisms do not allow cross-field transport to 
large lateral volumes. Here we attempt to correct this defi- 
ciency, detailing the underlying plasma processes needed and 
exploring their ramifications. Once radiated from the small 
— 100 m core, a spaghetti-like current-carrying strand, the ion 
cyclotron waves are slowly damped in the surrounding cocoon, 
heating it. Damping by ion-electron collisions sets a limit on 
the cocoon size, R, exceeding 1 km. Moreover, stochastic ion 
motions in the wave fields can heat ions and electrons. This 
constraint leads to R comparable to 1 kilometer for plausible 
turbulence levels. Thus sizable cocoons are unavoidable. The 
true value of R must take into account convective cooling 
along the field lines as well. 

We consider these matters qualitatively, although with some 
quantitative treatment of energy transfer rates in the nonlinear 
stages. The aim of our work is to argue that such effects gener- 
ally make lateral transport an important mechanism. Other 
schemes, such as electron cyclotron emission are also possible 
(Melrose and Dulk 1984). 

We depend on the general approach of Paper 1, adding in 
§ II of this paper a qualitative sketch of how important and 
widespread the effect can be. Basically, we show how trans- 
verse wave loss sets an upper limit on the size, R, of any heated 
cocoon. This size is much larger than, say, the ion gyroradius 
because ordinary Coulomb processes are ineffective for solar 
coronal conditions. Section III explores in detail the under- 
lying microscopic behavior of the ion cyclotron instability 
which supports the transverse energy transport. We summarize 
our conclusions in § IV. 

ii. size of the lukewarm cocoon 

Consider a cylinder representing a straightened magnetic 
loop (Fig. 1). Currents flowing with drift speed u, exceeding the 
critical onset drift for ion cyclotron instability u c , heat a core to 
temperature T c , as described in Paper I. The system quickly 
reaches a “ percolating ” steady state characterized by marginal 
stability of the ion cyclotron modes. Convection down the field 
lines (z-axis) offsets wave deposition of energy into both ions 
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The activity of the Sun occurs in the solar 
atmosphere and is driven and confined by 
the Sun's magnetic field. The plasma atmo- 
sphere comprises the cool (S 10*K) and 
dense chromosphere, an intermediate transi- 
tion region, and the hot (T, 2 10*K) and dif- 
fuse (a,. 2 10* cm ’) corona; the solar atmo- 
sphere has been well observed from the ultra- 
violet through hard X rays, by the Skylab and 
Solar Maximum Mission spacecraft, among 
others. The atmospheric magnetic field has 
its source in the solar interior and is driven 
and energised by the global differential rota- 
tion and local turbulent motions at the visible 
photospheric surface; this field is mostly 
known from Zee man -effect and microwave 
measurements. The magnetic field lines re- 
turn to the surface in active regions and are 
open in coronal holes that provide the source 
for the solar wind. 

The primary concern of the Solar Plasma 
Theory Croup at the University of California. 
Irvine, involves the dynamics and energet ics 
of magnetic activity. The coronal field is 
strong (A 1 /2uoP >> I), nonuniform, and 
stressed by currents driven from the surface. 
This field provides a source of stored energy 
and an anisotropic medium for the channel- 
ing of mass and energy flows in the atmo- 
spheric plasma. 

There are a number of unsolved problems 
connected with observed solar activity. They 
are briefly stated (to be amplified in the fol- 
lowing sections) here. How is the corona 
heated: that is. how is the necessary energy 
transported up the adverse temper a ture gra- 
dient from the photospheric source surface 
and deposited in the corona, apparently roost 
effectively in strong-field regions? How does 
the atmospheric magnetic field channel and 
confine the ambient heat and mass flows so 
that the former are suppressed, allowing run- 
away radiation losses to occur, and the latter 
are abetted so that the cool condensation of a 
solar prominence is formed? How can the 
stressed magnetic field of an active region be 
reconnected quickly enough, in the highly 
conducting corona, to explain the short time 
scales of a solar flare? How are coronal mass 
ejections launched against the gravitational 
and magnetic forces of the tower atmosphere 
so that they can escape into the solar wind? 

The UCI group has made significant pro- 
gress on each of these problems, as will be 
described in the following sections. The key 
to our attack on the coupled dynamics and 
energetics of such nonuniform, anisotropic, 
nonlinear, active phenomena is the applica- 
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tion of large-scale numerical simulations, bal- 
anced and supported by analytic calculations 
and perturbation computations. 

Coronal Heating 

One of the outstanding mysteries of solar 
physics is how the corona is able to sustain its 
high temperature while the chromosphere 
below, separated only by a narrow transition 
region, has a temperature two orders of mag- 
nitude lower. Energy must be transported 
from the solar surface to the corona by mech- 
anisms other than thermal conduction, or 
large-scale convection, which is not present in 
the transition region. 

Because the corona is threaded by magnet- 
ic fields of various intensities. Alfven waves 
provide a promising medium of energy trans- 
port since they are naturally excited by pho- 
tospheric fluid motions driven by subsurface 
c onvection [Ionian, 1978]. However , the ener- 
gy-deposition rate of classic Alfven waves in a 
uniform magnetic field is known to be insuf- 
ficient by many orders of magnitude; thus we 
have explored the characteristics of these 
waves in various non uniform magnetic struc- 
tures b eli eved to be present in the solar atmo- 
sphere and have recalculated their energy-de- 
position rates. The effects of rlassiral dissipa- 
tion mechanisms such as resistivity [Met and 
Emandi. 1985; Emmudi and Mo*. 1985, 1987] 
are enhanced by several orders of magnitude 
by the field non uniformities existing in the 
corona. This is because fluid motions change 
sharply over a narrow resonant layer located 
in the magnetic inhomogenehy. The large 
magnetic-perturbation gradients, associated 
with the sharp velocity changes and main- 
tained by the wave motions outside the layer, 
induce a strong current that dissipat es the 
wave energy through resistivity. The dissipa- 
tion of kinetic energy by viscosity is also en- 
hanced in a similar manner. As a result, the 
energy deposition due to Alfven waves is now 
believed to be significant in the process of 
coronal heating. 

A second mechanism proposed for coronal 
heating [Pan krr, 1983] involves the generation 
of very narrow current layers (“sheets," in 
Parker's view) as a result of the convect i ve/ 
vortical motions of the surface foot poulu of 
the coronal field lines. We have performed a 
three-dimensional MHD simulation of this 
process [MiJde it aL, 1989], which evolves the 
coronal magnetic field from Parker's uniform 
initial state. Late-stage contours of the cur- 
rent density halfway along the field structure 
are shown in Figure 1 . A diagnosis of this 
randomly driven evolution demonstrates that 
the magnetic field gradients grow exponen- 
tially, rather than disconunuoualy. This result 


verifies an analytic prediction by van Balie- 
goaejen ( I988|. As the current densities in- 
crease because of the convective driver, their 
ohmic losses inj 1 ) also grow, until the limn at 
which the loul Poynting-fiux surface input is 
dissipated in the corona. 

We have also studied the problem of the 
cross-field transport of the coronal heat input 
from the typical narrow energy-deposiuon fil- 
aments (as in the example just described! to 
the wider atmosphere. One promising proc- 
ess involves ion-cyclotron turbulence, which 
has the lowest drift-current threshold (and 
thus is also a contributor to the primary heat- 
ing) and propagates substantially across the 
field. Amend! and Benford [1989] have shown 
that ion-cyclotron waves can effectively 
spread the heat input over 10 2 - 10 *x the ini- 
tial heat-input volume. 

Thermal and Filamentation 
Instabilities 

Prominences are relatively dense, gaseous 
structures extending many thousands of 
kilometers into the Sun's atmosphere. They 
typically form in a few hours and can endure 
for many weeks. Although prominences were 
obse r v e d in ancient times, the physical proc- 
ess responsible for their formation has re- 
mained a mystery. It has been proposed that 
prominences result from a radiative-conden- 
sation instability [Field, 1965] that occurs 
when heat flow is inhibited by solar magnetic 
fields [Cfcuden on d Van Haven, 1979). When 
the Adds are twisted (as indicated by in- __ 
creased fibril inclination), the dominant par- 
allel heat flow (kII/k. > 10‘°) is suppressed by 
the increased boundary distance along the 
field. The condensation instability is driven 
by temperature-dependent, optically thin ra- 
diation [« p 2 HT)] from highly ionized, heavy 
ion lines. When a dense cooling perturbation 
loses more energy through radiation [d&dT 
< 0, initially] than it gains through thermal 
conduction or adiabatic and nonadiabatic 
heating processes, the local temperature and 
pressure drop. This results in additional mass 
inflow along the magneuc field, leading to in- 
creased density and runaway cooling. 

The UCI group has been engaged in a 
long-term study of the linear and nonlinear 
condensational instability in the presence of 
magnetic shear. Linear studies (see Sparks and 
Van Haven, 1985. 1988. and references cited 
therein) have distinguished between dynamic 
condensations, modes whose spatial structure 
a determined primarily by force balance, and 
kinematic condensations, where consider- 
ations of energy balance (especially perpen- 
dicular conduction over small scales) are par- 
amount in setting characteristic scales of spa- 
tial variation. It is the latter, which do not 
exist when ■ 0, that are of greatest inter- 
est in explaining prominence formation. For 
such short- wavelength *.. v* 0 modes ( Drake et 
aL, 1988] plasma flow is primarily parallel to 
the magnetic field. Hence the field does not 
impede plasma compression, and significant 
density enhancement can occur. Kinematic 
modes also have nearly degenerate growth 
rates and quite small scales in both cross-field 
directions [Van Haven et al., 1986]. Thus they 
provide a natural explanation for the com- 
plex fine structure observed in prominences. 

To study the nonlinear evolution of the 
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We explore the long-term nonlinear evolution of a tearing-mode-unstable 
sheared-field plasma in a turbulent environment. Two different physical 
configurations are modeled, and a different computational system is used for 
each. Results of both sets of calculations show that magnetic tearing arises 
spontaneously provided that the initial turbulence energy level is below the 
natural saturation level of the tearing instability. 

We discuss briefly the relationship between our results and those of pre- 
vious calculations, concluding that there are no significant unexplainable 
disagreements. 


1. Introduction 

There has been considerable interest, in various communities, in turbulence 
and its relation to magnetic reconnection or tearing. Turbulence is perceived as 
a fundamental state of a wide variety of physical systems, with recent studies 
in ‘ chaos ' providing some insight into the structure and development, as well 
as the description, of turbulent behaviour. In short, we can expect just about 
any realistic complex system to exist in a turbulent environment ; thus we need 
to understand the consequences of this condition. From a more immediate point 
of view, many hydromagnetic phenomena in plasma physics are found in 
circumstances where turbulence is known to occur : for example, in tokamaks 
and other confinement devices, geomagnetic substorms, the solar corona (and 
presumably also the less observable regions of the sun). Quite reasonably, there 
have been efforts to understand puzzling effects, such as disruptions of 
confinement and anomalously fast solar-flare magnetic-energy release rates, in 
terms of turbulent reconnection and its consequences. This interest in turbulent 
effects on reconnection has been manifested in theory and in numerical 
simulation. Among the analytical work of interest in the present context has 
been that of Biskamp & Welter (1983), who claimed that under certain 
assumptions it could be shown that an anomalous resistivity and viscosity arose 
from excesses of magnetic or kinetic energy in the turbulent spectrum. 
Diamond el al. (1984) have also derived formal expressions for turbulentlv 
modified resistivity and viscosity, and Strauss (1986) has identified a hyper- 
resistivity, depending on the spectral composition of turbulence and having 
spatial extent, as well as an anomalous resistivity also dependent on the 
turbulence spectrum but effective only at the tearing layer. Working geometries 
are either slab or toroidal, the latter usually being approximated at some stage 
as locally cylindrical. 
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1. INTRODUCTION 

The temperature of the solar atmosphere rises from 6000 K at the 
photosphere to 1-2 million degrees in the corona. This paradoxical 
temperature increase implies the existence of a coronal heating 
source. The energy loss mechanism (principally conduction) in the 
corona requires a heating source of between 10® and 10 1 erg/cm^-sec, 
depending on the region considered, to sustain the corona. Several 
explanations have been advanced to account for this heat source, 
including the dissipation of waves incident from below the 
photosphere, and the direct dissipation of magnetic field energy. 
Detailed studies of wave dissipation have shown that it is difficult 
to heat the upper corona by this mechanism. He consider a mechanism 
in which the magnetic field acts as an Intermediary in converting 
mechanical energy present in photospheric convection into heat 
energy in the corona . 

The convective motions in the photosphere can twist and braid the 
large-scale solar magnetic field. This induces electric currents to 
flow in the corona. It has been suggested that resistive 
dissipation of these currents can heat the corona 1 . The 
photospheric motions have a characteristic length scale l v - 10 6 m. 
Since the high-temperature corona is highly conducting, the magnetic 
field needs to acquire fine-scale spatial structure (on the length 
scale l r - 10-100 m> to provide an effective heat source through 
resistive dissipation. 

Parker^' ^ has described a scenario in which footpoint motions 
produce a loss of equilibrium and lead to the formation of current 
sheets. Here we investigate how photospheric flows may induce fine- 
scale structure in the magnetic field as a result of a cascade of 
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THE EFFECT OF MASS LOADING ON THE TEMPERATURE OF A FLOWING PLASMA 
Jon A. Linker 1 ' 2 , Margaret G. Kivelson 1 and Raymond J. Walker 
Institute of Geophysics and Planetary Physics, University of California, Los Angeles 


Abstract. We have investigated how the addition of 
ions at rest (mass loading) affects the temperature of a 
flowing plasma in a magnetohydrodynamic (MHD) ap- 
proximation, using analytic theory and time dependent, 
three-dimensional MHD simulations of plasma flow past 
Io. The MHD equations show that the temperature can 
increase or decrease relative to the background, depend- 
ing on the local sonic Mach number (Ms) of the flow. For 
flows with Ms > -y/9/5 (when 7 = 5/3), mass loading 
increases the plasma temperature. However, the simula- 
tions show a non-linear response to the addition of mass. 
If the mass loading rate is large enough, the temperature 
increase may be smaller than expected, or the tempera- 
ture may actually decrease, because a large mass loading 
rate slows the flow and decreases the thermal energy of 
the newly created plasma. 

Introduction 

The term “mass loading” refers to the addition of mass 
to a flowing plasma by ionization of neutrals, lo, Jupiter’s 
innermost Galilean satellite, is embedded in a flowing 
plasma (the Io torus) and is also surrounded by a cloud 
of neutral atoms and molecules [Brown et al., 1983, and 
references therein]. In this paper we present results from 
three-dimensional MHD simulations of plasma flow past 
Io, including the effect of ionization of neutrals in the 
equations. Ionization may modify the plasma tempera- 
ture significantly. We demonstrate these effects by exam- 
ining the MHD equation for temperature and the results 
of simulation runs for three different values of Ms- 

Mass loading plays an important role not only in Io’s 
interaction with the plasma torus, but also in the solar 
wind interaction with comets and with Venus, in plasma 
flow past Saturn’s satellite Titan, and possibly in other 
plasma-satellite interactions. Although the simulation 
parameters we use in this paper correspond most closely 
to the Io torus parameter regime, our remarks regard- 
ing the plasma temperature (in an MHD approximation) 
apply to the plasma in these other systems as well. 


1987; McGrath and Johnson 1987]. The torus plasma 
flows at 57 km/s relative to Io, so on average, the neu- 
trals can be assumed to be at rest relative to Io. With 
that assumption, the normalized MHD equations can be 
written in the form: 

% + v-(pi r) = p. (i) 

|( pv) + V-(pvS+(p+~ > jl-BB > j=0 ( 2 ) 

,3, 

§f“Vx (vxB) = -Vx ( V jf) (4) 

where j> is the density, v the velocity, P the plasma pres- 
sure, B the magnetic field, 17 the resistivity, and J = V x B 
the current density, p, is the mass added per unit volume 
per unit time. We take 7 = 5/3. (Viscous terms are 
added to (2) and (3) in the simulation runs for numerical 
stability purposes, but they are irrelevant to the present 
discussion.) Equation (3) can be derived in a straightfor- 
ward manner, using (1) and an equation for the conserva- 
tion of energy (note that d/dt is a convective derivative) 
[Linker, 19871. Because we consider that ionization of 
neutrals results from electron-impact, there should also 
be a heat loss term in (3) representing the energy given 
up by the electrons. However, the flow energy gained by 
an oxygen ion created in the background flow (280 eV) 
is large compared to the energy lost by the electron (10- 
12 eV). The inclusion of the electron energy loss term did 
not alter the simulation results significantly; therefore, we 
have ignored the term. 

From (1-3) we can also derive equations for the plasma 
velocity and plasma pressure similar to those used by 
Schmidt and Wegmann [1982] and Ogino et al. [1988]: 


The MHD Equations and the Plasma Temperature 

We use the time-dependent MHD equations, including 
the effect of neutrals ionized near Io. When Jupiter’s 
gravity is accounted for, the escape velocity for neutrals 
from the surface of Io is about 2.3 km/s [Linker et al., 
1985] and the neutral density in Io’s exosphere is dom- 
inated by lower velocity neutrals on non-escape trajec- 
tories [Watson, 1981; Sieveka and Johnson, 1984; Linker 
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"I --*'(( ,,+ t ) 1 (5) 

^ + V.(Pv) = ( 7 -l)(-P(V-v) 

+ TjJ 2 + \psv' 1 ) (6) 


We see from (1), (5), and (6) that mass loading causes 
increases in the plasma density and pressure, and de- 
creases in the plasma velocity. In the context of the Io 
torus, it has been noted that mass loading can increase 
the plasma temperature [e.g., Goertz, 1980; Brown et al., 
1983; Shemansky, 1988]. However, if we assume an ideal 
gas law (P = pRT), (6) implies that the plasma temper- 
ature, T, satisfies 


dT 


pR 


V J 2 


- T(V ■ v) + ■— J - (7) 
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THE NONLINEAR EVOLUTION OF MAGNETIZED SOLAR FILAMENTS 
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ABSTRACT 

Thermal instability driven by optically thin radiation is believed to initiate the formation of plasma fila- 
ments in the solar corona. The fact that filaments are observed generally to separate regions of opposite, line- 
of-sight, magnetic polarity in the underlying, differentially routing photosphere suggests that filament forma- 
tion requires the presence of a highly sheared, local magnetic field. We have performed two-dimensional, non- 
linear, magnetohydrodynamic simulations of the local genesis and growth of solar filaments in a force-free, 
sheared, magnetic field. To clarify the essential physics of the nonlinear behavior, we have traced the evolution 
of generic perturbations possessing broad spatial profiles. We find that many of the physical considerations 
that govern the behavior of linear condensations play a similar role nonlinearly. Simulations of the evolution 
of initial random-noise perturbations produce filamentary plasma structures that exhibit densities and tem- 
peratures characteristic of observed solar filaments. Furthermore, in each of these simulations, the filament 
axis lies at a finite angle with respect to the local magnetic field, consistent with solar observations. 

Subject headings: hydro magnetics — Sun: magnetic fields — Sun: prominences 


L INTRODUCTION 

Observations of the solar disk in Ha often show thin, dark 
lines extending distances of 200,000 km or more. Such struc- 
tures are called solar filaments. They consist of cool, dense 
plasmas suspended in the solar atmosphere above the photo- 
sphere. When viewed above the limb, they appear as promi- 
nences. Magnetograms and associated Ha observations of 
fibrils reveal that solar filaments generally arise in regions 
where the local magnetic field is twisted and strongly varying. 
Observations have revealed that the axis of a fully formed 
prominence lies typically at a small angle (whose distribution is 
centered on 22! 5) with res pea to the local magnetic field 
(Leroy 19891 

The processes responsible for the formation of solar fila- 
ments are not well understood. It is widely believed that fila- 
mentary structures form initially as a consequence of a local, 
thermal instability in the corona (Parker 1953; Field 1965). 
This mechanism, designated “ the condensational instability," 
is driven by optically thin radiation. In thermal equilibrium, 
radiative energy losses are balanced by adiabatic and non- 
adiabatic heating and thermal conduction. Perturbing the 
plasma can upset this energy balance. When a dense plasma 
perturbation allows radiative cooling to dominate, the local 
pressure drops. This results in additional mass inflow and 
runaway cooling. 

Thermal conduction tends to inhibit the condensational 
instability (Orrall and Zirker 1961). Indeed, it is well known 
that heat flow (mediated largely by electron collisions) in an 
unmagnetized coronal plasma is so rapid that the conden- 
sational instability is completely suppressed for all temperature 
(and density) perturbations whose transverse scales are less 


1 Promt addros: Jet Propulsion Laboratory, California Institute of Tech- 
nology. Pasadena. California. 


than a limit somewhat larger than the size of a typical filament 
Consequently, the ability of a magnetic field to reduce 
(electron) heat flow by many orders of magnitude (in directions 
perpendicular to the field) is essential for the formation of 
observed coronal condensations. 

The presence of a magnetic field, however, affects the plasma 
dynamics as well as the energetics. Due to the high electrical 
conductivity and the low 0-value of the ambient corona, 
plasma fluid elements are tied to magnetic field lines, and mass 
flow tends to be directed parallel to the magnetic field — the 
direction in which heat flow is not affected by the field. How is 
it possible for a magnetic field to inhibit heat flow in a coronal 
plasma without simultaneously restricting the mass flow required 
for the growth of density condensations ? This is, perhaps, the 
central question facing attempts to attribute the initiation of 
fil ame nt formation to the condensational instability, and this is 
the principal question we wish to address in this paper. 

It has been suggested that filament formation may require 
the presence of a magnetic field B(x) that is locally sheared 
(Chiuderi and Van Hoven 1979). Coronal condensations with 
densities and spatial scales characteristic of solar filaments 
cannot form in a uniform magnetic field. Such a field prevents 
transverse mass flow from contributing to a local growth in 
density, while heat flow parallel to the field ensures that paral- 
lel mass flow will not give rise to condensational instability. 
Observational evidence supports the conclusion that shear in 
the magnetic field is a necessary local condition for filament 
formation. Magnetograms reveal that filaments typically 
separate regions of the Sun possessing opposite, line-of-sight, 
magnetic polarity, and Ha observations demonstrate that fila- 
ment formation is often initiated after a period of divergent 
fibril inclination (Martin 1973). 

The anisotropic influence of local, force-free magnetic shear 
on the linear dynamics of the condensation instability has been 
examined extensively in a series of studies (see Sparks and Van 
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MHD Simulations of Coronal Mass Ejections: 
Importance of the Driving Mechanism 

J. A. Linker and G. Van Hoven 
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We have investigated the importance of the form of the driving mechanism in MHD 
simulations of coronal mass ejections. Previous authors have performed simulations using 
a thermal driving mechanism, and have found that this mechanism cannot reproduce the 
observed features of mass ejections unless an elaborate model of the initial corona is used. 
We have devised a model simulation problem and have found that the use of a simple form 
for the initial corona, with an upward moving parcel of cold, dense plasma as the driving 
mechanism, can produce results that are consistent with many of the features observed 
by coronagraphs. Our results imply that the nature of the driving mechanism may play 
an important role in determining the dynamical evolution of mass ejections. 


Introduction 

Coronal mass ejections (CMEs) are important but 
poorly understood phenomena in solar terrestrial phy- 
sics. CMEs have been observed in space-based white- 
light coronagraphs since the early 1970s, yet fundamen- 
tal questions about their initiation and propagation in 
the solar corona remain unanswered. 

About 1/3 of the CMEs seen in the Skylab data, ami 
as many as 80% of those observed in the Solar Max- 
imum Mission (SMM) data below 3 R 3 (solar radii) 
can be classified as "looplike” [Wagner. 1984], Sime et 
al. [1984] identified three general characteristics of the 
looplike mass ejections observed with the Skylab coro- 
nagraph: (1) the sides of the loop (often referred to as 
the “legs” of the mass ejection) have a greater density 
enhancement than the top of the loop; (2) a depletion of 
density occurs between the legs of the mass ejection; (3) 
the legs of the mass ejection exhibit very little lateral 
motion throughout the time evolution of the mass ejec- 
tion. In this paper our references to CME observations 
are specifically to the looplike ejections of the form that 
Sime et al. describe. 

Coronal mass ejections are complex phenomena, and a 
realistic model must involve at least two spatial dimen- 
sions. time dependence, and magnetic forces. One ap- 
proach for studying such phenomena are time-dependent 
magnetohydrodynamic (MHD) simulations. The earli- 
est attempts at numerical simulation of CMEs explored 
the possibility that mass ejections could be initiated by 
the rapid injection of thermal energy at the base of the 
corona from a solar flare [e.g. Steinolfson et al.. 1978; 
Dryer et al.. 1979; VVu el al.. 1982]. In these simu- 
lations. the initial corona was modeled with a potential 
(current-free) magnetic field, and a pressure and density 
satisfying hydrostatic equilibrium. A pressure pulse was 
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introduced to model the thermal energy released in the 
flare, and the subsequent time evolution was followed. 
The resulting expanding fast-mode shock wave was iden- 
tified as the CME. However. Sime et al. [1984] pointed 
out that these simulations failed to replicate the three 
(previously listed) observed characteristics of looplike 
ejections. Because of these discrepancies. Sime et al. 
questioned the idea that CMEs could be regarded as 
compressional disturbances initiated by thermal energy 
release. We note that Dryer and Wu [1985] disputed 
this conclusion, as well as the generality of the CME 
characteristics identified by Sime et al. 

Independent of the observations of Sime et al.. another 
difficulty with the idea that CMEs are driven by thermal 
energy release from solar flares arises from the observa- 
tion that mass ejections are more often associated with 
eruptive prominences than with flares [Kahler, 1987. and 
references therein]. The eruption of a prominence and 
the deduced departure time of an associated CME have 
been found to be roughly simultaneous, whereas in flare- 
associated ejections, extrapolation of the observations 
suggests that the flare is initiated by conditions occur- 
ring after the launch of the mass ejection [Wagner et al., 
1983, Harrison. 1986]. 

As one response to the problems with the thermal en- 
ergy driving mechanism. Steinolfson [1988] and Steinolf- 
son and Hundhausen [1988] have taken an approach that 
emphasizes the role of the initial corona. They first ver- 
ified that simulations with a thermal driving mechanism 
and an initially hydrostatic, current-free corona with a 
closed magnetic t opology could not reproduce any of the 
observed features of mass ejections. They then went on 
to show that a thermal driving mechanism can simu- 
late many of the observed features of coronal mass ejec- 
tions. if one is able to specially tailor the initial corona. 
A streamer-like configuration was used, and a heating 
term was added to the energy equation to increase the 
background sound speed in the calculation. 

The work by Steinolfson and Hundhausen [1988] was 
a useful step in understanding CME evolution, in that 
they showed one way of qualitatively reproducing the 
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The ideal and resistive magnetohydrodynamic ( MHD) model is used to examine the dynamics 
and structure of the solar corona. When the coronal magnetic held is deformed by 
photosphenc flow it can evolve to states that become unstable to ideal MHD modes. The 
nonlinear evolution of these instabilities can lead to the generation of current sheets, held line 
reconnection, and energy release. The disruption of an arcade held and the kinking of coronal 
loops is described. The braiding of the large-scale coronal field by convective photosphenc 
mouons develops fine-scale structure in the magnetic field and leads to the development of 
intense current filaments. The resistive dissipation of these currents can provide an efficient 
coronal heating mechanism. 


I. INTRODUCTION 

The solar corona presents an environment that chal- 
lenges the ingenuity of the plasma physicist. The combina- 
tion of the slow evolution of magnetic fields in three-dimen- 
sional geometry and the occurrence of impulsive energetic 
events that are observed during solar flares makes the analy- 
sis of coronal plasmas difficult The hot highly conducting 
corona has widely disparate Alfven and resistive time scales: 
the Lundqum number 5 = r,/r A . where r A is the Alfven 
time and r M is the resistive diffusion time, is of order 10** . 
The simultaneous large separation of length scales (global 
field structures — 10* m. current filaments — 10* m) com- 
plicates the analysis further. A description of the corona 
based on heuristic, severely approximate qualitative models 
gives necessarily vague conclusions about its p rop e r ties. Re- 
cent developments in algorithms for the numerical solution 
of the magnetohydrodynamic (MHD) equations have en- 
abled us to follow the dynamical evolution of id e alized cor- 
onal magnetic fields in realistic geometry. Such MHD simu- 
lations have become indispensable in developing an 
understanding of coronal physics. We describe the applica- 
tion of three-dimensional, time-dependent, ideal and resis- 
tive MHD simulations to global coronal behavior. In partic- 
ular. we describe the disruption of an arcade field that is 
sheared by photosphenc flow, we investigate the nonlinear 
development of the kink instability in coronal loops, and we 
show how the braiding of the coronal magnetic field by 
smooth, random photosphenc flows can introduce fine-scale 
structure into the coronal field. The dissipation of the elec- 
tric current filaments so introduced can act as an effective 
coronal heating mechanism. The application of the MHD 
model to coronal magnetic fields is discussed by Parker, 1 
Priest.* and Bim and Schindler.* Also see a recent article by 
Low.* 

The crucial difference between coronal and laboratory 
plasmas is that field lines in the corona are rooted in the 
photosphere (line tied). Photosphenc motions cause cor- 
ona) field structures to evolve quasistaucally through se- 
quences of neighbonng equilibria. Cenain mouons can drive 
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the system across a stability threshold, leading to impulsive 
dynamical behavior. The practical simulation of this slow 
evolution, coupled with impulsive nonlinear behavior, re- 
quires the use of efficient numerical methods. The semi-im- 
plicit numerical method 1,4 is suitable for this problem. This 
technique allows a large time step to be used dunng the 
quasistatic evolution phase, which involves long-wavelength 
phenomena. The time step is reduced automaucally dunng 
the more dynamic phases of the nonlinear evolution ( cur- 
rent filament formation, reconnection, and front propaga- 
tion) for accuracy. 

The MHD model provides an accurate description of 
global coronal physics because coronal phenomena are char- 
acterized by targe length scales and long time scales (com- 
pared to the ion gyrondius and gyropenod. respectively). 
Kinetic models are still necessary, though, to describe the 
phenomena of x-ray and particle radiation that accompany 
flares. Although the plasma resistivity is generally negligibly 
small, it becomes crucially important in regions of concen- 
trated electric current density. Thus we use the resistive 
MHD model. 


TXB= (4ir/c)J, 

(1) 

VXE --JL2* 

c dt 

(2) 

E 4- ( l/c)vxB =3 » 7 J. 

(3) 

p ( ~ -t- v»Vv^ =5 ~-JxB — Vp + vpV : v. 

(4) 

+ T*(pv ) =0. 
dt 

(5) 

+ vVp = - ypV-v, 

(6) 

where B is the magnetic field. J is the electric 

current den- 


sity, ▼ is the plasma velocity, E is the electric field, p is the 
plasma pressure, p is the mass density, rj is the resistivity, v is 
the viscosity, and y = j. The ideal MHD model is obtained 
by setting 17 = 0. Line tying at the photosphere is implemen- 
ted by ustng an ideal Ohm's law at the boundaries corre- 
sponding to the photosphere, 

E 4- ( l/f)r,XB »0. (7) 
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PROMINENCE FORMATION IN A CORONAL LOOP 
Y. Mok. J. F. Drake. D. D. Schnack. and G. Van Hoven 
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ABSTRACT 

We present a model, which depends on the preferential deposition of heating in the legs of a coronal loop, 
that produces a stable prominence-scale condensation at the loop top. Dynamic stability is attained by the 
subsequent adjustment of local parallel gravity by a magnetic inversion at the loop (or arcade) apex. We 
describe a nonlinear numerical simulation of this process, which includes a deep chromosphere, a heating rate 
with a fixed dissipation length, and full solar gravity. 

Subject headings: hydrodynamics — Sun: atmosphere — Sun: corona — Sun: prominences 


L INTRODUCTION 

One of the unsolved problems of solar physics is the charac- 
ter of the energy-balance and mass-transfer mechanisms by 
which a prominence can cool and condense at the apex of a 
magnetic arcade or loop. The elements involved include suffi- 
cient local densificabon or suppression of parallel heat conduc- 
tion to allow radiative cooling to dominate the energy 
transport (Field 1963L and sufficient pressure gradients to pull 
material along the field from the chromospheric source 
(Hildner 1971), against solar gravity, to the coronal apex. One 
attempt to provide these conditions in a multistep process 
using a fixed gravitational well was described by Poland and 
Mariska ( 1986). A related, continuous, nonlinear process was 
recently detailed by Sparks. Van Hoven. and Schnack (1990), 
in the absence of gravity and with the chromosphere modeled 
by inflow boundaries. 

In this paper we describe a single-step or continuous mecha- 
nism for providing a cool condensation which is then confined 
in a gravitational well located at the apex of a coronal loop or 
arcade. 


IL MODEL AND EQUATIONS 


We begin by assuming that the system is confined by a 
strong magnetic field, so that the mass flow and heat flux are 
essentially one-dimensional. This assumption is valid in most 
coronal loops and magnetic arcades observed in the corona, 
which have field lines anchored in the photosphere. The evolu- 
tion of this system is described by the equations of hydrody- 
namics. supplemented by an energy equation to take into 
account the effects of radiative loss, plasma heating and paral- 
lel thermal conduction. The energy equation can be written as 




(1) 


where v is the fluid velocity, T the temperature, n the electron 
number density, H the heating function. n J #7") the average 
radiative loss rate per unit volume, y the ratio of specific heats. 
k Boltzmann's constant and id T) the (parallel) thermal con- 
ductivity. In our model *d — T >n ) is taken to be the classical 
value, although other nondassical models can easily be imple- 
mented in the present calculation. 


The radiation law <HT) is that given by Athay (1986). modi- 
fied by an imposed cut-off below 2 x I0 4 K due to Lya and 
other atomic lines which make the plasma optically thick at 
lower temperatures. The heating function is based on a model 
in which energy is supplied to the chromospheric and coronal 
plasma from the photosphere through various channels, such 
as Alfven waves, acoustic waves, etc Since the detailed heating 
mechanism and its characteristics are not yet fully understood, 
we will treat the present heating function in a phenomenologi- 
cal way as a working modeL Accepting the fact that photo- 
spheric fluid motions are the source of energy, we argue that 
this supply naturally decays with altitude as energy is depos- 
ited in the plasma along the upward path (see, for example 
Priest and Smith 1979). The rate of decay, however, depends on 
the poorly understood Jheating mechanism and is treated 
as a variable parameter in the present calculation. In the 
actual computation, the heating function has the form 
H 0 exp ( - x/rj, where H 0 is a constant, x is the distance from 
the base, and x M is the spatial decay-length of the heating 
source 

The magnetically collimated (1 — D) momentum equation is 
given by 



where p 3 nm p is the mass density, p ( — 2 nkT) is the pressure, 
and g(x) is the component of the gravitational force along the 
magnetic field. In the present calculation. g(x) contains two 
contributions. The first one has the form ~g 0 cos (itz/f), 
assuming that the confining magnetic field has a semicircular 
shape, where g 0 =» 2.74 x 10* cm s' 1 is the solar-surface 
gravity, and f is the total length of the loop. The second contri- 
bution is introduced to simulate the fact that the magnetic field 
will be distorted at the location where the condensation occurs 
(Kippenhahn and Schliiter 1937; Poland and Mariska 1986). 

This distortion of the projected gravitational force is 
modeled (as the condensation forms) by the quantity 
3 , sin 2n([x - /J/[/ - 2/J), with g x the amplitude of the dis- 
tortion in the range l f < x < 1 — centered near the top of the 
loop structure, where condensation is expected because the 
heating reaches a minimum. In particular, this term is not 
turned on until a substantial condensation (with density > 3 
times the nearby coronal value) has formed. This contrasts 
with the approach used by Poland and Mariska (1986X whose 
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DYNAMICAL EVOLUTION OF TWISTED MAGNETIC FLUX TUBES. I. EQUILIBRIUM AND 
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ABSTRACT 

The three-dimensional dynamical evolution of twisted magnetic flux tubes is studied using a time-dependent 
magnetohydrodynamic (MHD) model. The flux tubes are intended to model solar coronal loops, and include 
the stabilizing effect of photospheric line tying. The model permits the complete evolution of flux tubes to be 
followed self-consistently, including the formation, equilibrium, linear instability, and nonlinear behavior. 
Starting from an initial uniform background magnetic field, a twisted flux tube is created by the application of 
slow, localized photospheric vortex flows. The flux tube evolves quasi-statically through sequences of 
equilibria with increasing twist, until it becomes linearly unstable to an ideal MHD kink mode. In this paper 
the equilibrium properties and the linear stability behavior are discussed. The nonlinear behavior will be 
treated in the second paper of this series. The application of the method to the uniform-twist, Gold-Hoyle 
field confirms the previous stability threshold for kink instability and provides estimates of the resulting 
growth rate. When the method is applied to a loop with a localized twist profile, it is found that the loop 
becomes kink unstable when the twist exceeds 4.8n on axis. 

Subject headings: hydromagnetics — Sun: corona 


I. INTRODUCTION 

Twisted magnetic flux tubes can be detected in numerous 
observations of space plasmas. They are observed in the solar 
corona, and are likely to be present in the majority of plasma 
environments, including stellar coronae, planetary 
atmospheres, and the magnetosphere. It is believed that 
magnetic flux ropes influence the form and evolution of the 
Earth’s magnetotail. Loop-like magnetic structures are 
ubiquitous in EUV and X-ray images of active regions of the 
solar corona. They play a crucial role in^The structure, 
dynamics, and heating of the solar corona. In this paper we 
emphasize the behavior of solar coronal loops, although many 
of the properties we investigate are generally applicable to 
twisted magnetic flux tubes. For a review of recent theoretical 
and observational results on the behavior of magnetic flux 
ropes, see Russell, Priest, and Lee (1990). 

Coronal loops consist of magnetic flux tubes which can 
result either from the emergence of bundles of twisted field 
structures from beneath the photosphere, or can be formed 
from the ambient coronal magnetic field by the twisting 
motions inherent in subphotospheric convection. The large 
photospheric density and the small value of the plasma 
resistivity create a situation of almost “ frozen flow,” in which 
the footpoints of the coronal magnetic field lines which 
penetrate the photosphere are convected by the photospheric 
flow. It has long been recognized (Gold and Hoyle 1960) that 
coronal loops are likely to be unstable to kink modes if they 
become too twisted. It is well known from laboratory-plasma 
theory and experiments that pinch and tokamak discharges 
become unstable to ideal magnetohydrodynamic (MHD) kink 
modes if the toroidal plasma current (which is proportional to 
the twist of the field) exceeds a threshold amount (Kruskal ei 
at. 1958; Shafranov 1957). 

■ .V'. •• •. 


It was, however, recognized that on the Sun photospheric 
line tying provides a stabilizing influence which allows coronal 
loops to be stable up to higher levels of twist than 
corresponding laboratory plasmas (Raadu 1972; Giachetti, 
Van Hoven, and Chiuderi 1977; Einaudi and Van Hoven 
1981). Hood and Priest (1979, 1981) and Einaudi and Van 
Hoven (1983) have computed stability thresholds for the ideal 
MHD kink mode in coronal loops which explicitly 
demonstrate the stabilizing effect of photospheric line tying. 
Recently, Velli, Einaudi, and Hood (1990a) have analyzed the 
linear, ideal MHD stability of kink modes in coronal loops, 
and have presented a technique for the investigation of 
resistive instabilities (Velli, Einaudi, and Hood 19906). 

It is of interest to determine the dynamical consequences of 
kink instability in coronal loops, especially as they affect 
observable coronal phenomena. We describe a method which 
is employed in the investigation of the temporal evolution of 
twisted flux tubes. We advance the three-dimensional MHD 
equations numerically in time, using ideal and resistive models, 
to simulate the dynamics of the coronal plasma. Our model 
allows us to self-consistently follow the complete history of the 
evolution of flux tubes, including the formation, equilibrium, 
linear instability, and nonlinear behavior. Starting from an 
initial uniform background magnetic field, a twisted flux tube 
is created by the application of slow, localized photospheric 
vortex flows. For twists beyond a critical threshold, flux tubes 
become linearly unstable to ideal or resistive kink modes. The 
nonlinear evolution of these kink instabilities may generate 
concentrated current filaments. The resistive dissipation of 
these current filaments may provide a heating source for the 
corona. In this paper we discuss the equilibrium properties and 
the linear stability behavior. The nonlinear behavior will be 
treated in the second paper of this series. 


FILAMENT COOLING AND CONDENSATION 
IN A SHEARED MAGNETIC FIELD 
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Abstract Thermal instability driven by optically thin 
radiation in the corona ia b e ttere d to initiate the forma* 
tion of aoiar filaments (Parker. 1953). The fact that fila- 
menu are observed generally to sep ar a t e r e gion a of oppo- 
site. liae-of-sight, magnetic polarity in the diffe re nt ially 

that eiMWMWt wm. 

quirea the pre s e n ce of a highly aheared magnetic held. 
In thia paper we riiaraaa the coopled energ eti ca and dy- 
namics of the most important mniimaatimi modes, those 
due to p etp endint la r t herm al co n d uct ion at abort wave- 
lengths. We describe their linear structure in the aheared 
field and their growth rates. We have also performe d 
two-d u n ma ional. nonlinear, m a gnstnh y dr o dyu amic sim- 
ulations of the evolution of these modes in a foroe -f re s 


field. To clarify the rwasntial p h y sics of the — K -" be- 
havior, we have traced the evolution of generic pertuib a- 
tioas paeaseemg breed rpatial profiles. The shnulatiaas 
achieve the fine thermal s tructur es, «■■■»<■■■■■■» tempera- 
tures and maximum densities characteristic of o bs er ved 


onlae 


The anisotropic influence of shear on the 

linear dynamics of the condensation instability has been 
examine d extensively in a series of studies (Chhaderi and 
Van Hoven, 1979; Van Horen and Mok, 1984; Van Horen 
et aL, 1986; and Sparks and Van Horen, 1988]. Several 
conclusions have em er g ed from them undies; 

1. A condensation will form p — *— «*;«iiy ia regio ns 
in which kg s k-B/fi u 0 where k represents the (local) 
wave vector of the per turba tion. 

2. A condensation may be »>»— accordin g to 
whether iu spatial structure is determined primarily by 
considerations of force balance or e n er g y balance. The 
former is identified as a dynamic donation and the 
latter a hmgnatic 
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3. The modes that exhibit the fastest growth exist 
only in the presence of anisotropic heat flow. 

4. Plasma mass flow within kinematic 
is directed parallel to the magnetic field. 

3. Kinematic condensa t ions are most compressible 
when sound w a ve s traveling parallel to the magnetic field 


To study the evolution of solar filaments, we adopt 
a ana-fluid, transport model [Van Hoven et aL, 1987). 
The dynamical be h avior of the plasma is g ov er ned by the 
nonlin e ar equations of resistive magnetohydrodynaaucs, 
specifically including the continuity equation, the mo- 
mentum aquation. Maxwell's equations, and Ohm’s law. 
This set of equations is closed by including an energy 
equation of the farm 

* " 7 5d +h_1)[7 *T-c + si (i) 

where a is the thermal conductivity tensor. C = Rp t T’ 
depot ss due to radiative loss [Hildner. 1974). and 

H is a hsatiag function. We specify the magnetic Add 
by the model farm Be - B% (sech (y/a) e, + tanb (y/a) 
«,] w h ere a is the shear scale, which is ia dynamic equi- 
librium with a unifor m temperature and density, with an 
m o jpgjyifl g 4D0£8ftsc bftliDfic of uoifais 
radiation The MHO system can be red u ced to asmmha- 

tbe ta mp er utu re T, the three fluid- velocity componen ts 
Ur, tty, and u„ the magnetic field component B„ sad the 
magnetic flux function #(y,x). 

Tit uii nf*— ♦»*»** the “h-T-***” of « K * initial p lM »» exo- 
tations that lead to «c** i ** Mt ff oobng and condensation, 
— - — co ns ide r » modd rJ the ■man.MwpKtnde 

b eh a vior of this system, lb do so. we first linearise the 
MUD equations around the f orce-fre e , isobaric equilib- 
rium just described by uringT *» T 0 + Ti(y) exp (vt+ikx) 
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Abstract: We present results of time-dependent MHD 
simulations of mass ejections in the solar corona. Pre- 
vious authors have shown that results from simulations 
using a thermal driving mechanism are consistent with 
the observations only if an elaborate model of the initial 
corona is used. Our first simulation effort, using a sim- 
ple model of a plasmoid as the driving mechanism and a 
simple model of the initial corona, produces results that 
are also consistent with many observational features, sug- 
gesting that the nature of the driving mechanism plays 
an important role in determining the subsequent evolu- 
tion of mass ejections. Our first simulations are based on 
the assumption that mass ejections are driven by mag- 
netic forces: we are now developing simulations where the 
initial corona is perturbed magnetically by introducing 
a “plasmoid-like" current perturbation. The preliminary 
results from these simulations show some features that are 
consistent with the observations, others that are not. The 
discrepancies may be caused by the lack of internal force 
balance in the initial plasmoid structure. In the future, 
we plan to perform simulations where plasmoid formation 
occurs self-consistentlv. 

Introduction 

Coronal mass ejections (CMEs) have been observed 
in space-based white light coronagraphs since the early 
1970s. Because of the complex nature of CMEs, their ini- 
tiation and propagation in the solar corona is still not well 
understood. An important class of CME events are the 
looplike CMEs [Wagner, 1984]. Sime et al. [1984] iden- 
tified three general characteristics of the looplike CMEs 
observed with the Skvlab coronagraph: (1) the sides of 
the loop (the “legs" of the mass ejection) have a greater 
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density enhancement than the top of the loop, (2) a deple- 
tion of density occurs within the loop; (3) the legs exhibit 
very little lateral motion throughout the time evolution 
of the mass ejection. 

One approach for investigating CMEs is the time-dep- 
endent magnetohydrodynamic (MHD) simulation. The 
first MHD simulations of CMEs were based on the as- 
sumption that mass ejections are initiated by the rapid 
injection of thermal energy at the base of the corona from 
a solar flare [e.g.. Dryer et al., 1979; Wu et al., 1982], 
In these simulations, the initial corona was modeled hy- 
drostatically with a current-free magnetic field. A pres- 
sure pulse was introduced to model the thermal energy 
released in the flare, and the subsequent time evolution 
was followed. The resulting expanding fast-mode shock 
wave was identified as the CME. Although these simula- 
tions did give a coherent explanation for the initiation and 
propagation of CMEs. Sime et al. [1984] pointed out that 
these simulations failed to replicate the observed charac- 
teristics of looplike ejections. Because of these discrep- 
ancies, Sime et al. questioned the idea that CMEs could 
be regarded as compressional disturbances initiated by 
thermal energy release. 

In order to obtain simulation results that are compat- 
ible with the observations, Steinolfson and Hundhausen 
[1988] examined the role of the initial corona. They first 
verified that simulations with a thermal driving mecha- 
nism and an initially hydrostatic, current-free corona with 
a closed magnetic topology could not reproduce any of the 
observed features of mass ejections. They then showed 
that if one used a specially tailored initial condition and 
a thermal driving mechanism, many of the observed fea- 
tures of CMEs could be reproduced. Their model assumes 
that the evolution of CMEs in the solar corona is primar- 
ily influenced by the pre-event corona, and that the na- 
ture of the driving mechanism does not greatly affect the 
dynamics. 
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ABSTRACT 

Recent observations of intense, impulsive y-ray- and X-ray-emitting solar flares underline the suddenness of 
these events. The simultaneous emission of X-rays >40 keV from electron bremsstrahlung and y-rays 
requiring several MeV protons shows that all particles must be accelerated in less than 5 s. This paper pro- 
poses a simple model to explain such events, using the energy stored in the stretched held lines of a coronal 
arch. When reconnection occurs at the top of the arch, field lines retract like stretched rubber bands, sweeping 
up plasma and acting like a piston or slingshot. When the slug of plasma caught in the magnetic fields strikes 
the photosphere, it deposits its considerable kinetic energy, heating and compressing the intruding slug. Ten 
slugs of 100 km radius striking the photosphere may account for the 10” ergs radiation from loop flares. 
Subject headings: gamma rays: bursts — hydromagnetics — particle acceleration — Sun: flares — 

Sun: X-rays 


1. INTRODUCTION 

There is ample evidence that intense y-ray- and X-ray- 
emitting solar flares require sudden processes, capable of yield- 
ing X-rays of 40 keV or more and y-rays demanding protons of 
10 MeV energies. Recent observations underline these con- 
straints (Forrest & Chupp 1983; Nakajimaet al. 1983; Ramaty 
et al. 1980). Earlier work assumed that electrons were acceler- 
ated first in solar flares (Svestka 1977;. Brown & Smith 1980; 
Ramaty et al. 1980), yielding hard X-rays ~ 100 keV. Slow 
acceleration processes then gave second-phase energetic ions, 
and heat conduction down from the reconnection site heated 
the feet of the arch. 

Here we briefly propose that direct kinetic transport may be 
an important contributor to the overall energetics. If the relax- 
ing magnetic field lines in a loop can directly accelerate plasma 
(Spicer 1977; Van Hoven 1979). this will lead to energy transfer 
which is quicker and more efficient than several-phase models. 
It can yield simultaneous emission throughout the X-ray and 
y-ray spectrum from the same location at the arch feet. 

The basic physical picture is of a cold, beamlike slug, ejected 
from the reconnection region. The slug thermal speed is less 
than its directed velocity, which is of the order of the Alfven 
speed. This plasma strikes the photosphere, exciting ion-ion 
two-stream instabilities. These modes heat the slug to energies 
approaching a MeV, much faster than Coulomb collisions. 

This paper sketches the basic energetics argument. Colliding 
plasmas may well suffer turbulent stopping processes, and 
these we discuss briefly, using a marginal instability estimate. 
Future analysis will take up the detailed problems of deceler- 
ation and thick-target emission. 

2. MODEL 

Figure 1 shows a coronal loop which has suffered some 
reconnection between field lines near the top. This may occur 
through convective twisting, or when buoyant field lines of 
opposite sign rise into the underside of an existing loop, 
reversing toroidal fields locally, so that reconnection occurs. A 
new arch can also emerge inside a preexisting arch of opposite 
polarity, although here we discuss primarily a single loop. 
Some heat is liberated near the reconnection, and the notable 
effect far from the site is an outflow of both plasma and mag- 


netic fields at velocities near the Alfven speed, v A (Petschek 
1964; Mok & Van Hoven 1982; Forbes & Priest 1987). This 
work showed that reconnection cannot occur without reversal 
of local toroidal field lines. We reflect this in the reversed con- 
figuration of Figure 1. which we assume can apply to a variety 
of distance scales. In loops, reversals at the footpoints are of 
size ~ 100 km and not readily observable (Svestka 1981). 

While the immediate reconnection energy release is con- 
siderable, much more available energy is stored in the field 
lines extended from the photosphere below. These can now 
contract, shortening the local field-plasma tube (Hewaerts 
1981). 

This retracting tube will act like a slingshot, since magnetic 
forces dominate in the arch; the equilibrium ratio of plasma 
pressure to magnetic pressure, /?,. is 3.5 x 10' > n B TJB\. where 
the plasma density n 8 is in units of 10 8 cm' 3 . T 6 = (T 10'° K), 
and B 2 = (B 10' 2 G). These are typical values taken from 
observations (Svestka 1981). The vast magnetic field energy 
can convert to particle kinetic energy. 

The force acting on a plasma slug is j x B„, where j is the 
current perpendicular to the plane of Figure 1 and B„ is the 
magnetic field perpendicular to the main toroidal field of the 
figure. We shall assume this component is considerably less 
than the equilibrium field, which is of order 100 G. B„ forms 
the cup of the slingshot. 

What is the average ratio of thermal energy to magnetic 
energy in the slug, fj, as it leaves the reconnection region? We 
assume an essentially force-free loop configuration, so there 
will be a field B, perpendicular to the plane of Figure 1 and 
probably comparable in strength to the outside toroidal field. 
Reconnection will not annihilate this component, and it will 
resist compression of the slug plasma, setting a limit on heating 
from this compression. We shall take /? < 1. but not necessarily 
very small, and certainly ji > /!,, the equilibrium loop value. 

The compressibility of the slug plasma is governed by 

8 nnT + B} 

P * = IT- 

rather than simply fj. We expect 0, = 0.1-1. because B, < B for 
the equilibrium, but B, can be of order B. This implies that 
contraction of field lines can increase kinetic energies of 


rb?T 
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The application oi supercomputers and advanced numerical techniques to problems of coronal structure and dynamics is 
described. Numerical methods appropriate for the long time scale simulation of nonlinear magnetohydrodynarruc systems are 
discussed. Three specific examples of the application of these techniques to the solar corona are given. These are magnetic 
energy storage and conversion, a model for steady coronal heating, and calculation of stable force-free equilibria from given 
boundary data, such as that obtained with a vector magnetograph, h is suggested that the continued applicauon of these 
methods will result in substantial advances in the understanding of coronal dynamics and structure. 


1. Introduction 

The solar corona, or outer atmosphere of the 
sun. consists of a hot. tenuous plasma permeated 
with a magnetic field. This magnetized plasma 
displays a rich variety of dynamical activity. The 
coronal magnetic field is rigidly attached to the 
photosphere, the visible solar surface. The foot- 
pomts of these field lines are constantly being 
moved about slowly by the motion of the photo- 
sphere. which is driven from below bv convective 
processes. This motion produces a continual stir- 
ring of the coronal plasma. New magnetic flux 
emerges from below the photosphere and dis- 
places plasma above it. Sunspots and their associ- 
ated strong magnetic fields interact with each 
other, causing the magnetic field lines to slowly 
twist into loops and arcades resembling croquet 
hoops. Immense filaments of dense, cool plasma 
mvsienouslv form within the corona, and persist 
for periods of days or weeks. Occasionally these 
slow, churning processes yield to rapid, dramatic 
events. The local magnetic structure of the corona 
may be suddenly altered, with the release of vast 
amounts of energy in a very short period of time. 


A broad band of radiative energy is released in the 
form of a solar flare. A cool, dense filament may 
erupt from the low corona. Rapid ejections of 
mass into interstellar space may accompany topo- 
logical changes in the coronal magnetic field. These 
events may occur in concert or individually, and 
often seem to emerge spontaneously and without 
warning from the relatively quiescent corona. 

A quantitative theoretical desenpuon of the 
dynamical corona is one of the main challenges of 
modern solar physics. It is generally believed that 
the source of the energy released in disruptive 
events such as solar flares, erupting prominences, 
and coronal mass ejections is the coronal magnetic 
field. The energy in this field has two components. 
One is due to electric current flowing below the 
photosphere. This energy can generally not be 
tapped as this would require the alteration of 
distant current sources. The other arises from elec- 
tric currents flowing locally in the corona itself. 
This energy can be tapped by the local rearrange- 
ment of magnetic field. The origin of this free 
energy is likely the work done on the field by the 
slow motion of the photosphenc surface. The 
central question is how large amounts of magnetic 
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Abstract We have used a time-dependent magnetohy- 
drodynamic (MHD) simulation to investigate the mag- 
netic topology of an idealized streamerlike configuration 
in three dimensions. Starting from an initially current- 
free multipole field and a transonic flow, a quasi-steady, 
streamerlike configuration forms. The simulated streamer 
exhibits a closed magnetic field region bounded by a cur- 
rent layer and surrounded by open fields. This helmet re- 
gion is surmounted by a current layer, with an enhanced 
density in the closed field and current sheet regions. The 
simulated streamer has a finite longitudinal extent, with 
the closed field region primarily confined to the center of 
the structure. 

Introduction 

Coronal or helmet streamers are large-scale, relatively 
long-lived structures observed in the solar corona. They 
encompass spatial scales of a flew solar radii, and exist 
on time scales of a few days to as long as a few months. 
Coronal streamers appear in the corona (in white light) 
as structures that are brighter than the background (e.g. 
Poland, 1978; Priest, 1982, and references therein) and 
have been generally interpreted as outlining coronal mag- 
netic fields. As the white-fight in the corona near the 
sun is predominantly scattered by electrons in the coro- 
nal plasma, the larger relative brightness of these objects 
indicates that they are denser than the background. 

Streamers are often associated with coronal mass ejec- 
tions (CMEs), and observations have shown specific in- 
stances where a streamer was disrupted by a CME [e.g., 
Bildner et al., 1975; filing and Bundhansen, 1986). Cur- 
rent theoretical ideas suggest that CMEs may result from 
the driven instability of large-scale plasma and magne- 
tic field configurations [e.g.. Priest, 1988); therefore, the 
magnetic structure of coronal streamers may play a role 
in CME initiation. 

One of the first numerical studies of coronal streamers 
was by Pneumann and Kopp (1971). They used an itera- 
tive procedure to solve the steady-state MHD equations. 
Their solution stipulated that the plasma be hydrostatic 
in the closed field region, with flow parallel to the magne- 
tic field in the open field region. A cusplike neutral point 
was assumed to occur at the boundary between the open 
and closed field lines. Using this model, Pneumann and 
Kopp obtained a two-dimensional (axisymmetric) stream- 
erlike configuration. Their solution exhibited the general 
properties expected of a coronal streamer; a region of 
closed magnetic field, bounded by a current sheet and 
->pen magnetic fields. 
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Another approach that has been used to investigate 
the structure of coronal streamers is the time-dependent 
MHD simulation [e.g., Steinolfson et oL, 1982; Waskimi et 
a L, 1987). In this approach, one solves the the MHD equa- 
tions as an initial-value problem, starting from a given 
initial condition. The system is evolved in time until a 
steady-state (in practice only a quasi-steady) configura- 
tion is reached This approach has the advantage that 
the only assumption imposed is the initial condition, and 
assumptions about the properties of the solution, such as 
the nature of the currents, etc., are unnecessary. Like 
the Pneumann and Kopp calculation, these earlier MHD 
simulations of streamerlike configurations were also 2-D 
( axisymmetry assumed). Howe v er, observations of coro- 
nal streamers (e.g., Poland, 1978] indicate that they have 
a finite longitudinal extent, and are thus 3-D objects. No 
models have attempted to depict the 3-D magnetic topol- 
ogy of a streamer. In this paper, we use time-dependent 
MHD simulations to examine, for the first time, an ideal- 
ized, three-dimensional streamerlike configuration. 

Description of the Simulation 

Our simulation code solves the following form of the 
MHD equations (in normalized units) in spherical coordi- 
nates a two-step Lax-Wendroff method (similar to 
Linker et aL, 1990): 

g + V(pr) = 0 ( 1 ) 

(pv) + V • (pv v + PI) = J x B 

-Mr/r l -i-VW (2) 

| + 7-(P»>(l-l)(-P(V.v) 

+ ^ J -^W : Vu) (3) 

|-7x(«xl) = -Vx(,j) (4) 

In the above equations, W is a viscous stress tensor used 
for numerical stability purposes. Rt is the fluid Reynolds 
number (as 50 — 100), 7 is the resistivity, and 7 is the 
polytropic index. The simulation grid is 71 x 17 x 17 (in 
r, 8, d) with Ar = 0.1, A0 = 6°, and Ad — 6° (this in- 
cludes extra grid points used in the 8 and d directions for 
formulating the boundary conditions). When expressed 
in spherical coordinates. Equations (1-4) have a singular- 
ity at the poles of the coordinate system (e.g., 9 = 0 °). 
We calculate at the poles by deriving and differencing a 
form of the equations valid at the poles [ Linker , 1987). 


